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FOREWORD 


Acid  deposition  occurs  when  acidifying  pollutants  emitted  from  anthropogenic  and  other 
processes  undergo  chemical  reactions  in  the  atmosphere  and  fall  to  the  earth  as  wet 
deposition  (rain,  snow,  cloud,  fog)  or  dry  deposition  (dry  particles,  gas).  Acidic 
pollutants  can  be  transported  long  distances  in  the  atmosphere  from  their  sources  and 
eventually  be  deposited  in  ecosystems  over  broad  regional  scales  and  in  locations  far 
from  the  emission  sources. 

Dry  deposition  is  generally  more  a  local  problem  than  wet  deposition.  Direct 
measurement  of  dry  deposition  rates  is  difficult.  Dry  deposition  depends  on  many 
factors,  including:  meteorological  conditions,  characteristics  of  the  pollutants  being 
deposited  (e.g.  different  gaseous  chemical  and  particle  size),  and  characteristics  of  the 
surface  on  which  deposition  occurs. 

The  most  accepted  and  common  method  for  estimating  dry  deposition  is  the  so-called 
"inference  method."  The  inferential  method  is  a  combination  of  measurement  and 
modeling  that  involves  indirect  estimation  of  dry  deposition  rates  on  the  basis  of 
routinely  measured  air  concentrations  and  meteorological  parameters.  The  method  is 
based  on  an  assumed  steady-state  relationship  F  =  Vd  C,  where  the  dry  deposition  flux  or 
rate  (F)  is  a  product  of  the  dry  deposition  velocity  (Vd)  and  the  concentration  (C)  of  an 
airborne  pollutant. 

A  series  of  studies  have  been  initiated  by  AENV  to  evaluate  the  inference  method  and 
search  for  the  most  suitable  and  simple  model  for  deposition  rate  estimations  in  Alberta. 
This  report  documents  the  third  study  in  the  series.  Titles  for  the  reports  of  the  other 
studies  are:  ''Review  and  Assessment  of  Methods  for  Monitoring  and  Estimating  Dry 
Deposition  in  Alberta  ",  and  ''Dry  Deposition  Monitoring  Method  in  Alberta  ".  It  is 
anticipated  that  once  all  necessary  information  is  gathered,  an  Alberta  protocol  for  dry 
deposition  measurement  will  be  prepared. 


Lawrence  Cheng,  Ph.  D. 
Air  Policy, 

Climate  Change,  Air  and  Land  Policy  Branch 
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SUMMARY 


Environment  Canada  developed  a  dry  deposition  inference  model  to  predict  deposition  of 
gaseous  and  particulate  species  onto  land  or  water  in  Canada.  Alberta  Environment  has 
developed  a  dry  deposition  inference  model  to  similarly  predict  deposition  of  gaseous  and 
particulate  species  onto  land  and  water  in  Alberta.  While  numerous  dry  deposition 
inference  models  are  used  in  practice  in  North  America,  uncertainty  exists  in  results  from 
these  models  because  of  an  inability  to  take  direct  measurements  to  test  their  validity. 

The  objectives  of  this  study  were  to  examine  the  performance  of  the  Alberta  Environment 
inference  model  for  predicting  dry  deposition  in  Alberta;  and  to  identify  and  recommend 
changes  to  a  number  of  assumptions  and  boundary  conditions  in  the  model  to  improve  its 
performance  relative  to  the  Environment  Canada  inference  model.  The  following 
findings  and  recommendations  are  noted: 

1.  The  three  most  important  species  contributing  to  potential  acid  input  among 
eleven  species  considered  in  the  Alberta  Environment  model  were  SO2,  NO2,  and 
HNO3. 

2.  Sensitivity  testing  of  resistance  factors  (Ra,  Rb,  and  Rc)  used  in  Alberta 
Environment's  model  showed  that  step  changes  in  the  magnitude  of  Rc  are  much 
more  important  in  influencing  total  deposition  compared  to  step  changes  in  Ra  or 
Rb  for  all  depositing  species  examined.  This  is  generally  consistent  with  scientific 
literature. 

3.  AENV  model  performance  in  terms  of  overall  deposition  flux  of  1 1  species  was 
observed  to  be  closely  dependent  on  deposition  velocity  modeling  for  SO2  and 
NO2;  particularly  with  respect  to  selection  of  Rc  values  used  to  parameterize  these 
gases.  Relative  to  performance  of  the  ENVC  model,  values  used  to  represent  Rc 
in  the  AENV  model  were  less  conservative  for  SO2  in  winter  and  spring  and  more 
conservative  in  summer  and  autumn.  For  NO2,  values  used  to  represent  Rc  in  the 
AENV  model  were  less  conservative  in  autumn  and  winter  and  more  conservative 
in  spring  and  summer.  For  HNO3,  deposition  in  the  AENV  model  was  driven  by 
Rb.  These  differences  do  not  necessarily  imply  problems  in  the  AENV  model. 
The  ENVC  inference  model  does  not  necessarily  represent  the  best  set  of 
assumptions  and  conditions  to  parameterize  dry  deposition.  In  addition,  AENV's 
parameterization  of  Rc  originates  from  previously  published  values  that  are 
seasonal-  and  species-specific  and,  therefore,  can  be  considered  reasonable. 

4.  Ra  is  an  indication  of  the  mixing  conditions  in  the  atmosphere  such  that  a  very 
low  Ra  implies  maximum  mixing  conditions  (experienced  during  daytime). 
Setting  a  minimum  value  of  5  s/m  for  Ra  in  the  Alberta  Environment  model  had 
only  a  small  positive  change  in  model  performance  relative  to  the  Environment 
Canada  model.  Although  the  effect  of  this  change  is  small,  it  is  recommended  for 
future  use  because  a  minimum  limit  for  Ra  acknowledges  that  there  is  always 
some  limit  to  mixing  possible  in  an  atmosphere. 
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5.  Setting  a  maximum  value  of  1000  s/m  for  Ra  in  the  Alberta  Environment  model 
resulted  in  better  agreement  in  overall  deposition  flux  of  species  relative  to  the 
Environment  Canada  model  for  the  conditions  tested.  Very  high  R.,  values  would 
be  associated  with  no  mixing  in  the  atmosphere  (ideal  stable  night  time 
conditions).  An  extremely  high  Ra  value  would  not  necessarily  offer  a  reasonable 
representation  of  natural  conditions.  From  a  practical  point  of  view,  large  or 
moderately  large  Ra  values  (and  correspondingly  small  or  moderately  small 
deposition  velocity  and  deposition  fluxes)  are  unimportant.  This  upper  limit  is 
recommended  for  future  use  in  the  Alberta  Environment  model. 

6.  The  Monin-Obukhov  length  is  used  in  calculation  of  the  stability  function  which 
is  ultimately  used  to  calculate  Ra.  Maximum  and  minimum  values  of  -5  and  5  m 
for  negative  and  positive  limits  are  recommended  for  future  use  in  the  Alberta 
Environment  model  as  small  absolute  Monin-Obukhov  length  values  correspond 
to  extremes  which  are  not  representative  of  natural  conditions. 

7.  A  minimum  wind  speed  value  of  1  m/s  is  recommended  for  future  use  in  the 
Alberta  Environment  model. 

8.  As  part  of  the  calculation  steps  for  Ra,  minimum  requirements  discussed  in  more 
detail  in  this  report  are  recommended  on  how  to  calculate  surface  roughness 
length  (Zo)  as  a  monthly  average  using  hourly  data  where  the  wind  speed  is  >6  m/s 
in  the  Alberta  Environment  model. 

9.  Wetness  alters  the  characteristics  of  a  surface  to  act  a  sink  for  depositing  species 
(primarily  gases).  The  Alberta  Environment  model  requires  an  hourly 
measurement  to  represent  the  amount  of  wetness  on  a  surface  in  determination  of 
surface  resistance  (Rc).  The  model  is  set  up  to  use  field  measurement  data  from 
surface  wetness  sensing  equipment.  However,  past  experience  has  indicated  that 
this  field  equipment  does  not  provide  reliable  measurement  data  in  an  Alberta 
setting  due  to  the  wide  range  of  climatic  conditions  experienced.  Alternatively, 
the  model  allows  a  relative  humidity  measurement  cutoff  (87%)  to  be  used  as  an 
approximation  to  represent  a  surface  as  being  wet.  It  is  recommended  that  this 
criterion  be  maintained  in  the  model  for  future  use. 
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1.0  INTRODUCTION 


1.1  Background 

Environment  Canada  (ENVC)  developed  a  dry  deposition  inference  model  to  predict 
deposition  of  gaseous  and  particulate  species  onto  land  or  water  in  Canada  (Zhang  et  al., 
2003a,  2003b,  2002,  2001;  Brook  et  al.,  1999).  Environment  Canada^ s  model  uses 
pollutant  and  meteorological  data  collected  from  the  Canadian  Air  and  Precipitation 
Monitoring  Network  (CAPMoN)  (Brook  et  al,  1999).  Alberta  Environment  (AENV)  has 
developed  a  dry  deposition  inference  model  (Cheng  et  al.,  2001)  to  similarly  predict 
deposition  of  gaseous  and  particulate  species  onto  land  and  water  in  Alberta.  Both 
approaches  are  a  form  of  dry  deposition  inference  modeling. 

Dry  deposition  inferential  models  evolved  from  methods  used  by  meteorologists  and 
foresters  to  characterize  H2O,  CO2,  heat,  and  momentum  exchange  between  the 
atmosphere  and  plant  canopies.  While  numerous  dry  deposition  inference  models  are 
used  in  practice  in  North  America  and  elsewhere,  uncertainty  exists  in  results  from  these 
models  because  of  an  inability  to  take  direct  measurements  to  test  their  validity  (Lovett, 
1994).  Dry  deposition  inference  models  should  have  the  capability  to  represent  dry 
deposition  of  atmospheric  pollutants  in  a  manner  that  reflects  current  scientific 
understanding.  Also  important,  results  from  any  one  approach  should  compare 
reasonably  well  with  other  approaches  used  to  represent  dry  deposition.  It  has  however 
been  noted  by  Wesely  and  Hicks  (2000)  that  further  work  on  parameterizing  inferential 
models  is  necessary  before  resultant  dry  deposition  estimates  can  be  considered  reliable 
and  this  needs  to  be  considered  when  interpreting  and  applying  model  outputs. 

Alberta  Environment  is  interested  in  developing  their  dry  deposition  inference  model  for 
wider  use  in  Alberta.  WBK  &  Associates  Inc.  (WBK,  2008)  undertook  initial  testing  of 
the  ENVC  and  AENV  models  using  a  common  set  of  meteorological  and  pollutant  data 
to  examine  how  they  compared  for  predicting  dry  deposition  of  eleven  gaseous  and 
particulate  species.  Initial  findings  of  this  work  indicated  both  models  gave  comparable 
results  for  amounts  of  dry  deposition  for  gaseous  and  inorganic  species  on  a  monthly  and 
annual  basis.  One  of  the  recommendations  from  this  study  was  to  undertake  further 
comparative  testing  of  these  two  models  to  better  understand  reasons  why  model  results 
departed  from  each  other. 

1.2  Objectives 

Further  investigation  of  the  AENV  and  ENVC  dry  deposition  inference  models  was 
undertaken  as  a  follow-up  to  recommendations  of  WBK  (2008).  The  objectives  of  this 
study  were:  i)  to  examine  the  performance  of  the  Alberta  Environment  inference  model 
relative  to  the  ENVC  model  for  predicting  dry  deposition  in  Alberta;  and  ii)  to  use  this 
performance  comparison  to  recommend  changes  for  a  number  of  assumptions  and 
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boundary  conditions  in  the  AENV  inference  model  to  improve  its  performance.  The 
ENVC  inference  model  does  not  necessarily  represent  the  best  set  of  assumptions  and 
conditions  to  characterize  dry  deposition.  However,  results  from  any  one  approach 
should  compare  reasonably  well  with  other  approaches  used  to  characterize  dry 
deposition. 

1.3  Report  Organization 

Dry  deposition  is  initially  discussed  and  general  modeling  approaches  of  acid  deposition 
are  presented  based  upon  current  scientific  literature.  The  ENVC  modeling  approach  is 
presented  in  detail  along  with  the  AENV  model.  The  discussion  is  split  into  gaseous  and 
particulate  deposition. 

In  general,  the  inference  method  involves  indirect  estimation  of  dry  deposition  rates  on 
the  basis  of  routinely  measured  air  concentrations  and  meteorological  parameters  (WBK, 
2006).  The  method  is  based  on  an  assumed  steady-state  relationship  F  =  Vd  C,  where  dry 
deposition  flux  or  rate  (F)  is  a  product  of  the  dry  deposition  velocity  (Vd)  and  the 
concentration  (C)  of  an  airborne  pollutant.  Vd  is  estimated  on  the  basis  of  resistance 
models  and  is  defined  as  -  using  an  Ohm's  law  analogy  -  the  inverse  of  the  sum  of 
multiple  resistance  factors  (aerodynamic  resistance  (Ra),  boundary-layer  resistance  (Rb), 
and  surface  resistance  (Rc))  (Wesely  and  Hicks,  2000,  1977).  Gaseous  deposition  is 
discussed  according  to  the  three  major  resistance  pathways  and  particulate  deposition  is 
discussed  according  to  important  acidic  species  (sulfate,  nitrate,  and  ammonium)  and 
base  cations  (calcium,  magnesium,  potassium  and  sodium). 

The  AENV  model  was  run  using  different  boundary  conditions  for  selected  model 
parameters  and  the  results  were  compared  with  the  model  run  under  original  conditions 
for  similar  inputs.  A  number  of  the  different  boundary  conditions  were  based  on 
conditions  set  in  the  ENVC  model.  A  summary  of  conditions  evaluated  is  presented  in 
Table  1 .  Finally,  the  AENV  model  was  run  using  different  boundary  conditions  and 
compared  with  the  ENVC  model  run  using  similar  inputs.  Based  upon  comparisons  of 
results  for  these  sets  of  model  runs,  recommendations  were  made  about  how  best  to 
proceed  with  using  assumptions  and  boundary  conditions  for  operating  the  AENV  model 
in  the  future. 
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Table  1         Boundary  conditions  evaluated  during  study. 


Parameter 

Units 

Original  Default 

Revised  DefauU 
Evaluated 

Basis  for 
Revision 

Ra 

s/m 

upper  limit 
undefined 

upper  limit  capped  at 
1000  s/m  (at  2000  s/m 
for  water) 

limit  set  in 
ENVC  model 

Ra 

s/m 

lower  limit 

lower  limit  capped  at  5 

limit  set  in 

undefined 

s/m 

ENVC  model 

Monin- 

m 

undefined  limit 

upper  limit  of  -5  m  on 

limit  set  in 

Obukhov  length 

negative  side  and 
lower  limit  of  +5  m  on 
positive  side 

ENVC  model 

wind  speed 

m/s 

set  to  0  m/s  if  less 

set  to  1  m/s  if  less  than 

limit  set  in 

than  0.09  m/s 

1  m/s 

ENVC  model 
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2.0 


ACID  DEPOSITION  BACKGROUND 


2.1  Introduction  to  Dry  Deposition 

Acid  rain  may  contribute  to  acidification  of  terrestrial  and  aquatic  systems  (AENV, 
1999).  A  portion  of  the  deposited  nitrogen  may  also  contribute  to  terrestrial  fertilization 
and/or  aquatic  eutrophication.  Acid  rain  involves  three  processes  including:  transport, 
transformation,  and  deposition  (Wesely  and  Hicks,  1977).  Acid  deposition  can  occur  by 
wet  and  dry  processes.  Dry  deposition  is  reported  to  represent  about  40%  of  total 
deposition  (Smith  et  al.,  2000). 

Increased  industrialization  in  Alberta,  particularly  oil  sands  extraction  and  processing  in 
north  eastern  Alberta  and  petroleum  upgrading  and  refining  activities  in  central  Alberta, 
is  raising  awareness  about  potential  acid  deposition  due  to  industrial  activities.  A  critical 
load  method  is  used  in  Alberta  as  a  receptor-based  approach  for  managing  acid 
deposition  (AENV,  1999).  Since  it  is  very  difficult  to  direcdy  measure  dry  deposition, 
models  are  essential  for  the  estimation  of  the  amount  of  dry  deposition  that  a  particular 
land  type  will  experience  given  an  atmospheric  concentration  of  a  depositing  species. 

Dry  depositing  species  can  be  divided  into  sulfur  compounds,  nitrogen  compounds,  and 
base  cations.  Sulfur  compounds  of  interest  include  a  gaseous  sulfur  dioxide  (SO2)  and 
particulate  sulfate  (S04^").  Nitrogen  compounds  include  several  gases  -  nitrogen  dioxide 
(NO2),  nitric  acid  (HNO3),  nitrous  acid  (HNO2),  nitric  oxide  (NO),  and  ammonia  (NH3). 
Particulate  forms  of  nitrogen  include  nitrate  (NO3")  and  ammonium  (NH4^).  Base  cations 
include  particulate  calcium  (Ca"^),  sodium  (Na^),  potassium  (K^),  and  magnesium 
(Mg'^-^). 

2.2  Inferential  Methods 

The  inferential  method  refers  to  a  method  of  approximating  deposition  processes  with  a 
resistance  analogy  (Hicks  et  al.,  1987).  The  inferential  method  reflects  the  concept  that 
different  surfaces  will  experience  different  deposition.  Early  inference  models 
approximated  the  entire  canopy  as  one  large  surface  and  were  referred  to  as  a  "big-leaf 
method  (Hicks  et  al.,  1987).  A  two-big-leaf  model  evolved  to  emphasize  separate 
treatment  of  sunlit  and  shaded  portions  of  a  canopy  (Zhang  et  al.,  2001a,  2001b). 

Multi-layer  models  (MLMs)  or  second  generation  models  have  evolved  that  include 
integration  of  deposition  along  the  canopy  height  of  vegetation  (Meyers  et  al.,  1998). 
Peters  et  al.  (1995)  reported  that  third  generation  models  will  include  interactions  among 
chemistry,  deposition  processes,  and  meteorology;  and  may  include  physiological 
phenomena  such  as  photosynthesis  and  CO2  production  (Niyogi  et  al.,  2003). 
Consequently,  third  generation  models  will  be  highly  complex.  A  schematic  of  possible 
interrelations  showing  this  complexity  is  included  in  Figure  1 . 
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Figure  1        Possible  interactions  for  potential  inclusion  in  third-generation  dry 
deposition  inference  models  (Peters  et  ah,  1995). 


2.3  Deposition  Velocity 

Deposition  velocity  is  a  parameter  multiplied  by  the  concentration  of  a  chemical  in  air  in 
order  to  calculate  chemical  flux  or  deposition  rate  (Wesely  and  Hicks,  2000).  Values  for 
deposition  velocities  are  reported  to  range  from  0  cm/s  for  inert  substances  to  about 
greater  than  2  cm/s  for  reactive  substances  such  as  nitric  acid  (HNO3)  (Wesely  and  Hicks, 
2000).  Deposition  velocities  are  quoted  in  several  units  including  m/s,  mm/s,  or  cm/s,  the 
latter  being  the  most  commonly  used.  The  behavior  of  gaseous  ammonia  near  a  soil 
surface  is  bidirectional  as  both  deposition  to  and  emission  from  soil  occurs,  depending 
upon  ambient  concentrations.  The  concept  of  deposition  velocity  should  not  be  used  for 
ammonia  at  very  low  ambient  concentrations  because  it  may  originate  from  soil  and  it 
violates  an  assumption  that  the  flux  is  unidirectional  (towards  the  surface)  (Hicks  et  al., 
1987;  Wesely  and  Hicks,  2000). 


Refinement  Study  of  Dry  Deposition  Inference 
IVIethods  Used  in  Alberta 


5 


2,3,1        Gaseous  Deposition 


Deposition  velocity  is  commonly  modeled  as  three  resistances  in  series  (Hicks  et  al., 
1987).  These  resistances  include:  aerodynamic  resistance  (Ra)  associated  with  turbulent 
transfer,  boundary  layer  resistance  across  the  near-surface  layer  (Rb),  and  surface 
resistance  involving  the  pollutant-surface  interaction  (Rc)  (Voldner  et  al.,  1986). 
Parameterization  methods  for  representing  these  resistances  in  the  Environment  Canada 
and  Alberta  Environment  inference  methods  are  presented  in  detail  in  Appendix  A  and  B, 
respectively. 

Aerodynamic  Resistance  (Ra) 

Parameters  used  for  modeling  aerodynamic  resistance  are  meteorological-specific, 
including:  wind  speed,  wind  direction,  and  air  temperature  at  two  heights  above  the 
ground  (usually  2  and  10  m).  Modeled  Ra  values  vary  depending  upon  stability  and 
surface  roughness  conditions,  which  are  determined  using  these  parameter  inputs.  Most 
inference  methods  represent  surface  roughness  by  using  default  values  for  different  sets 
of  surface  roughness  conditions  (e.g.,  Voldner  et  al.,  1986).  Because  the  input 
parameters  are  meteorological-  and  not  chemical-specific,  the  aerodynamic  resistance 
term  is  independent  of  chemical  species. 

Boundary  Layer  Resistance  (Rb) 

The  boundary  layer  is  a  thin  layer  next  to  the  surface  (refer  to  Figure  2).  Molecular 
diffusion  determines  gas  transfer  while  Brownian  diffusion  and  inertial  impact  determine 
particle  transfer  (WBK,  2008).  Parameters  needed  for  this  resistance  term  include:  wind 
speed,  roughness  length,  and  temperature.  Rb  is  usually  smaller  than  aerodynamic 
resistance;  which  means  it  is  usually  less  very  important  for  deposition  velocity 
calculations  (Walcek  et  al,  1986). 

Total  Surface  Resistance  (Rc) 

Rc  includes  all  processes  involving  the  actual  surface  and  the  particular  chemical  of 
interest  (Hicks  et  al.,  1987).  Most  of  the  difference  between  inference  models  is  in  the 
modeling  of  surface  resistance.  As  models  evolve  to  include  more  and  more 
physiological  phenomena,  variables  such  as  soil  pH,  surface  wetness,  and  photo 
synthetically  active  radiation  are  required  (Wesely  and  Hicks,  2000).  These  variables  can 
be  difficult  or  expensive  to  obtain  (measure). 

Sophisticated  parameterizations  of  Rc  have  represented  it  as  a  series  of  multiple 
resistances  and  leaf  area  index  of  the  canopy,  for  example  (Zhang  et  al.,  2003):  stomatal 
resistance,  mesophyll  resistance,  in-canopy  aerodynamic  resistance,  ground  resistance, 
and  canopy  cuticle  resistance.  Stomatal  resistance  involves  a  gas  entering  the  leaf 
through  the  stomata.  Mesophyll  resistance  involves  resistance  to  a  gas  dissolving  in  the 
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mesophyll  cells  in  the  interior  of  the  leaf.  In-canopy  aerodynamic  resistance  refers  to 
meteorological  issues  within  the  canopy,  including  temperature  gradients.  Ground 
resistance  includes  resistance  of  soil,  ground  litter,  and  any  water  surface  under  the 
canopy. 


Rc 

Figure  2        Relative  locations  where  dry  deposition  resistance  factors  Ra,  Rb,  and 
Rc  apply. 

For  a  water  surface  or  even  water  on  the  surface  of  soil  or  vegetation,  Rc 
parameterizations  can  include  aspects  of  the  nature  of  the  water  itself  in  terms  of 
parameters  such  as:  pH,  concentration  of  the  chemical  already  dissolved,  and  surface 
characteristics  of  the  water  (Zhang  et  al.,  2003).  Canopy  cuticle  resistance  includes 
interaction  with  the  cuticular  surface  of  the  plant.  Alternatively,  simpler 
parameterizations  of  Rc  have  represented  it  as  a  series  of  default  values  for  different 
surface  conditions  represented  by  season,  time  of  day,  and  surface  wetness  condition 
(Voldneret  al.,  1986). 

2J,2  Particles 

Particles  in  the  atmosphere  containing  chemical  species  or  elements  tend  to  have  a  range 
of  sizes  that  can  be  represented  as  a  log-normal  distribution  (Ruijgrok  et  al.,  1995). 
Larger  mass  median  elements  tend  to  have  a  unimodal  distribution  while  smaller  mass 
median  elements  (<2  |im  in  diameter)  tend  to  have  a  bimodal  distribution  (Milford  and 
Davidson,  1985).  The  general  size  distribution  of  a  particle  typically  depends  upon  its 
origin.  Particles  can  be  divided  into  three  size  modes  (Holsen  and  Noll,  1992): 

•  The  largest  mode  (>2.5|am)  results  from  mechanical  abrasion  and  wind  erosion. 

•  Fine  mode  particles  (0.1  to  2.5  |jm)  originate  through  coagulation  processes  and 
are  primary  aerosol  emissions. 

•  The  smallest  particles  (<0.1  |jm)  are  secondary  products  of  emissions. 
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Gases  may  also  react  with  larger  particles  and  be  present  in  larger  particles  (Ruijgrok  et 
al.,  1995).  This  duality  of  origin  may  explain  the  bimodal  shape  of  the  size  distribution. 
Atmospheric  chemical  reactions  involving  gaseous  precursors  to  particulates  can  be 
affected  by  presence  of  other  chemicals  (Ruijgrok  et  al.,  1995).  Nitrate  may  reside  on 
large  particles  as  a  result  of  reactions  of  nitric  oxide  with  alkaline  species  of  crustal  origin 
(Ruijgrok  et  al.,  1995;  Wolff,  1984).  If  ammonia  gas  is  present,  which  is  common  in 
animal  feedlot  areas,  nitric  oxide  may  react  with  it,  resulting  in  nitrate  (and  ammonium) 
on  small  particles  (Ruijgrok  et  al.,  1995;  Wolff,  1984). 

When  modeling  particulate  deposition,  it  is  important  that  the  size  distributions  are 
reasonable  representations  of  conditions  for  an  area  being  modeled.  Neighboring  salt 
water  bodies  (i.e.,  near  coastal  areas),  animal  feedlots,  and  soil  types  can  all  strongly 
influence  size  distribution  of  particulate  matter.  Thus  particle  size  distributions  may  be 
different  for  different  areas  or  different  times  of  the  year  (Ruijgrok  et  al.,  1997;  Wiman 
and  Agren,  1985;  Whitby,  1978;  Milford  and  Davidson,  1985;  Hoff  et  al.,  1983). 
Arimoto  et  al.  (1997)  stated  that  mass  particle  size  distributions  may  never  reach 
equilibrium.  In  addition,  Hedin  et  al.  (1994)  stated  that  the  amount  of  atmospheric 
cations  may  be  declining  over  time. 

Models  that  represent  deposition  based  on  particle  size  are  considered  process-oriented 
models;  as  opposed  to  bulk  resistance-oriented  models  that  do  not  distinguish  for  particle 
size  (Ruijgrok  et  al.,  1995).  Particle  deposition  mechanisms  include  (Slinn,  1982): 
diffusion,  interception,  impaction,  gravitational  settling,  and  phoretic  mechanisms. 
Brownian  diffusion  dominates  small  particle  (<0.1  pm)  deposition;  whereas  gravitational 
sedimentation  dominates  large  particle  (>2.5  |am)  deposition.  Figure  3  illustrates  that 
deposition  velocities  are  highest  for  large  particles  and  are  controlled  by  sedimentation. 
Deposition  velocities  reach  a  minimum  for  particles  from  0.1  to  1  pm  in  size. 

Process-oriented  models  can  be  further  increased  in  complexity  by  inclusion  of  phoretic 
effects;  including  electrophoresis,  thermophoresis,  and  diffusiophoresis  (Ruijgrok  et  al., 
1995).  Hygroscopic  growth  of  particles  can  be  an  issue  for  high  relative  humidities 
(>80%)  (Draaijers  et  al.,  1997;  Ruijgrok  et  al.,  1995).  For  example,  an  inference  model 
used  by  Erisman  et  al.  (1997)  includes  cases  for  wet  and  dry  surface  conditions  as  well  as 
for  relative  humidity  above  and  below  80%.  Another  complication  not  commonly 
addressed  in  particle  deposition  models  is  that  Brownian  diffusion  will  decrease  if  a  gas 
becomes  attached  to  a  particle  (Sehmel,  1980).  Finally,  atmospheric  stability  issues  may 
be  a  concern  if  a  model  was  only  designed  for  neutral  stability  conditions  (Ruijgrok  et  al., 
1995). 

High  deposition  velocities  related  to  sedimentation  of  large  particles  result  in  models  that 
can  be  extremely  sensitive  to  the  method  used  to  represent  the  large  size  end  of  the  size 
distribution  (Holsen  and  Noll,  1992).  Solely  using  the  mass  median  diameter  (MMD)  to 
represent  particle  size  without  integrating  across  the  whole  size  distribution  can  cause 
large  differences  in  model  results  (Ruijgrok  et  al.,  1995).  Ideally,  it  is  very  important  to 
represent  the  distribution  as  divided  into  sufficiently  small  portions  so  that  the  full  effects 
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of  the  large  particle  end  of  a  distribution  can  be  properly  represented.  For  example,  when 
using  1,  6,  or  100  different  diameters  to  represent  the  distribution  for  sodium-bound 
particulates  (mean  MMD  of  6.5  pm);  Dulac  et  al.  (1987)  calculated  mean  dry  deposition 
fluxes  as  915,  700,  and  30,000  pg/m7day,  respectively. 


PARTICLE  DIAKTER  (|im) 

Figure  3        Deposition  velocity  of  particles  as  function  of  size  in  a  stable 

atmosphere  for  two  different  roughness  lengths  (Zo,  cm),  particle 
densities  (p,  kg/m^),  and  horizontal  wind  speed  of  5  m/s  using  model  of 
Sehmel  and  Hodgson,  1980  (as  cited  in  Ruijgrok  et  al.,  1995). 


Large  particles,  while  rarer,  can  be  responsible  for  more  than  90%  of  the  flux  of  a 
chemical  because  of  large  deposition  velocities  (Holsen  and  Noll,  1992).  For  large 
particles,  process-oriented  models  using  numerous  size  distributions  are  highly  sensitive 
to  selection  of  the  largest  particle  size  as  well  as  to  representativeness  of  the  lognormal 
distribution  (Holsen  and  Noll,  1992).  Particles  originating  from  crustal  material  tend  to 
be  large  (>60  to  100  pm)  (Ruijgrok  et  al.,  1995.)  Given  that  larger  particles  mostly 
contain  base  cations  which  act  to  reduce  the  impact  of  acidic  components, 
underestimating  the  amount  of  deposition  of  these  larger  particles  may  act  as  a  safety 
factor  in  estimation  of  net  total  acid  deposition  or  potential  acid  input  (PAI). 
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2.4 


Dry  Deposition  Model  Issues 


Dry  deposition  can  be  difficult  to  model  because  of  scientific  complexities  and  the 
limitations  in  approaches  used  to  represent  these  complexities.  Others  have  reported  that 
it  is  not  possible  to  incorporate  all  relevant  physiological  and  atmospheric  processes  into 
an  acid  deposition  model  for  several  reasons  (Lovett,  1994;  Venkatram  et  al.,  1988): 

•  Incomplete  understanding  of  the  processes  involved. 

•  Limitations  in  computational  capacity. 

•  Limitations  in  model  validation,  including  limitations  of  measuring  equipment. 

Wesely  and  Hicks  (2000)  state  that  there  is  an  extremely  large  number  of  processes 
involved  in  the  transfer  of  a  chemical  from  air  to  a  surface.  Phenomena  that  are  currently 
difficult  to  represent  in  a  deposition  model  include  (Wesely  and  Hicks,  2000;  Lovett, 
1994;  Venkatram  et  al.,  1988):  surface  wetness,  surface  roughness,  edge  effects,  in- 
canopy  chemical  reactions,  and  bi-directional  chemical  gradients. 

Computational  capacity  is  always  an  issue  in  modeling  and  models  are  usually  designed 
up  to  the  limits  of  the  applicable  computer  processing  speed  and  capacity  (Venkatram  et 
al.,  1988).  Modeling  natural  processes  always  involves  approximating  complicated 
relationships,  including:  interactions,  thresholds,  and  feedback  loops.  In  this  study,  the 
requirement  to  have  a  model  that  is  relatively  easy  to  use  was  an  important  priority  so 
that  others  could  readily  employ  it  to  estimate  acid  deposition  behavior  in  Alberta. 

Model  validation  is  of  particular  concern  in  acid  deposition.  Current  measurement 
equipment  is  not  able  to  detect  small  concentrations  rapidly  enough  to  detect  or  to 
measure  deposition  phenomena.  Eddy  correlation  techniques  are  reported  to  be  the  most 
sophisticated  approaches  for  measuring  dry  deposition  (Meyers  et  al.,  1996;  Lovett, 
1994),  however,  this  method  requires  sensitive  anemometry  and  rapid-response  pollutant 
monitoring  equipment  (Lovett,  1994).  In  addition,  errors  in  measurement  techniques  are 
reported  to  frequently  be  higher  than  40%  (Businger,  1986).  Finally,  there  is  uncertainty 
as  to  whether  deposition  information  gathered  at  one  location  in  one  time  period  can  be 
applied  to  represent  different  locations  or  at  different  times. 

Wesely  and  Hicks  (2000)  reported  that  quasi-comparisons  of  model  results  with 
deposition  measurements  have  shown  differences  as  large  as  +/-30%.  These  comparisons 
refer  only  to  converting  measured  ambient  concentrations  into  deposition  estimates  and 
not  the  element  of  computer  dispersion  modeling  that  translates  emission  estimates  into 
ambient  concentrations;  which  introduces  another  source  of  uncertainty  into  modeled 
deposition  estimates. 

In  addition,  Wesely  and  Hicks  (2000)  state  that  the  Ohm's  law  analogy  for  representing 
deposition  velocity  is  likely  an  oversimplification  and  that  the  role  of  plant  physiological 
processes  are  not  sufficiently  incorporated  into  these  models.  For  modeling  gaseous 
deposition,  possible  modeling  limitations  can  be  due  to  a  failure  to  adequately  represent 
edge  effects,  heterogeneity  of  the  landscape,  and  terrain  complexity  (Hicks  et  al.,  1987). 
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Difficulties  associated  with  modeling  particle  deposition  mechanisms  may  result  in 
underestimating  particle  deposition  by  as  much  as  an  order  of  magnitude  (Holscn  and 
Noll,  1992). 
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3.0 


MODEL  DESCRIPTIONS 


3.1 


Environment  Canada  (ENVC)  Model 


3.1.1  Gases 

The  gas  module  of  the  ENVC  model  calculates  deposition  velocity  (V^)  as  the  inverse  of 
the  sum  of  three  resistances       R^,  and  Re)  as  follows  (Zhang  et  al.  2003): 

— —  (1 

Aerodynamic  resistance  {Ra)  is  calculated  as: 


using  von  Karman's  constant  (A:),  friction  velocity  (w*),  reference  height  (z),  roughness 
length  (zo),  and  the  stability  function  (<//)  (Brook  et  al.  1999).  Further  details  regarding 
calculation  of  these  parameters  can  be  found  in  Appendix  A. 

Sublayer  or  boundary  layer  resistance  {Rb)  is  calculated  as: 


using  friction  velocity  (w*),  thermal  diffusivity  (V,),  and  molecular  diffusivity  (D,)  (Brook 
etal.  1999). 

Total  surface  resistance  {Re)  of  the  ENVC  model  is  represented  by  stomatal  and  non- 
stomatal  resistances  (Zhang  et  al.  2003).  The  stomatal  resistance  components  include 
stomatal  {Rst)  and  mesophyll  {Rm)  resistance,  adjusted  for  stomatal  blocking  caused  by 
wet  conditions  {Wst)  (Zhang  et  al.  2003).  Non-stomatal  resistance  components  include  in- 
canopy  aerodynamic  resistance  {Rac),  ground  resistance  {Rg),  and  cuticle  resistance  {Rent) 
(Zhang  et  al.  2003).  The  relationships  between  these  resistances  are  illustrated  in  Figure 
4. 

Mathematically,  the  Re  portion  of  the  model  can  be  stated  as  follows  (adapted  from 
Zhang  et  al.,  2003): 


R^=±-[0nA\n{zlz,)-YH] 


(2) 


(3) 
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Pollutant  Conceuriaticn 


Soil 


Figure  4        Scheme  of  resistance  analogy  in  Environment  Canada  dry  deposition 
inference  model  (after  Zhang  et  al.  2003). 


Stomatal  resistance  (Rst)  includes  four  components  (Zhang  et  al.  2003):  unstressed  leaf 
stomatal  conductance  {Gs(PAR)),  conductance  reducing  effects  of  air  temperature  (f(T)), 
water-vapor  pressure  deficit  (f(D)X  and  leaf  water  potential  (f(ii/)).  Mathematically, 
stomatal  resistance  can  be  stated  as  (Zhang  et  al.  2003): 


"  G^{PAR)fiT)f{D)f{y/)DJD^, 

Mesophyll  resistance  (/?,„)  is  species  dependent  and  Environment  Canada  represents  this 
term  as  a  series  of  default  values  for  different  chemical  species  (refer  to  Table  A. 3  in 
Appendix  A)  (Zhang  et  al.  2003). 

Non  stomatal  resistances  include  (Appendix  A):  in-canopy  aerodynamic  resistance  {Rac), 
ground  resistance  {Rg),  and  canopy  cuticle  resistance  {Rent)  (Zhang  et  al.  2003).  In- 
canopy  aerodynamic  resistance  is  represented  as  a  series  of  default  values  adjusted  for 
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leaf  area  indices  (LAIs)  for  different  seasons.  Ground  resistance  is  represented  as  a  series 
of  default  values  adjusted  for  temperature,  rain  or  dew  conditions,  and  snow  cover 
fraction  (Zhang  et  al.  2003).  Cuticle  resistance  is  represented  as  a  series  of  default  values 
according  to  land  use  category  and  adjusted  for  temperature,  rain  or  dew  conditions,  leaf 
area  indices,  and  snow  cover  fraction  (Zhang  et  al.  2003). 


3,1,2  Particulates 

The  ENVC  model  uses  a  separate  module  for  modeling  particle  deposition  velocity  (Vj); 
which  includes  components  of  gravitational  settling  (V^),  aerodynamic  resistance  (/?«), 
and  surface  resistance  {Rs)  in  the  following  relationship  (Zhang  et  al.  2001): 

Gravitational  setthng  velocity  (V^)  includes  the  following  variables:  particle  density  (/?), 
particle  diameter  {dp),  gravitational  acceleration  (g),  viscosity  (41),  and  the  Cunningham 
slip  correction  factor  (Q  in  the  following  relationship  (Zhang  et  al.  2001): 

(7) 

^  18// 

Aerodynamic  resistance  (/?«)  is  calculated  the  same  way  as  for  gaseous  deposition  in  the 
ENVC  model  (Zhang  et  al.  2001).  Surface  resistance  {Rs)  includes  collection  efficiency 
of  Brownian  diffusion  {Eb),  impaction  {Em),  interception  {Ein),  a  stickiness  factor  {Rj), 
friction  velocity  (m*),  and  a  constant  {co)  as  illustrated  in  the  following  relationship 
(Zhang  et  al.  2001): 

R,  =  ^   (8) 

£qU,{Ei^  +Ej^  +E,^)R^ 

For  actual  apphcation  of  the  model,  size  distributions  used  for  a  particular  particle  group 
is  divided  into  40  incremental  sizes  using  the  following  formula,  where  I  is  from  1  to  40: 

size  _  increment  =  0.00 1  *  /  (9) 

The  resulting  particle  sizes  range  from  0.001  |nm  for  I  =  1 ;  to  10  |Lim  for  I  =  40.  Once  a 
size  increment  has  been  generated,  it  is  used  to  calculate  the  deposition  velocity 
associated  with  that  size.  The  resulting  deposition  for  each  size  increment  is  weighted 
according  to  the  area  under  the  curve  for  that  size  increment  and  then  summed  across  the 
forty  increment  sizes  for  the  total  deposition  velocity  (Zhang  et  al.  2001).  A  similar 
process  is  used  for  the  distribution  of  larger  particles. 
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3.1.3      Alberta  En  vironment  (AENV)  Model 


3.1.4  Gases 

A  description  of  the  AENV  model  is  provided  by  Cheng  et  al.  (2001 )  and  WBK  (2008) 
and  briefly  described  here.  The  AENV  model  uses  the  same  resistance  analogy  approach 
for  calculation  of  deposition  velocity  (V^)  as  the  ENVC  model,  but  calculates  each  of  the 
three  components  differently  (Appendix  B)  (WBK,  2008;  Cheng  et  al.  2001 ).  The  AENV 
calculation  for  Ra  involves  calculating  a  value  for  the  surface  roughness  length  (zo)  as  a 
monthly  average  using  hourly  data  where  the  wind  speed  is  >6  m/s: 


The  AENV  model  calculates     based  upon  von  Karman's  constant  (/:),  the  friction 
velocity  (w*),  the  Schmidt  number  iSc)  for  a  chemical,  and  the  Prandt  number  (Pr)  for  air 
in  the  following  relationship  (Cheng  et  al.,  2001): 


The  AENV  model  represents  Rc  as  default  values  that  are  seasonally  adjusted  for  the 
gases  SO2  and  NO2.  For  HNO2  and  HNO3,  Rc  is  always  10  s/m.  The  default  values  for 
SO2  and  NO2  are  drawn  from  published  literature  as  a  function  of  surface  type,  surface 
wetness,  and  incident  radiation;  and  as  such  address  some  of  the  variation  in  R,  caused  by 
stomatal  physiology.  The  AENV  model  incorporates  differences  for  wet  and  dry  surface 
conditions,  weighted  based  on  surface  wetness  (Appendix  B).  If  surface  wetness 
measurement  data  are  unavailable,  then  relative  humidity  is  used  a  surrogate  for  surface 
wetness.  The  model  allows  a  minimum  relative  humidity  (RH)  measurement  cutoff 
(87%)  to  be  used  as  an  approximation  to  represent  a  surface  as  being  wet. 

3.1.5  Particulates 

The  AENV  model  uses  the  same  deposition  velocity  relationship  for  particulates  as  used 
by  the  gaseous  model.  Initially,  use  of  this  approach  may  appear  to  be  quite  different 
from  that  of  the  ENVC  model.  However,  the  two  models  are  actually  quite  similar  as  the 
gravitational  settling  velocity  component  of  the  ENVC  model  is  considered  to  be  quite 
small. 


(10) 


(11) 


(12) 
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For  particulates,  Rc  is  considered  negligible  and  this  resistance  is  assigned  a  value  of  0. 
Ra  is  determined  independent  of  chemical  species.  Rt  is  represented  as  default  values  for 
different  species  and  are  drawn  from  scientific  literature  and  adjusted  for  dry/wet  and 
seasonal  conditions  as  a  function  of  surface  type  and  weighted  according  to  average  day 
length  for  each  month  at  a  mid- Alberta  latitude  location  (54°N  latitude). 
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4.0        DRY  DEPOSITION  MODELING  FOR  ENVC  AND  AENV 
MODELS 


4.1  Model  Methods  for  Coniferous  Forest  (Base  Case) 

The  two  dry  deposition  inference  models  in  their  original  form,  described  in  Appendices 
A  and  B,  were  run  to  establish  a  starting  point  from  which  to  evaluate  effects  of  changes 
made  to  the  AENV  model,  hiitially,  results  were  shown  for  monthly  deposition  flux  to 
permit  a  comparison  of  the  overall  performance  of  the  two  models  for  a  coniferous  forest 
land  use  category  (LUC)  and  considering  1 1  species  (SO2,  S04~ \  NO2,  HNO3,  HNO2, 
NOs',  NH/,  Ca""",  Na^  K^  and  Mg"^). 

Results  were  further  examined  using  deposition  velocities  (Vd)  for  each  species  involved 
in  the  calculation  of  monthly  deposition  flux.  Since  deposition  velocities  should  be 
relatively  constant  for  all  chemical  species  for  a  given  land  use  and  set  of  meteorological 
conditions,  comparisons  with  deposition  velocity  literature  values  were  made  for 
individual  species.  The  behavior  of  the  three  most  important  individual  species  making 
up  the  largest  difference  in  monthly  deposition  flux  between  the  two  models  was  further 
examined  to  better  understand  the  importance  of  these  species.  These  species  included 
SO2,  NO2,  and  HNO3.  Using  the  "contribution  to  total  deposition  flux"  rather  than 
"deposition  velocity"  of  an  individual  species  allowed  an  understanding  of  how  important 
each  species  was  in  contributing  to  change  in  total  deposition  flux  when  modifications  to 
the  AENV  model  were  made. 

The  original  model  form  for  the  AENV  model  was  also  used  to  generate  results  for  seven 
other  LUCs  for  which  the  AENV  model  was  designed.  These  results  were  compared 
with  original  results  for  a  coniferous  forest  LUC.  These  comparisons  were  done  for 
monthly  deposition  flux;  as  well  as  for  surface  resistance  (Rc)  values  for  gases  and 
deposition  velocity  (Vd)  values  for  particles.  It  is  important  to  point  out  that  inputs  for 
meteorological  parameters  used  in  comparing  the  two  dry  deposition  inference  models 
were  obtained  from  a  setting  consistent  with  a  coniferous  forest  LUC. 

These  parameter  inputs  were  used  to  make  comparisons  for  the  seven  LUCs  to  better 
understand  their  general  behavior  for  the  AENV  model.  These  parameter  inputs  are  not 
necessarily  representative  of  the  settings  for  the  seven  LUCs  because  some 
meteorological  parameters  -  such  as  roughness  length  and  relative  humidity  -  are 
affected  by  the  respective  surface  conditions.  As  such,  deposition  velocities  and  surface 
resistance  values  cannot  be  fairly  compared  with  literature  values  for  these  seven  LUCs. 

4,1.1        Input  Parameters  for  Modeling 

All  of  the  modeling  in  this  study  was  performed  using  measurements  of  hourly 
pollutant  concentrations  for  1 1  species  and  hourly  values  for  meteorological  parameters 
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made  at  the  Wood  Buffalo  Environmental  Association  Fort  Mac  Kay  air  monitoring 
station  for  the  year  2003.  Gaseous  species  measured  included  SO2,  NO2,  HNO3,  and 
HNO2.  Particulate  ions  included  SO4"",  NO3',  and  NH/.  Base  cations  included  K^,  Na^, 
CaT^,  and  Mg"^.  Measured  meteorological  values  included  wind  speed,  standard 
deviation  of  wind  direction,  air  temperature  at  2  m  and  10  m  heights,  and  relative 
humidity.  Inputs  for  solar  zenith  angle  and  solar  radiation  values  used  in  the  ENVC 
model  were  computed  using  a  Visual  Basic  program  developed  during  a  previous  study 
(WBK,  2008). 

Quality  assurance  and  quality  control  was  performed  on  the  input  data.  Input  values 
found  to  be  very  uncommon  were  left  in  the  input  data  set  to  allow  performance  of  the 
models  to  be  evaluated  using  actual  data.  Missing  hourly  data  were  replaced  with  the 
average  of  the  value  for  the  hours  before  and  after  the  hour  of  the  missing  value  if  only 
one  hour's  value  was  missing.  If  two  or  more  consecutive  hours  were  missing,  then  the 
month's  median  value  for  each  of  those  hours  was  used  to  represent  the  missing  data  for 
those  hours.  For  species  measured  intermittently  (24  hours  every  6  days),  the  hourly 
average  was  used  for  all  of  the  144  hours  (6  days)  in  that  time  period.  Complete  details 
of  data  collection,  processing,  and  assumptions  are  described  in  WBK  (2008). 

4,1,2        A  tmospheric  Concentrations  of  Chemical  Species 

Deposition  of  a  particular  species  is  the  product  of  species  concentration  and  deposition 
velocity.  From  a  practical  point  of  view,  it  is  important  to  understand  those  species 
having  the  greatest  influence  on  acid  deposition  in  terms  of  hydrogen  ion  equivalents. 
Monthly  averages  species  concentrations  are  provided  in  Appendix  C. 

Concentrations  of  nitrogen  compounds  (NO2,  HNO3,  HNO2,  NO3',  and  NH4^)  tend  to  be 
higher  in  winter  than  summer  months.  Deposition  velocities  for  NO2,  HNO3,  and  HNO2 
tend  to  be  lower  in  winter  months.  Species  and  respective  concentration  ranges  were  as 
follows  for  the  2003  data  set: 

•  NO2  was  the  species  with  the  highest  recorded  hourly  concentration  (range  3.7  to 
21.7  ^gW). 

•  SO2  was  the  second  most  important  species  although  hourly  levels  observed  were 
much  lower  (range  1.0  to  4.9  pg/m^). 

•  Calcium  ion  (Ca^^)  levels  ranged  from  0.2  to  1.9  lugW.  Calcium  concentration 
maxima  occurred  in  April,  July,  and  November. 

•  Sulphate  (S04^")  ions  ranged  from  0.5  to  2.4  |ug/ml 

•  HNO3  concentration  ranged  from  0.2  to  1.5  |igW. 

•  Levels  of  all  other  species  were  less  than  1.0  |ug/m^  for  all  months  of  the  year. 
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4.2  Model  Results  -  Coniferous  Forest 


4,2 A        Deposition  Flux 

The  first  modeling  was  done  using  AENV  and  ENVC  models  in  their  original  forms. 
Total  monthly  deposition  flux  values  ranged  from  near  0  to  0.016  kg  hydrogen  ion  per 
hectare  (kg  H'^/ha)  for  the  AENV  model  and  from  0.004  to  0.018  kg  H^ha  for  the  ENVC 
model  (Table  2).  Total  deposition  flux  for  the  year  was  0.141  kg  HVha  for  the  AENV 
model  and  0.153  kg  H^/ha  for  the  ENVC  model. 

For  total  annual  deposition,  the  difference  between  the  two  models  was  small  at  0.012  kg 
H^/ha,  representing  a  difference  of  about  8%.  As  shown  in  Figure  5,  months  having  the 
largest  differences  were  November,  July,  and  September  -  all  of  which  the  AENV  model 
underestimated  the  ENVC  model.  May,  October,  April,  and  December  showed  moderate 
differences.  Total  deposition  for  the  months  of  January,  February,  March,  June,  and 
August  were  similar  for  both  models. 


Table  2         Comparison  of  monthly  summed  deposition  flux  (kg  H^/ha)  for  1 1 
species  to  a  coniferous  forest  cover  using  the  AENV  and  ENVC 
models. 


Month 

AENV 

ENVC 

Difference 

Month  1 

0.016 

0.016 

0.000 

Month  2 

0.011 

0.011 

-0.001 

Month  3 

0.016 

0.016 

0.000 

Month  4 

0.016 

0.017 

-0.001 

Month  5 

0.014 

0.012 

0.002 

Month  6 

0.007 

0.006 

0.001 

Month  7 

0.007 

0.011 

-0.003 

Month  8 

0.013 

0.013 

0.000 

Month  9 

0.010 

0.012 

-0.003 

Month  10 

0.016 

0.018 

-0.002 

Month  1 1 

0.000 

0.004 

-0.004 

Month  12 

0.016 

0.017 

-0.001 

Total 

0.141 

0.153 

-0.012 

Total  %  difference 

8 
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Figure  5        Comparison  of  monthly  summed  deposition  flux  for  11  species  to  a 
coniferous  forest  cover  using  the  AENV  and  ENVC  models. 


4,2,2        Deposition  Ve/odties  -  Coniferous  Forest 

The  initial  model  results  were  broken  out  to  show  monthly  average  deposition  velocities 
for  each  species  in  Figures  6  to  11.  Tabular  values  associated  with  Figures  6  to  1 1  can  be 
found  in  Appendix  D. 

4.2.2.1       Deposition  Beliavior  of  Gases  (SO2,  NO2,  HNO3,  HNO2) 
SO2 

Resulting  average  deposition  velocities  (Vd)  for  SO2  for  the  two  models  indicate  that  SO2 
deposition  ranges  from  0.09  to  0.49  cm/s  across  the  year.  The  AENV  model  under 
represents  the  ENVC  model  in  winter  and  spring  months,  but  over  represents  it  in 
summer  months  (Figure  6). 
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Figure  6        Comparison  of  monthly  average  SO2  deposition  velocity  for  a 
coniferous  forest  cover  using  the  AENV  and  ENVC  model. 


In  a  study  from  literature,  a  U.S.  National  Dry  Deposition  Network  (NDDN)  model  - 
which  is  reported  to  be  similar  to  the  ENVC  model  -  was  used  to  simulate  acidic 
deposition  at  Ann  Arbor,  Michigan,  for  a  land  use  combination  of  grass,  water,  and  forest 
(Brook,  1996).  Their  model  results  indicated  SO2  monthly  average  deposition  velocities 
ranged  from  0.20  to  0.55  cm/s,  with  the  peak  in  June  (Brook,  1996). 

NO2 

NO2  monthly  average  deposition  velocities  are  shown  in  Figure  7  and  ranged  from  0.04 
to  0.25  cm/s  with  the  AENV  model  under  representing  the  ENVC  model  in  winter  and 
spring  months  while  over  representing  it  in  summer.  Generally,  ENVC  results  were  quite 
similar  to  AENV  results. 
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Figure  7 


Comparison  of  monthly  average  NO2  deposition  velocity  for  a 
coniferous  forest  cover  using  the  AENV  and  ENVC  model. 


HNO3 

Monthly  average  HNO3  deposition  velocities  are  shown  in  Figure  8  and  they  range  from 
0.8  to  2.1  cm/s  with  the  AENV  model  over  representing  the  ENVC  model  in  all  months 
except  December.  The  two  models  produced  similar  trends  with  monthly  average  AENV 
deposition  velocities  being  higher  than  those  of  the  ENVC  model.  In  the  previously 
mentioned  Ann  Arbor  study  (Brook,  1996),  the  NDDN  model  produced  HNO3  deposition 
velocities  ranging  from  1.4  to  2.2  cm/s  with  peaks  in  April  and  December  and  troughs  in 
August. 
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Figure  8        Comparison  of  monthly  average  HNO3  deposition  velocity  for  a 
coniferous  forest  cover  using  the  AENV  and  ENVC  model. 


HNO2 

Monthly  average  HNO2  deposition  velocities  are  shown  in  Figure  9  and  they  ranged  from 
0.3  to  2.2  cm/s.  The  AENV  model  consistently  over  represented  the  ENVC  model  by 
two  to  three  times.  While  this  difference  may  seem  substantial,  the  low  concentration  of 
HNO2  in  the  atmosphere  results  in  a  relatively  small  contribution  to  total  monthly 
deposition  flux  for  both  of  the  models. 
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Figure  9        Comparison  of  monthly  average  HNO2  deposition  velocity  for  a 
coniferous  forest  cover  using  the  AENV  and  ENVC  model. 


4.2.2.2       Fine  Particle  (S04^",  NH/)  Deposition  Behavior 

Deposition  velocities  for  smaller  particles  (S04^'  and  NH/)  are  shown  in  Figure  10  and 
they  ranged  from  0.09  to  0.35  cm/s  with  the  AENV  model  over  representing  the  ENVC 
model  in  every  month.  Particle  deposition  behavior  in  both  models  is  mostly  dominated 
by  aerodynamic  resistance  (Ra)  based  on  aerodynamic  properties  (diameter)  of  the 
particle  and  independent  of  individual  species. 

The  AENV  results  showed  a  peak  in  summer  months  that  was  absent  in  the  ENVC 
model.  In  the  literature  study  in  Ann  Arbor  (Brook,  1996),  the  NDDN  model  reported 
small-particle  monthly  average  deposition  velocities  ranging  from  0.07  to  0.15  cm/s  with 
slightly  higher  values  in  May  and  a  minimum  in  winter  months. 
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Figure  10      Comparison  of  monthly  average  small  particle  deposition  velocity  for 
a  coniferous  forest  cover  using  the  AENV  and  ENVC  model. 


4.2.2.3       Coarse  Particle  (NO3",  Na*,  K*,  Ca  ,  Mg  )  Deposition  Behavior 


Again,  particle  deposition  behavior  in  both  models  is  mostly  dominated  by  aerodynamic 
resistance  (Ra)  based  on  aerodynamic  properties  (diameter)  and  independent  of  individual 
species.  Monthly  average  deposition  velocities  for  larger  (coarse)  particles  shown  in 
Figure  1 1  and  they  ranged  from  0.29  to  0.64  cm/s.  Similar  to  Figure  10  (small  particle 
deposition),  monthly  average  deposition  velocities  for  large  particles  showed  a  peak  in 
sunamer  months  for  the  AENV  model  that  did  not  occur  in  the  ENVC  model. 


Refinement  Study  of  Dry  Deposition  Inference 
Methods  Used  in  Alberta 


25 


0.70 


Figure  11      Comparison  of  monthly  average  large  particle  deposition  velocity  for 
a  coniferous  forest  cover  using  the  AENV  and  ENVC  model. 


Monthly  averaged  diurnal  deposition  velocities  for  the  three  most  dominant  species  - 
SO2,  NO2,  and  HNO3  -  were  calculated  and  compared  for  each  model.  These 
comparisons  are  presented  in  Appendix  E.  hi  all  cases,  both  models  were  similar  and 
behaved  as  expected  in  that  higher  monthly  averaged  diurnal  deposition  velocities 
occurred  during  daytime  heating  hours  and  lower  monthly  averaged  diurnal  deposition 
velocities  occurred  during  nighttime. 

4.3  Model  Methods  for  Other  Land  Use  Categories 

The  AENV  model  was  rerun  for  seven  land  uses  that  represent  other  land  use  categories 
(LUCs)  for  which  the  model  was  designed,  hiput  data  were  the  same  as  that  used 
previously,  the  only  change  was  the  selection  of  LUC.  Because  calculation  of  roughness 
length  (Zo)  is  based  on  the  meteorological  input  data,  roughness  lengths  used  for  these 
seven  LUCs  were  the  same  as  that  used  for  a  coniferous  forest  LUC. 


Refinement  Study  of  Dry  Deposition  Inference 
Methods  Used  in  Alberta 


26 


4.4  Model  Results  -  Coniferous  Forest 


4.4.1        Deposition  Flux 

Monthly  summed  deposition  flux  results  for  all  of  LUCs  are  compared  with  a  coniferous 
forest  surface  cover  in  Figures  12  through  18  and  detailed  for  individual  species  in  Tables 
3  through  9.  A  negative  value  in  the  tables  indicates  that  the  monthly  average  result  for 
an  individual  species  for  the  alternative  LUC  is  less  than  for  a  coniferous  forest  surface 
(the  value  indicating  the  magnitude  of  the  difference)  and  vice  versa. 

LUC  1  -  Deciduous  Forest  Compared  with  Coniferous  Forest 

The  differences  in  monthly  summed  deposition  flux  for  a  deciduous  forest  cover  are 
shown  in  Figure  12.  From  Table  3,  it  is  observed  the  largest  differences  are  for  NO2  and 
Ca^^  deposition  in  winter  and  Ca'^^  deposition  in  summer. 
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Figure  12      Comparison  of  monthly  summed  deposition  flux  for  11  species  to  a 
deciduous  forest  and  coniferous  forest  cover  using  the  AENV  model. 
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Table  3         Breakdown  of  magnitude  of  differences  in  monthly  summed 

deposition  flux  in  kg  H"^/ha  for  individual  species  between  deciduous 
forest  and  coniferous  forest  cover  using  the  AENV  model. 


Month 

SO2 

NO2 

HNO2 

HN03 

S04^" 

NO3 

Mg2^ 

Total 

Month  1 

-0.5 

-2.5 

0.0 

0.0 

-0.7 

-1.2 

-0.3 

0.4 

0.1 

0.6 

0.2 

-3.9 

Month  2 

-0.7 

-2.5 

0.0 

0.0 

-0.8 

-1.0 

-0.2 

0.3 

0.0 

0.9 

0.2 

-3.7 

Month  3 

-0.7 

-1.8 

0.0 

0.0 

-1.3 

-1.0 

-0.1 

0.2 

0.0 

0.8 

0.1 

-3.6 

Month  4 

-0.6 

-1.8 

0.0 

0.0 

-0.6 

-0.8 

-0.1 

0.2 

0.0 

2.2 

0.2 

-1.3 

Month  5 

-0.3 

-1.1 

0.0 

0.0 

-0.4 

-0.5 

0.0 

0.1 

0.0 

1.0 

0.2 

-1.0 

Month  6 

0.0 

0.0 

0.0 

0.0 

0.4 

0.5 

0.0 

-0.1 

-0.1 

-1.4 

-0.2 

-0.8 

Month  7 

0.0 

0.0 

0.0 

0.0 

0.4 

0.4 

0.1 

-0.1 

-0.1 

-2.5 

-0.3 

-2.2 

Month  8 

0.0 

0.0 

0.0 

0.0 

0.5 

0.6 

0.1 

-0.1 

-0.1 

-1.8 

-0.4 

-1.1 

Month  9 

-0.7 

-1.5 

0.0 

0.0 

-0.2 

-0.3 

0.0 

0.1 

0.0 

0.4 

0.1 

-2.1 

Month  10 

-0.7 

-1.7 

0.0 

0.0 

-0.3 

-0.4 

0.0 

0.1 

0.0 

1.2 

0.1 

-1.6 

Month  1 1 

-0.3 

-2.5 

0.0 

0.0 

-0.3 

-0.6 

-0.2 

0.1 

0.1 

3.1 

0.1 

-0.6 

Month  12 

-1.2 

-3.9 

0.0 

0.0 

-0.8 

-1.0 

-0.6 

0.4 

0.2 

2.6 

0.4 

-3.8 

LUC  3  -  Wetland/Swamp  Compared  with  Coniferous  Forest 

The  differences  in  monthly  summed  deposition  flux  for  a  wetland/swamp  are  shown  in 
Figure  13.  From  Table  4,  it  is  observed  that  the  largest  differences  are  for  NO2 
deposition  being  low  across  the  year  because  of  wet  surface  conditions  representative  of 
wetland/swamp  and  SO2  deposition  being  higher  in  spring  and  summer  months  because 
of  wet  surface  conditions  representative  of  wetland/swamp  compared  to  winter  months. 

LUC  4  -  Grassland  as  Compared  with  Coniferous  Forest 

The  differences  in  monthly  summed  deposition  flux  for  grassland  are  shown  in  Figure  14. 
The  trend  in  monthly  summed  deposition  flux  for  this  surface  type  was  similar  to  that  for 
a  wetland/swamp  cover  (Figure  13).  Importance  of  individual  species  and  underlying 
conditions  contributing  to  these  differences  can  be  understood  from  Table  5.  Again,  the 
largest  differences  are  for  SO2  and  NO2  deposition. 
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Figure  13      Comparison  of  monthly  summed  deposition  flux  for  11  species  to 

wetland/swamp  and  coniferous  forest  cover  using  the  AENV  model. 


Table  4         Breakdown  of  magnitude  of  differences  in  monthly  summed 
deposition  flux  in  kg  H^/ha  for  individual  species  between 
wetland/swamp  and  coniferous  forest  cover  using  the  AENV  model. 


Month 

SO, 

NO. 

HNO2 

HN03 

S04- 

NH4^ 

NO3 

Mg-^ 

Total 

Month  1 

-0.3 

-4.2 

0.0 

0.0 

-0.7 

-1.2 

-0.3 

0.4 

0.1 

0.6 

0.2 

-5.4 

Month  2 

-0.4 

-4.2 

0.0 

0.0 

-0.8 

-1.0 

-0.2 

0.3 

0.0 

1.0 

0.2 

-5.1 

Month  3 

3.8 

-5.2 

0.0 

0.0 

-0.7 

-0.6 

-0.1 

0.1 

0.0 

0.5 

0.1 

-2.1 

Month  4 

11.2 

-6.8 

0.0 

0.0 

-0.3 

-0.5 

-0.1 

0.1 

0.0 

1.2 

0.1 

5.0 

Month  5 

12.8 

-4.7 

0.0 

0.0 

-0.2 

-0.3 

0.0 

0.1 

0.0 

0.6 

0.1 

8.3 

Month  6 

6.1 

-3.8 

0.0 

0.0 

-0.2 

-0.3 

0.0 

0.0 

0.0 

0.9 

0.1 

2.9 

Month  7 

12.4 

-5.3 

0.0 

0.0 

-0.2 

-0.3 

-0.1 

0.1 

0.1 

1.6 

0.2 

8.6 

Month  8 

10.5 

-6.9 

0.0 

0.0 

-0.3 

-0.4 

0.0 

0.1 

0.0 

1.1 

0.2 

4.3 

Month  9 

4.1 

-3.9 

0.0 

0.0 

-0.2 

-0.3 

0.0 

0.1 

0.0 

0.4 

0.1 

0.3 

Month  10 

3.7 

-4.5 

0.0 

0.0 

-0.3 

-0.4 

0.0 

0.1 

0.0 

1.2 

0.1 

-0.1 

Month  1 1 

0.2 

-4.8 

0.0 

0.0 

-0.3 

-0.6 

-0.2 

0.1 

0.1 

3.1 

0.1 

-2.4 

Month  12 

-0.7 

-6.4 

0.0 

0.0 

-0.8 

-1.0 

-0.6 

0.4 

0.2 

2.7 

0.4 

-5.7 

Refinement  Study  of  Dry  Deposition  Inference 
IVIethods  Used  in  Alberta 


29 


5  0.025 


Figure  14 


6  7 
Month 


Grassland  "^Coniferous  Forest 


Comparison  of  monthly  summed  deposition  flux  for  11  species  to 
grassland  and  coniferous  forest  cover  using  the  AENV  model. 


Table  5         Breakdown  of  magnitude  of  differences  in  monthly  summed 

deposition  flux  in  kg  H'^/ha  for  individual  species  between  grassland 
and  coniferous  forest  cover  using  the  AENV  model. 


Month 

SO2 

NO2 

HNO2 

HN03 

S04' 

NH4^ 

NO3 

Na^ 

Ca^^ 

Total 

Month  1 

-0.3 

-4.2 

0.0 

0.0 

-0.7 

-1.2 

-0.3 

0.4 

0.1 

0.6 

0.2 

-5.4 

Month  2 

-0.4 

-4.2 

0.0 

0.0 

-0.8 

-1.0 

-0.2 

0.3 

0.0 

1.0 

0.2 

-5.1 

Month  3 

1.8 

-2.5 

0.0 

0.0 

-1.3 

-1.0 

-0.1 

0.3 

0.0 

0.8 

0.1 

-1.9 

Month  4 

5.8 

-2.1 

0.0 

0.0 

-0.6 

-0.8 

-0.1 

0.2 

0.0 

2.3 

0.2 

4.8 

Month  5 

6.4 

-1.1 

0.0 

0.0 

-0.4 

-0.6 

0.0 

0.1 

0.0 

1.1 

0.2 

5.7 

Month  6 

2.9 

-1.4 

0.0 

0.0 

-0.5 

-0.6 

-0.1 

0.1 

0.1 

1.7 

0.3 

2.6 

Month  7 

5.9 

-2.0 

0.0 

0.0 

-0.5 

-0.5 

-0.1 

0.2 

0.1 

3.2 

0.4 

6.6 

Month  8 

5.1 

-2.6 

0.0 

0.0 

-0.6 

-0.8 

-0.1 

0.2 

0.1 

2.2 

0.4 

4.0 

Month  9 

1.8 

-2.5 

0.0 

0.0 

-0.4 

-0.6 

-0.1 

0.2 

0.1 

0.8 

0.2 

-0.6 

Month  10 

1.7 

-3.0 

0.0 

0.0 

-0.6 

-0.7 

-0.1 

0.1 

0.1 

2.3 

0.3 

0.2 

Month  1 1 

0.0 

-4.2 

0.0 

0.0 

-0.5 

-1.2 

-0.4 

0.3 

0.1 

5.8 

0.3 

0.0 

Month  12 
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LUC  5  -  Cropland  as  Compared  with  Coniferous  Forest 

Monthly  summed  deposition  flux  for  cropland  is  shown  in  Figure  15  and,  as  expected,  it 
shows  a  large  increase  in  the  spring  (growing)  season.  It  should  be  noted  that  the  scale  in 
this  figure  has  been  changed  to  fit  the  large  increase  associated  with  growing  season 
conditions.  From  Table  6,  it  can  be  seen  that  this  increase  is  solely  a  result  of  changes 
(increases)  in  SO2  deposition.  Monthly  summed  deposition  fluxes  are  comparable  to  that 
of  a  coniferous  forest  surface  for  other  months  of  the  year. 


=5 


too 


o 

Oh 
0) 
-O 

-a 

B 
B 


0.090 


ffi  0.080 


^  0.070 


0.060 


2  0.050 


0.040 
0.030 
0.020 
0.010 


5  0.000 


§  -0.010 


Figure  15 


6  7 
Month 


10  11 


12 


Cropland        "*~  Coniferous  Forest 


Comparison  of  monthly  summed  deposition  flux  for  11  species  to 
cropland  and  coniferous  forest  cover  using  the  AENV  model. 


LUC  6  -  Urban  as  Compared  with  Coniferous  Forest 

Conditions  representing  an  urban  surface  resulted  in  similar  monthly  summed  deposition 
fluxes  during  winter  and  lower  values  during  other  periods  compared  to  a  coniferous 
forest  cover  (Figure  16).  From  Table  7,  differences  in  monthly  summed  deposition 
fluxes  from  March  through  October  can  be  seen  to  be  largely  a  result  of  differences  in 
NO2  deposition. 
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Table  6         Breakdown  of  magnitude  of  differences  in  monthly  summed 

deposition  flux  in  kg  H^/ha  for  individual  species  between  cropland 
and  coniferous  forest  cover  using  the  AENV  model. 


Month 
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LUC  7  -  Water  as  Compared  with  Coniferous  Forest 

As  shown  in  Figure  17,  water  has  the  abiHty  to  receive  much  greater  acid  input  than  a 
coniferous  forest  cover.  From  Table  8,  it  is  observed  that  this  is  largely  a  result  of 
differences  in  deposition  of  SO2,  NO2,  SO4'",  NH4^  NO3",  Na^,  Ca-^  and  Mg'^.  In 
general,  water  is  far  more  susceptible  to  SO2  deposition  and  less  susceptible  to  NO2 
deposition  than  a  coniferous  forest  surface.  Small  particle  deposition  (S04^',  NH4^)  is 
characteristically  higher  over  water  compared  to  a  coniferous  forest  cover.  Because 
water  has  Httle  or  no  resistance  to  SO2  or  particle  deposition,  total  deposition  is  largely 
dependent  on  atmospheric  concentration  of  these  species. 
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Figure  16      Comparison  of  monthly  summed  deposition  flux  for  11  species  to 
urban  and  coniferous  forest  cover  using  the  AENV  model. 


Table  7         Breakdown  of  magnitude  of  differences  in  monthly  summed 

deposition  flux  in  kg  H"^/ha  for  individual  species  between  urban  and 
coniferous  forest  cover  using  the  AENV  model. 


Month 

SO2 

NO2 

HNO2 

HN03 

S04^ 

NH4^ 

NO3- 

Ca-^ 

Mg-" 

Total 

Month  1 

-0.5 

0.0 

0.0 

0.0 

-0.8 

-1.3 

-0.3 

0.5 

0.1 

0.6 

0.2 

-1.4 

Month  2 

-0.7 

0.0 

0.0 

0.0 

-0.9 

-1.1 

-0.2 

0.3 

0.0 

1.0 

0.2 

-1.2 

Month  3 

-1.0 

-3.3 

0.0 

0.0 

-1.9 

-1.5 

-0.2 

0.4 

0.0 

1.3 

0.2 

-6.1 

Month  4 

-1.7 

-5.5 

0.0 

0.0 

-0.9 

-1.3 

-0.2 

0.3 

0.1 

3.4 

0.2 

-5.6 

Month  5 

-1.5 

-4.3 

0.0 

0.0 

-0.7 

-0.8 

-0.1 

0.2 

0.1 

1.6 

0.3 

-5.1 

Month  6 

-1.3 

-3.6 

0.0 

0.0 

-0.2 

-0.3 

0.0 

0.0 

0.0 

0.9 

0.1 

-4.4 

Month  7 

-2.6 

-5.1 

0.0 

0.0 

-0.2 

-0.3 

-0.1 

0.1 

0.1 

1.6 

0.2 

-6.2 

Month  8 

-2.4 

-6.6 

0.0 

0.0 

-0.3 

-0.4 

0.0 

0.1 

0.0 

1.1 

0.2 

-8.3 

Month  9 

-0.3 

-3.5 

0.0 

0.0 

-0.6 

-0.7 

-0.1 

0.2 

0.1 

1.0 

0.3 

-3.5 

Month  10 

-0.5 

-4.0 

0.0 

0.0 

-0.7 

-0.8 

-0.1 

0.2 

0.1 

2.9 

0.3 

-2.6 

Month  1 1 

-0.3 

-1.4 

0.0 

0.0 

-0.7 

-1.5 

-0.6 

0.3 

0.1 

7.4 

0.4 

3.8 

Month  12 

-1.2 

0.0 

0.0 

0.0 

-0.8 

-1.1 

-0.6 

0.4 

0.2 

2.9 

0.4 

0.2 

Refinement  Study  of  Dry  Deposition  Inference 
Methods  Used  in  Alberta 


33 


^  0.20 

+ 


S  -0.05  ^  '  ^  ^  ^  ^  ^  ^  ^  '  ^  1 

1        2       3       4       5       6       7       8       9       10      11  12 
 Month  


^  Water  "^Coniferous  Forest 


Figure  17      Comparison  of  monthly  summed  deposition  flux  for  11  species  to 
water  and  coniferous  forest  cover  using  the  AENV  model. 


Table  8         Breakdown  of  magnitude  of  differences  in  monthly  summed 

deposition  flux  in  kg  H"^/ha  for  individual  species  between  water  and 
coniferous  forest  cover  using  the  AENV  model. 


Month 
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LUC  8  -  Snow/Ice  as  Compared  with  Coniferous  Forest 

For  a  snow/ice-covered  surface,  the  months  of  January  and  February  are  repeated  at  the 
end  of  the  year  in  Figure  18  to  allow  for  continuity  of  the  data.  Snow/ice  cover  resuhs  in 
lower  monthly  summed  deposition  fluxes  compared  to  a  coniferous  forest  cover.  From 
Table  9,  it  can  be  seen  that  difference  are  largely  driven  by  SO2,  NO2,  and  small  particle 
deposition. 
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Comparison  of  monthly  summed  deposition  flux  for  11  species  to 
snow/ice  and  coniferous  forest  cover  using  the  AENV  model. 


4.4,2        Resistance  Factors  and  Deposition  Velocities 


Aerodynamic  resistance  (Ra)  is  dependent  on  meteorological  conditions  and  independent 
of  deposition  species  for  all  LUCs  in  the  AENV  model.  Calculation  of  boundary  layer 
resistance  (Rb)  for  gases  is  dependent  on  meteorological  conditions  and  individual 
species  and  independent  of  the  surface  (i.e.  LUC);  consequently  there  are  no  differences 
in  Rb  values  from  one  LUC  to  another  for  individual  species.  Finally,  surface  (canopy) 
resistance  (Rc)  is  dependent  upon  individual  species  and  surface  type.  Model  results  for 
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Rc  for  major  depositing  gases  (SO2,  NO2,  and  HNO3)  are  examined  more  closely  in  the 
following  section. 


Table  9         Breakdown  of  magnitude  of  differences  in  monthly  summed 

deposition  flux  in  kg  H^/ha  for  individual  species  between  snow/ice 
and  coniferous  forest  cover  using  the  AENV  model. 
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4.4.2.1       Gaseous  Deposition  (SO2,  NO2,  HNO3) 
Surface  Layer  Resistance  (Rc)  for  SO2 

The  calculation  of  Rc  for  SO2  in  the  AENV  method  is  done  using  default  values  such  that 
it  represents  a  total  "dry"  condition,  total  "wet"  condition,  or  "weighted  wet"  condition 
using  relative  humidity  (RH)  and/or  surface  wetness  (SW)  measurement  data  for  a 
specific  LUC  and  season  (Appendix  A).  For  the  input  meteorological  data  set  used  in 
this  study,  the  absence  of  surface  wetness  data  resulted  in  the  wetness  condition  being 
based  upon  whether  the  relative  humidity  was  above  or  below  87%.  Seasons  are 
classified  as  winter  (December  to  February),  spring  (March  to  May),  summer  (June  to 
August),  and  autumn  (September  to  November).  The  exception  to  this  seasonal 
designation  is  that  if  the  temperature  for  any  hour  is  below  0°C  during  any  season,  then 
Rc  for  that  hour  is  represented  by  the  corresponding  wet/dry  value  for  winter  for  that 
species. 
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LUC  1  -  Deciduous  Forest 

Monthly  average  Rc  values  for  a  deciduous  forest  surface  cover  (shown  in  Figure  19)  are 
higher  in  winter  and  early  spring/late  autumn  than  for  a  coniferous  forest  surface  cover. 
This  would  be  expected  due  to  an  absence  of  leaves  during  these  periods.  Monthly 
average  Rc  values  for  both  covers  are  similar  during  the  leafy  season  (June,  July,  August). 


Figure  19 
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Comparison  of  SO2  monthly  average  surface  resistance  (Rc)  for  a 
deciduous  forest  and  coniferous  forest  (base  case)  cover  using  the 
AENV  model. 


LUC  3,  4,  and  5  -  Wetland/Swamp,  Grassland,  and  Cropland 

Monthly  average  Rc  values  for  wetlands/swamp  (Figure  20),  grassland  (Figure  21),  and 
cropland  (Figure  22)  are  all  similar  because  of  the  nature  of  the  growing  season  of 
vegetation  representing  these  cover  types  in  Alberta.  Monthly  average  Rc  values  are  the 
same  for  all  three  LUCs  for  the  winter  season;  and  they  are  all  greater  in  winter  compared 
to  that  for  a  coniferous  forest  cover  which  maintains  more  of  its  leaves  (needles)  during 
winter.  During  spring  and  sunmier  months,  monthly  average  Rc  values  for  a 
wetland/swamp  become  much  lower  than  those  of  the  coniferous  forest  cover  because  of 
the  growing  season  and  lushness  of  vegetation. 
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Figure  20      Comparison  of  SO2  monthly  average  surface  resistance  (Rc)  for 
wetland/swamp  and  coniferous  forest  (base  case)  cover  using  the 
AENV  model. 
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Figure  21  Comparison  of  SO2  monthly  average  surface  resistance  (Rc)  for  a 
grassland  and  coniferous  forest  (base  case)  cover  using  the  AENV 
model. 
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Figure  22  Comparison  of  SO2  monthly  average  surface  resistance  (Rc)  for  a 
cropland  and  coniferous  forest  (base  case)  cover  using  the  AENV 
model. 


LUC  6  -  Urban 

Monthly  average  Rc  values  for  urban  land  (shown  in  Figure  23)  are  higher  than  those  for 
a  coniferous  forest  cover;  which  may  be  a  result  of  building  and  road  surfaces  which  are 
not  as  interactive  to  atmospheric  SO2  as  is  vegetation.  The  dechning  resistance  as 
summer  progresses  to  autumn  is  mostly  associated  with  increasing  relative  humidity; 
which  is  illustrated  in  Table  10. 

LUC  7  -  Water 

Monthly  average  Rc  values  for  water  are  constant  and  set  at  zero  throughout  a  year 
(Figure  24). 

LUC  8  -  Snow/Ice 

Monthly  average  Rc  values  for  a  snow/ice  cover  are  shown  in  Figure  25.  These  values 
are  similar  as  those  for  a  snow-covered  grassland,  cropland,  or  wetland/swamp  surface. 
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Table  10       Monthly  average  relative  humidity  level  (%)  for  data  used  to  evaluate 
AENV  model. 
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Figure  23      Comparison  of  SO2  monthly  average  surface  resistance  (Rc)  for  urban 
and  coniferous  forest  (base  case)  cover  using  the  AENV  model. 
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Figure  24      Comparison  of  SO2  monthly  average  surface  resistance  (Rc)  for  a 

water  and  coniferous  forest  (base  case)  cover  using  the  AENV  model. 
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Figure  25      Comparison  of  SO2  monthly  average  surface  resistance  (Rc)  for 
snow/ice  and  coniferous  forest  (base  case)  cover  using  the  AENV 
model. 
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Surface  Layer  Resistance  (Rc)  for  NO2 

Similar  to  SO2,  Rc  for  NO2  is  dependent  upon  LUC,  season,  and  wetness  condition. 
Again,  the  wetness  condition  used  in  this  study  was  based  upon  whether  relative  humidity 
was  above  87%.  The  relationships  between  monthly  average  Rc  values  for  NO2  for 
different  land  uses  are  shown  in  Figures  26  to  32.  The  relationships  for  monthly  average 
Rc  values  for  NO2  were  similar  to  those  for  SO2  except  that  values  for  NO2  were 
generally  much  higher.  Monthly  average  Rc  values  for  NO2  are  more  sensitive  to 
differences  between  wet  and  dry  conditions  as  compared  to  SO2.  Increases  in  monthly 
average  Rc  values  for  NO2  across  spring  and  summer  months  generally  observed  for  all 
surface  types  reflects  increasing  relative  humidity  (as  shown  in  Table  10).  Even  though 
November  is  classified  as  an  autumn  month,  the  presence  of  numerous  hourly 
temperatures  less  than  0°C  results  in  November  behaving  as  a  winter  month. 
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Figure  26      Comparison  of  NO2  monthly  average  surface  resistance  (Rc)  for 
deciduous  forest  and  coniferous  forest  (base  case)  cover  using  the 
AENV  model. 
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Figure  27      Comparison  of  NO2  monthly  average  surface  resistance  (Rc)  for 
wetland/swamp  and  coniferous  forest  (base  case)  cover  using  the 
AENV  model. 
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Comparison  of  NO2  monthly  average  surface  resistance  (Rc)  for 
grassland  and  coniferous  forest  (base  case)  cover  using  the  AENV 
model. 
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Rc  cropland      ~*~  Rg  coniferous  forest 


Figure  29      Comparison  of  NO2  monthly  average  surface  resistance  (Rc)  for 
cropland  and  coniferous  forest  (base  case)  cover  using  the  AENV 
model. 
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Comparison  of  NO2  monthly  average  surface  resistance  (Rc)  for  urban 
and  coniferous  forest  (base  case)  cover  using  the  AENV  model. 
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Figure  31 
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Comparison  of  NO2  monthly  average  surface  resistance  (Rc)  for  water 
and  coniferous  forest  (base  case)  cover  using  the  AENV  model. 
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Comparison  of  NO2  monthly  average  surface  resistance  (Rc)  for 
snow/ice  and  coniferous  forest  (base  case)  cover  using  the  AENV 
model. 
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Surface  Layer  Resistance  (Rc)  for  HNO3 

Values  for  HNO3  surface  resistance  (Rc)  are  treated  as  a  constant  (10  s/m)  for  all  seasons 
and  all  surface  types. 

4.4.2.2       Fine  and  Coarse  Particle  Deposition  (S04^",  NH4'^,  NO3",  Na"^,  K\  Ca^*,  Mg^*) 

The  AENV  model  classifies  S04^'  and  NH/  as  fine  (small)  particle  species;  and  NO3", 
Na^,  K^,  Ca^^,  and  Mg^^  as  course  (large)  particle  species.  As  stated  previously,  particle 
species  deposition  behavior  is  mostly  dominated  by  aerodynamic  resistance  (Ra)  based  on 
aerodynamic  diameter  of  a  particle  and  independent  of  individual  species.  Boundary 
layer  resistance  (Rb)  is  treated  as  being  constant  within  a  season  and  a  LUC.  Surface 
(canopy)  resistance  (Rc)  is  treated  as  being  negligible  for  particulate  ions  and  is  always 
zero.  Therefore,  deposition  velocity  for  particles  is  treated  as  the  inverse  of  the  sum  of  Ra 
and  Rb  (Appendix  B).  Since  it  is  not  informative  to  examine  Ra  or  Rb  independently, 
deposition  velocity  for  particles  for  the  different  LUCs  was  compared: 

•  Deciduous  forest  cover  versus  coniferous  forest  cover  (base  case)  -  In  sunamer, 
monthly  average  deposition  velocities  for  small  and  large  particles  for  a 
deciduous  forest  cover  are  higher  than  those  for  a  coniferous  forest  cover  (base 
case)  as  shown  Figure  33. 

•  Wetland/swamp  versus  coniferous  forest  cover  -  In  Figure  34,  it  can  be  seen  that 
monthly  average  small  and  large  particle  deposition  velocities  are  lower  for 
wetland/swamp  compared  to  a  coniferous  forest  cover  (base  case). 

•  Grassland  versus  coniferous  forest  cover  -  In  Figure  35,  small  and  large  particle 
deposition  velocities  are  even  lower  than  those  for  wetland/swamp. 

•  Cropland  versus  coniferous  forest  cover  -  In  Figure  36,  a  summer  peak  remains  at 
the  same  level  as  that  for  grassland,  however  spring  and  autumn  deposition 
velocities  are  clearly  lower. 

•  Urban  surface  versus  coniferous  forest  cover  -  In  Figure  37,  a  summer  peak  is 
similar  and  higher  compared  to  cropland  and  the  spring  and  autumn  deposition 
velocities  are  lower. 

•  Water  versus  coniferous  forest  cover  is  shown  in  Figure  38. 

•  Snow/ice  versus  coniferous  forest  cover  -  In  Figure  39,  deposition  velocities  are 
low  and  similar  to  those  of  a  wetland/swamp,  grassland,  or  cropland  surface. 

Tables  summarizing  all  the  monthly  values  used  for  Figures  33  to  39  for  the  eight  LUCs 
are  included  in  Appendix  D. 
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Figure  33  Comparison  of  monthly  average  deposition  velocity  of  small  and  large 
particles  for  a  deciduous  forest  and  coniferous  forest  (base  case)  cover 
using  the  AENV  model. 
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Figure  34 


Comparison  of  monthly  average  deposition  velocity  of  small  and  large 
particles  for  a  wetland/swamp  and  coniferous  forest  (base  case)  cover 
using  the  AENV  model. 
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Figure  35      Comparison  of  monthly  average  deposition  velocity  of  small  and  large 
particles  for  grassland  and  coniferous  forest  (base  case)  cover  using 
the  AENV  model. 
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Figure  36 


Comparison  of  monthly  average  deposition  velocity  of  small  and  large 
particles  for  cropland  and  coniferous  forest  (base  case)  cover  using 
the  AENV  model. 
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Figure  37       Comparison  of  monthly  average  deposition  velocity  of  small  and  large 
particles  for  urban  and  coniferous  forest  (base  case)  cover  using  the 
AENV  model. 
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Figure  38      Comparison  of  monthly  average  deposition  velocity  of  small  and  large 
particles  for  water  and  coniferous  forest  (base  case)  cover  using  the 
AENV  model. 
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Figure  39  Comparison  of  monthly  average  deposition  velocity  of  small  and  large 
particles  for  snow/ice  and  coniferous  forest  (base  case)  cover  using  the 
AENV  model. 
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5.0 


MODEL  EVALUATION 


5.1 


Resistance  Factors 


5.1.1  Evaluation  Methods 

Deposition  velocity  (Vd)  is  calculated  as  the  inverse  of  the  sum  of  three  resistances  (Ra  + 
Rb  +  Rc)  on  an  hourly  basis.  In  terms  of  potential  acid  deposition  impacts,  greater 
concern  exists  for  species  with  higher  deposition  velocities.  Because  resistances  are 
added  before  the  sum  is  inverted,  resistance  factors  with  the  largest  magnitude  of  the 
three  will  dominate  determination  of  deposition  velocity.  Therefore,  influence  of  the 
largest  resistance  factor  is  the  most  important  in  terms  of  deposition  velocity. 

Two  methods  were  used  to  evaluate  behavior  of  the  three  resistance  factors  derived  in  the 
AENV  model.  In  the  first  method,  total  resistance  results  (s/m)  for  each  hour  for  the 
three  most  important  depositing  acidic  species  (SO2,  NO2,  and  HNO3)  derived  from  both 
models  were  divided  into  the  following  ranges:  i)  less  than  10,  ii)  10  to  100,  iii)  100  to 
300,  iv)  300  to  1000,  and  v)  >1000  s/m.  Within  each  range,  the  number  of  times  that 
each  of  the  three  resistances  was  largest  (or  dominant)  was  counted.  For  HNO3,  size 
ranges  less  than  100  s/m  were  further  broken  down  and  similarly  analyzed.  These  counts 
are  summarized  in  Appendix  F. 

In  the  second  method,  a  sensitivity  analysis  was  conducted  on  the  AENV  model  by 
separately  increasing  and  decreasing  each  of  the  three  resistances  by  a  magnitude  of  10% 
and  50%  and  observing  the  resulting  changes  in  monthly  summed  deposition  flux  for  1 1 
species. 

5.1.2  Results 


5.1.2.1       Dominant  (Influential)  Resistance  Factors 

The  number  of  hourly  occurrences  in  which  each  of  the  individual  resistances  was  the 
dominant  resistance  was  counted  for  SO2,  NO2,  and  HNO3  for  each  model  and  complete 
results  are  presented  in  Appendix  F.  For  SO2  in  the  ENVC  model,  Rc  dominated  in  every 
instance;  while  in  the  AENV  model,  Rc  usually  dominated,  but  there  were  instances  when 
Ra  and  Rb  dominated. 


For  example,  for  SO2  deposition  in  October,  there  were  96  occurrences  in  which 
total  resistance  was  between  10  and  100  s/m  (Appendix  F).  Of  these  96  occurrences,  Ra 


SO2 
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was  dominant  20  times,  Rb  was  dominant  in  28,  and  Rc  was  dominant  in  48.  In  contrast, 
for  the  ENVC  model  for  this  month,  there  were  only  12  occurrences  in  which  total 
resistance  was  between  10  and  100  s/m  and  all  of  which  were  Rc  dominated  (Appendix 
F).  This  difference  suggests  that  AENV's  approach  for  representing  Rc  values  has  a 
tendency  to  be  less  conservative  relative  to  the  ENVC  approach  for  October  conditions. 

In  January,  there  were  no  occurrences  in  which  total  resistance  was  less  than  300  s/m  for 
the  AENV  model  (Appendix  F),  but  there  were  7 1  occurrences  in  the  ENVC  model  in 
which  total  resistance  was  less  than  300  s/m  (three  of  which  were  between  10  s/m  and 
100).  This  difference  suggests  that  AENV's  approach  for  representing  Rc  values  has  a 
tendency  to  be  more  conservative  realtive  to  the  ENVC  approach  for  January  conditions 
and  in  the  winter  months  (December  through  February).  Rb  was  observed  to  be  dominant 
in  several  instances  in  summer  and  autumn  months  in  the  AENV  model. 

NO2 

For  NO2,  Rc  dominated  in  every  instance  for  the  ENVC  model  for  all  ranges  presented; 
while  in  the  AENV  model,  Rc  was  observed  to  usually  dominate,  but  there  were  instances 
where  Ra  dominated  (Appendix  F).  This  difference  suggests  that  AENV's  model  for 
representing  Rc  values  for  NO2  has  a  tendency  to  be  less  conservative  compared  to  the 
ENVC  approach.  In  spring  and  summer  months,  there  was  a  greater  occurrence  of 
smaller  total  resistance  values  in  the  AENV  model  than  in  the  ENVC  model  (Appendix 
F).  This  difference  suggests  that  again  AENV's  model  for  representing  Rc  values  for 
NO2  has  a  tendency  to  be  less  conservative  compared  to  the  ENVC  approach  spring  and 
summer  months. 

HNO3 

For  HNO3,  Rc  dominated  in  most  instances  in  the  ENVC  model,  but  there  were  instances 
in  which  Ra  or  Rb  was  more  dominant  (Appendix  F).  In  the  AENV  model,  Rb  dominated 
most  often,  followed  by  Ra  and  then  Rc  (Appendix  F).  Recall  that  Rc  for  HNO3  is  take  as 
a  constant  equal  to  10  s/m  for  all  seasons  and  all  surfaces. 

The  lowest  total  resistance  range  (10  to  100  s/m)  was  separated  into  smaller  sections 
(shown  in  Appendix  F)  to  allow  for  a  more  detailed  comparison  of  the  two  models.  In 
the  15  to  25  s/m  range,  Rc  was  always  the  dominant  resistance  for  both  models,  but  there 
were  far  more  occurrences  when  this  was  the  case  for  the  AENV  model  than  for  the 
ENVC  model. 

In  the  25  to  100  s/m  range,  Rb  was  the  most  common  resistance  factor  that  dominated  in 
the  AENV  model;  while  Rc  was  the  most  common  resistance  factor  that  dominated  in  the 
ENVC  model.  These  differences  suggest  that  AENV's  model  for  representing  Rc  values 
for  HNO3  has  a  tendency  to  be  less  conservative  compared  to  the  ENVC  approach  spring 
and  summer  months. 
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5.1.2.2       Sensitivity  Testing  of  Resistance  Factors 

The  other  method  to  examine  the  relative  importance  of  the  three  resistance  factors  was 
to  check  the  influence  of  increasing  and  decreasing  of  magnitude  of  each  of  the 
resistances  on  monthly  summed  deposition  flux. 

Percent  Change  in  Ra 

The  influence  of  a  ten  and  fifty  percent  change  (increase  and  decrease)  in  Ra  on  monthly 
summed  deposition  flux  (kg  H^/ha)  is  shown  in  Figure  40.  Months  that  were  most 
influenced  by  the  change  in  Ra  were  January,  April,  October,  and  December.  The  most 
important  species  that  were  responsible  for  these  changes  were  HNO2  (December)  and 
HNO3  (January,  April,  and  October).  Details  of  deposition  quantities  for  these  two 
species  are  tabulated  in  Appendix  G. 

Percent  Change  in  Rb 

The  influence  of  a  ten  and  fifty  percent  change  (increase  and  decrease)  in  Rb  on  monthly 
summed  deposition  flux  is  shown  in  Figure  41.  The  months  that  were  most  influenced  by 
the  change  in  Rb  were  January,  February,  March,  July,  and  November. 
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Figure  40      Monthly  summed  deposition  flux  for  11  species  using  AENV  model 
showing  influence  of  step  changes  in  Ra  for  a  coniferous  forest  cover. 
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A  change  in  Rb  values  has  more  impact  on  total  deposition  than  does  an  equivalent 
percent  change  in  Ra,  but  not  necessarily  in  the  same  months.  The  most  important 
species  that  were  responsible  for  changes  were  HNO3  (January  and  October),  S04^' 
(March),  NH4^  (March),  and  Ca^^  (July  and  November).  The  changes  observed  in 
February  were  caused  by  an  absence  of  the  neutralizing  effect  of  calcium.  Details  of 
deposition  quantities  for  these  species  are  tabulated  in  Appendix  G. 

Percent  Change  in  Rc 

The  influence  of  a  ten  and  fifty  percent  change  (increase  and  decrease)  in  Rc  on  monthly 
summed  deposition  flux  is  shown  in  Figure  42.  Figure  42,  when  compared  with  Figures 
40  and  41  shows  that  similar  step  changes  in  the  magnitude  of  Rc  are  much  more 
important  in  influencing  total  deposition  than  step  changes  in  either  Ra  or  Rb.  The 
influence  of  changes  in  Rc  is  independent  of  month  of  the  year  although  April  and  August 
had  the  greatest  changes.  The  most  important  species  that  were  responsible  for  changes 
were  SO2  (April,  July,  and  August)  and  NO2  (April,  August,  and  December).  Details  of 
deposition  quantities  for  these  two  species  are  tabulated  in  Appendix  G. 
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5.1 .2.3       Performance  Indicators  for  AENV  Model 


Results  presented  in  the  previous  sections  suggest  that  AENV  model  performance  is 
closely  dependent  on  deposition  velocity  modeling  for  NO2  and  SO2;  particularly  with 
respect  to  selection  of  surface  (canopy)  resistance  (Rc)  values  used  to  parameterize  NO2 
and  SO2  deposition.  Relative  to  performance  of  the  ENVC  model,  there  are  indications 
that  values  used  for  parameterizing  Rc  in  the  AENV  model  are  less  conservative  for  SO2 
in  winter  and  spring  and  more  conservative  in  summer  and  autumn.  For  NO2,  Rc  values 
used  for  parameterizing  Rc  in  the  AENV  model  are  less  conservative  in  the  autumn  and 
winter  and  more  conservative  in  the  spring  and  summer.  For  HNO3,  deposition  in  the 
AENV  model  is  driven  by  Rb. 

It  was  stated  previously  that  the  ENVC  inference  model  does  not  necessarily  represent 
the  best  set  of  assumptions  and  conditions  to  parameterize  dry  deposition,  hi  addition, 
AENV's  parameterization  of  Rc  originates  from  previously  published  values  that  are 
seasonal-  and  species-specific  and,  therefore,  can  be  considered  reasonable.  Thus  while 
there  are  apparent  differences  in  how  both  models  parameterize  Rc,  the  evaluation  results 
presented  above  only  show  where  the  models  depart  in  their  treatment  of  surface 
(canopy)  resistance. 
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5.2 


Parameter  Boundary  Conditions 


One  of  the  recommendations  put  forward  by  WBK  (2008)  was  that  performance  of  the 
AENV  model  may  be  improved  relative  to  performance  of  the  ENVC  model  by  adopting 
some  of  the  boundary  conditions  used  for  specific  parameters  in  the  ENVC  model  and 
including  them  in  the  AENV  model.  Specific  boundary  conditions  in  the  ENVC  model 
that  are  currently  not  used  in  the  AENV  model  include: 

•  A  defined  range  limit  for  Ra  (minimum  of  5;  maximum  of  1000  s/m  (2000  s/m 
over  water). 

•  A  defined  range  limit  for  the  Monin-Obukhov  length  (L)  used  as  a  measure  of 
instabihty  (L  =  -5  if  -5  <  L  <  0  and  L  =  5  if  0  <  L  <  5;  as  specified  in  the  ENVC 
model,  Appendix  A). 

•  Establishing  a  minimum  wind  speed  of  1  m/s  for  wind  speeds  recorded  below  this 
level. 


5.2.1  Evaluation  Methods 

The  AENV  model  was  run  several  times,  each  time  with  a  modification  to  include  a 
separate  boundary  condition  used  by  the  ENVC  model  using  a  coniferous  forest  cover 
and  hourly  pollutant  concentration  and  meteorological  data  from  the  Fort  MacKay  air 
monitoring  station.  Results  were  compared  to  original  AENV  model  results  (base  case). 

5.2.2  Results 


5.2.2.1       Aerodynamic  Resistance 
Minimum  of  5  s/m  for  Ra 

In  the  ENVC  model,  aerodynamic  resistance  (Ra)  is  limited  to  a  minimum  of  5  and  a 
maximum  of  1000  s/m  for  all  land  use  categories  except  for  water;  which  has  an  upper 
limit  of  2000  s/m.  Ra  is  an  indication  of  the  mixing  conditions  of  atmosphere  such  that  a 
very  low  Ra  imphes  maximum  mixing  conditions.  The  purpose  of  a  minimum  limit 
acknowledges  that,  inferentially,  there  is  always  some  limit  to  mixing  possible  in  a 
natural  atmosphere.  Because  Ra  is  added  to  Rb  and  Rc  in  the  determination  of  deposition 
velocity  for  a  pollutant,  a  minimum  Ra  value  will  only  have  an  important  effect  on 
deposition  velocity  if  values  for  Rb  and  Rc  for  that  hour  are  also  small. 

The  number  of  times  hourly  Ra  values  were  within  certain  ranges  within  a  maximum  of 
1000  s/m  were  counted  and  summarized  in  Table  1 1  for  the  AENV  model  modified  to 
run  Ra  with  an  allowable  upper  limit  of  1000  s/m.  Model  conditions  were  for  a 
coniferous  forest  cover  and  using  hourly  pollutant  concentration  and  meteorological  data 
from  the  Fort  MacKay  air  monitoring  station.  From  Table  1 1,  it  can  be  seen  that  May 
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was  the  month  with  the  greatest  number  of  hourly  occurrences  with  the  smallest 
computed  Ra  value  (i.e.,  0  to  10  s/m);  which  implies  that  it  should  be  the  month  with  the 
largest  difference  in  corresponding  deposition  flux. 

The  AENV  model  was  modified  and  rerun  using  an  allowable  lower  value  of      5  s/m 
for  Ra  and  results  (monthly  summed  deposition  flux  for  all  1 1  species)  were  compared  to 
the  model  run  under  base  case  conditions  and  to  ENVC  model  results.  Results  are 
presented  in  Table  12.  As  shown  in  Table  12,  May  did  in  fact  show  the  largest  change  of 
0.06  g  H^/ha  in  deposition  flux,  although  this  difference  only  represented  less  than  a 
0.5%  change  in  total  deposition.  This  indicates  that  the  overall  effect  of  setting  this 
minimum  limit  on  Ra  in  the  AENV  model  is  anticipated  to  be  small. 


Table  11        Counts  of  hourly  Ra  values  within  a  maximum  of  1000  s/m  for  the 
AENV  model  run  with  an  allowable  upper  limit  of  1000  s/m  for  a 
coniferous  forest  cover. 


Ra  (s/m) 

<  10 

<20 

<30 

<40 

<50 

<60 

<70 

<80 

<90 

<  100 

Jan 

58 

105 

58 

28 

16 

9 

7 

4 

8 

8 

Feb 

72 

150 

55 

18 

11 

8 

4 

2 

6 

2 

Mar 

209 

113 

35 

10 

8 

4 

5 

3 

2 

2 

Apr 

102 

221 

82 

41 

25 

15 

0 

6 

2 

2 

May 

293 

91 

22 

12 

15 

9 

5 

6 

1 

Jun 

212 

167 

52 

14 

12 

5 

9 

5 

5 

3 

Jul 

164 

199 

70 

31 

12 

9 

12 

3 

3 

2 

Aug 

140 

210 

70 

18 

14 

16 

4 

3 

3 

7 

Sep 

260 

128 

31 

21 

8 

6 

5 

5 

3 

3 

Oct 

139 

135 

45 

37 

19 

15 

13 

5 

5 

2 

Nov 

29 

182 

92 

55 

20 

18 

11 

13 

5 

7 

Dec 

7 

89 

142 

94 

35 

23 

22 

9 

12 

1 1 

Ra  (s/m) 

<200 

<300 

<400 

<500 

<600 

<700 

<800 

<900 

<  1000 

Jan 

26 

10 

9 

6 

6 

5 

5 

3 

372 

Feb 

11 

10 

7 

6 

1 

1 

1 

1 

330 

Mar 

12 

5 

8 

2 

2 

2 

1 

6 

291 

Apr 

16 

7 

4 

3 

2 

1 

0 

0 

191 

May 

19 

12 

4 

4 

11 

2 

1 

2 

233 

Jun 

25 

13 

6 

6 

5 

4 

1 

6 

170 

Jul 

15 

11 

8 

6 

3 

2 

4 

3 

187 

Aug 

21 

8 

4 

6 

2 

5 

4 

2 

207 

Sep 

22 

13 

9 

2 

11 

2 

2 

2 

187 

Oct 

31 

19 

14 

5 

5 

7 

5 

3 

240 

Nov 

29 

14 

7 

5 

3 

3 

2 

5 

220 

Dec 

54 

16 

10 

9 

15 

2 

0 

2 

192 

Note:  Ranges  presented  include  the  following  increments  -  <10  =  0  to  10  s/m,  <20  =  10  to  20  s/m.  etc. 


Maximum  of  1000  s/m  for  Ra 

Just  as  a  very  low  Ra  value  would  not  necessarily  offer  a  reasonable  representation  of 
natural  conditions  during  daytime;  neither  would  an  extremely  high  Ra  value.  A  very 
high  Ra  value  would  be  associated  with  no  mixing  in  the  atmosphere  (ideally  stable  night 
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time  conditions).  From  a  practical  point  of  view,  very  large  aerodynamic  resistances  are 
not  much  different  than  moderately  large  aerodynamic  resistances  since  deposition 
velocity  involves  the  inverse  of  the  sum  of  all  three  resistances  (Ra  +  Rb  +  Rc)-  When 
hourly  deposition  values  for  a  month  are  summed,  it  is  the  hourly  conditions  with  small 
summed  values  for  the  three  resistance  factors  (and  hence  correspondingly  large 
deposition  velocity  and  deposition  flux  values)  that  are  important  and  drive  the  overall 
monthly  summed  deposition  flux.  Hence  very  large  or  moderately  large  aerodynamic 
resistance  values  (and  correspondingly  very  small  or  moderately  small  deposition 
velocity  and  deposition  flux  values)  are  unimportant. 


Table  12  Comparison  of  monthly  summed  deposition  flux  (g  H'^/ha)  for  AENV 
model  with  minimum  Ra  limit  of  5  s/m  to  AENV  base  case  model  and 
ENVC  model  results  for  coniferous  forest  cover. 


AENV  w/ 

AENV  w/out 

difference  w/limit 

closer  to  ENVC 

limit 

Hmit 

ENVC 

in  AENV  model 

model 

Month  1 

15.5 

15.5 

15.5 

-0.02 

+ 

Month  2 

10.6 

10.6 

11.4 

0.00 

0 

Month  3 

16.1 

16.1 

16.2 

-0.02 

+ 

Month  4 

15.5 

15.5 

17.0 

0.00 

0 

Month  5 

13.6 

13.7 

12.0 

-0.06 

+ 

Month  6 

7.4 

7.4 

6.5 

-0.02 

+ 

Month  7 

7.5 

7.5 

10.7 

0.00 

0 

Month  8 

13.4 

13.4 

12.9 

0.00 

0 

Month  9 

9.6 

9.6 

12.2 

-0.02 

Month  10 

16.3 

16.3 

18.0 

-0.03 

Month  1 1 

-0.2 

-0.2 

3.6 

0.00 

0 

Month  12 

15.6 

15.6 

16.8 

0.00 

0 

total 

140.9 

141.0 

152.8 

-0.2 

An  upper  limit  on  Ra  might  be  expected  to  have  an  effect  where  a  monthly  summed 
deposition  flux  results  from  the  sum  of  a  large  number  of  small  hourly  deposition  flux 
values  (corresponding  to  large  number  of  large  hourly  deposition  velocity  values  perhaps 
influenced  by  large  Ra  values).  As  can  be  seen  from  Table  1 1,  January  had  a  large 
number  of  high  Ra  values  between  900  and  1000  s/m  (372)  and  few  small  Ra  values  (58 
between  0  and  10  s/m).  This  is  in  contrast  to  the  month  of  September;  which  had  a  small 
number  of  high  Ra  values  between  900  and  1000  s/m  (187)  and  numerous  small  Ra  values 
(290  between  0  and  10  s/m).  It  would  be  expected  that  the  month  of  January  would  show 
a  greater  response  to  the  presence  of  a  maximum  limit  set  to  Ra  compared  with 
September. 

This  is  demonstrated  to  be  the  case  in  Table  13.  The  AENV  model  was  rerun  using  a 
maximum  value  of  1000  s/m  for  Ra  and  results  (monthly  summed  deposition  flux)  for  all 
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1 1  species)  were  compared  to  the  model  run  under  the  base  case  conditions  and  to  ENVC 
model  results.  January  shows  the  greatest  change  relative  to  the  base  case  conditions.  In 
addition,  the  direction  of  change  caused  by  using  a  maximum  value  of  1000  s/m  for  Ra 
was  to  bring  AENV  model  results  closer  to  ENVC  model  results.  In  comparison  with  the 
AENV  base  case,  using  a  maximum  value  of  1000  s/m  for  Ra  in  the  AENV  model 
resulted  in  an  overall  increase  of  7  g  H^/ha  in  total  annual  deposition  flux.  This 
represents  a  5%  difference.  Monthly  deposition  is  shown  graphically  in  Figure  43. 


Table  13        Comparison  of  monthly  summed  deposition  flux  (g  H^/ha)  for  AENV 
model  with  maximum  Ra  limit  of  1000  s/m  to  AENV  base  case 
conditions  and  ENVC  model  results  for  coniferous  forest  cover. 


AENV  w/ 

AENV  w/out 

difference  w/limit 

closer  to  ENVC 

limit 

limit 

ENVC 

in  AENV  model 

model 

Month  1 

17.0 

15.5 

15.5 

1.5 

Month  2 

11.3 

10.6 

11.4 

0.7 

+ 

Month  3 

17.1 

16.1 

16.2 

1.0 

Month  4 

16.0 

15.5 

17.0 

0.5 

+ 

Month  5 

14.3 

13.7 

12.0 

0.6 

Month  6 

7.7 

7.4 

6.5 

0.2 

Month  7 

7.5 

7.5 

10.7 

0.1 

+ 

Month  8 

13.7 

13.4 

12.9 

0.3 

Month  9 

9.9 

9.6 

12.2 

0.3 

+ 

Month  10 

17.0 

16.3 

18.0 

0.7 

+ 

Month  1 1 

0.3 

-0.2 

3.6 

0.5 

+ 

Month  12 

16.5 

15.6 

16.8 

0.9 

+ 

total 

148.3 

141.0 

152.8 

7.2 

+ 

5.2.2.2       Monin-Obukhov  Length 

In  the  ENVC  model,  the  Monin-Obukhov  length  (L)  is  limited  to  a  maximum  and 
minimum  value  of  -5  and  5  for  negative  and  positive  values,  respectively.  L  is  used  in 
the  calculation  of  the  stability  function  which  is  ultimately  used  in  calculation  of  Ra.  L 
represents  the  height  above  the  surface  in  which  convectively  driven  turbulence 
dominates  over  mechanically  driven  turbulence  and  it  is  indirectly  a  measure  of  the 
convective  instability  generated  by  the  vertical  heat  flux  through  the  surface  layer.  L  is 
only  meaningful  in  daytime  convectively  driven  boundary  layers. 
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Figure  43      Monthly  summed  deposition  flux  for  11  species  using  AENV  model 
showing  changes  with  and  without  (base  case)  minimum  5  s/m  and 
maximum  1000  s/m  limits  to  hourly  Ra  values  for  a  coniferous  forest 
cover. 


A  negative  value  for  L  represents  atmospheric  instability  and  a  positive  value  represents 
stability.  Small  absolute  values  for  L  represent  extremes  which  are  not  representative  in 
nature.  In  general,  as  L  approaches  zero,  Ra  becomes  larger  and  corresponding 
deposition  velocity  is  smaller.  The  result  of  applying  the  limit  to  L  is  a  slightly  lower  Ra 
and  a  slightly  higher  corresponding  deposition  velocity.  Since  it  is  higher  Ra  values  that 
are  affected,  the  nature  of  the  change  should  be  similar  to  that  of  setting  an  upper  limit  on 
Ra.  The  effects  of  using  limits  to  L  are  shown  in  Table  14  and  Figure  44. 

Annual  summed  deposition  flux  was  increased  by  5  g  H^/ha  which  represents  a  4% 
difference  in  total  deposition.  The  direction  of  change  caused  by  using  limits  to  L  was  to 
bring  AENV  model  results  closer  to  ENVC  model  results  (the  same  general  effect  caused 
by  setting  an  upper  hmit  or  1000  s/m  to  Ra). 


5.2.2.3       Minimum  Wind  Speed 

Wind  speed  values  in  the  AENV  model  are  allowed  to  vary  to  a  minimum  of  0.09  m/s, 
below  which  they  are  assigned  a  value  of  0  m/s.  In  the  ENVC  model,  wind  speed  values 
are  allowed  to  vary  to  a  minimum  of  1  m/s,  below  which  they  are  assigned  a  value  of  1 
m/s.  The  effects  of  using  a  1  m/s  hmit  on  wind  speed  in  the  ENVC  model  is  shown  in 
Table  15  and  Figure  45.  Annual  deposition  flux  was  increased  by  6  g  H^/ha  which 
represents  a  4%  difference  in  total  deposition.  The  direction  of  change  caused  by  using  a 
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1  m/s  limit  to  wind  speed  was  to  bring  AENV  model  results  closer  to  ENVC  model 
results. 


Table  14        Monthly  summed  deposition  flux  (g  H^/ha)  for  1 1  species  showing 
changes  with  and  without  limits  set  on  L  with  AENV  base  case 
conditions  and  ENVC  model  results  for  a  coniferous  forest  cover. 


AENV  w 
limit 

AENV  w/out 
limit 

rilN  V  L 

difference  w/limit 
in  AHNV  model 

closer  to  ENVC 
model 

Month  1 

16.4 

15.5 

15.5 

0.9 

Month  2 

1 1 .0 

10.6 

1 1 .4 

+ 

Month  3 

16.6 

16.1 

16.2 

0.4 

Month  4 

15.9 

15.5 

17.0 

0.4 

+ 

Month  5 

14.4 

13.7 

12.0 

0.7 

Month  6 

7.7 

7.4 

6.5 

0.3 

Month  7 

7.6 

7.5 

10.7 

O.I 

+ 

Month  8 

13.6 

13.4 

12.9 

0.2 

Month  9 

9.8 

9.6 

12.2 

0.2 

+ 

Month  10 

16.9 

16.3 

18.0 

0.6 

+ 

Month  1 1 

0.3 

-0.2 

3.6 

0.5 

+ 

Month  12 

16.1 

15.6 

16.8 

0.5 

+ 

total 

146.2 

141.0 

152.8 

5.1 

+ 

0.020 


"  0.015 

X 
D 


0.010 


o.ooo 


-0.005 
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Figure  44      Monthly  summed  deposition  flux  for  11  species  using  AENV  model 
showing  influence  with  and  without  (base  case)  limits  on  L  for  a 
coniferous  forest  cover. 
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Table  15        Comparison  of  AENV  total  deposition  for  1  m/s  minimum  limit  on 

wind  speed  with  AENV  base  case  conditions  and  ENVC  model  results 
for  a  coniferous  forest  cover. 


AFMV  with 

/AJJiiN  V  Willi 

1  m/s  wind  speed 
limit 

AENV  without 
limit 

ENVC 

UliiCl  CULC 

w/limit 
in  AENV  model 

Closer  LO 
ENVC 
model 

Month  1 

17.9 

15.5 

15.5 

2.4 

0 

Month  2 

11.2 

10.6 

11.4 

0.6 

+ 

Month  3 

16.6 

16.1 

16.2 

0.5 

- 

Month  4 

15.6 

15.5 

17.0 

0.1 

+ 

Month  5 

13.8 

13.7 

12.0 

0.2 

Month  6 

7.5 

7.4 

6.5 

0.1 

Month  7 

7.6 

7.5 

10.7 

0.1 

+ 

Month  8 

13.7 

13.4 

12.9 

0.3 

Month  9 

10.0 

9.6 

12.2 

0.4 

+ 

Month  10 

16.9 

16.3 

18.0 

0.6 

+ 

Month  1 1 

-0.3 

-0.2 

3.6 

-0.1 

Month  12 

16.5 

15.6 

16.8 

0.9 

+ 

total 

147.3 

141.0 

152.8 

6.2 

+ 

Wind  speed  is  involved  in  the  calculation  of  both  Ra  (for  gases  and  particulates)  and  Rb 
(for  gases).  Using  a  minimum  wind  speed  limit  would  be  expected  to  affect  months  that 
have  a  large  number  of  low  hourly  wind  speeds,  resulting  in  low  Ra  and  Rb.  Similar  to  a 
maximum  limit  set  for  Ra,  as  discussed  previously,  January  would  be  expected  to  show 
the  largest  difference.  Table  15  indicates  that  January  shows  the  greatest  change  relative 
to  the  base  case  conditions. 


5.3  Non-Boundary  Condition  Components 

Three  non-boundary  issues  were  examined  during  this  study.  In  the  AENV  model,  two 
issues  were  further  analyzed:  i)  calculation  of  surface  roughness  length  (Zq),  and  ii)  use  of 
relative  humidity  (RH)  as  a  surrogate  parameter  to  define  wetness  conditions  of  a  surface 
instead  of  a  surface  wetness  measurement.  In  the  ENVC  model,  treatment  of  the  particle 
distribution  for  particle  deposition  was  further  analyzed. 
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Figure  45      Monthly  summed  deposition  flux  for  11  species  using  AENV  model 
showing  the  influence  of  with  and  without  (base  case)  a  1  m/s 
minimum  limit  on  wind  speed  for  a  coniferous  forest  cover. 


5.3.1        Evaluation  Methods 


5.3.1.1       Surface  Roughness  Length 

One  potential  limitation  of  the  AENV  method  is  a  requirement  to  calculate  surface 
roughness  length  (Zo)  as  a  monthly  average  using  hourly  data  where  the  wind  speed  is  >6 
m/s.  At  issue  is  a  situation  where  only  a  small  number  of  hourly  values  exist  during  the 
month  in  which  the  wind  speed  is  >6  m/s.  This  results  in  greater  uncertainty  in  the 
calculated  average  Zo  value  when  computed  from  a  small  number  of  hourly  values  in 
which  the  wind  speed  is  >6  m/s  during  the  month. 

The  AENV  model  was  rerun  using  modified  criteria  to  establish  representative  values  for 
Zo  based  on  a  minimum  requirement  that  at  least  seven  hourly  data  points  must  be  used  in 
the  calculation  of  an  average  zq  for  that  month: 

•  If  one  consecutive  month  did  not  have  at  least  seven  hourly  data  points  to  be  used 
in  the  calculation  of  the  Zo  (e.g.,  November  data  in  this  case);  then  an  average 
monthly  Zo  computed  from  for  the  two  adjacent  months  (October  and  December 
in  this  case)  are  to  be  used. 

•  If  two  or  more  consecutive  months  did  not  have  at  least  seven  hourly  data  points 
to  be  used  in  calculation  of  Zo  (e.g.,  February  and  March  in  this  case);  then  an 
average  annual  value  computed  from  Zo  values  for  all  months  that  had  at  least 
seven  hourly  data  points  was  used  to  represent  these  months.  This  could  be 
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applicable  to  surface  cover  characteristics  that  are  not  expected  to  undergo 
substantial  change  from  season  to  season  (e.g.,  a  coniferous  forest  canopy). 
•    Resulting  deposition  flux  values  were  then  compared  with  original  results. 

Changing  the  minimum  wind  speed  criterion  (>6  m/s)  in  which  to  calculate  a  monthly  Zo 
value  was  also  examined  by  evaluating  the  relationship  between  different  minimum  wind 
speed  criteria  above  and  below  6  m/s  and  corresponding  calculated  monthly  Zo  values. 

5.3.1 .2  Surface  Wetness  and  Relative  Humidity 

Wetness  alters  the  characteristics  of  a  surface  to  act  a  sink  for  depositing  species 
(primarily  gases).  The  AENV  model  requires  an  hourly  measurement  to  represent  the 
amount  of  wetness  on  a  surface  in  determination  of  surface  resistance  (Rc).  The  model  is 
set  up  to  use  field  measurement  data  from  surface  wetness  sensing  equipment.  However, 
past  experience  has  indicated  that  this  field  equipment  does  not  provide  reliable 
measurement  data  in  an  Alberta  setting  due  to  the  wide  temperature  ranges  experienced. 
Alternatively,  the  AENV  model  can  use  relative  humidity  (RH)  as  a  surrogate  criterion 
for  surface  wetness  (Bates,  1996).  The  AENV  method  uses  a  lower  limit  for  RH  of 
87.4%  in  classifying  a  surface  as  being  wet  for  the  full  hour  after  Bates  (1996). 

The  effect  of  using  different  limits  for  RH  was  examined.  The  original  AENV  model 
was  modified  to  set  the  RH  value  as  80%  for  the  cutoff  in  determination  of  a  wet/dry 
condition  for  selection  of  the  Rc  value  for  SO2  and  NO2  for  each  of  land  use  category. 
Similarly,  the  original  AENV  model  was  modified  to  use  a  RH  value  of  92%  as  the 
cutoff.  Results  in  terms  of  monthly  summed  deposition  flux  for  all  species  were 
compared  to  original  results  for  each  land  use  category. 

5.3.1 .3  Settling  Properties  of  Large  Particle  Size  Distributions 

The  ENVC  model  has  an  inherent  hmitation  in  the  way  that  it  treats  settUng  properties  of 
species  associated  with  large  particles  (i.e.,  NO3",  Na^,  K^,  Ca^^,  and  Mg'^"^).  The  model 
forces  the  treatment  settling  properties  for  >24%  of  these  particles  (by  mass)  as  being 
similar  to  settling  properties  of  fine  particles.  This  results  in  an  underestimate  of  the 
neutralization  capacity  offered  by  depositing  particulates  containing  base  cations;  and  an 
overestimate  in  total  deposition  flux  to  a  surface  when  all  depositing  species  are  taken 
into  account.  This  Hmitation  was  further  investigated  to  show  how  large  this  effect  was 
for  the  ENVC  model,  hi  essence,  because  the  ENVC  model  tends  to  estimate  higher  total 
deposition  fluxes  compared  to  the  AENV  model,  better  agreement  among  these  two 
models  would  be  achieved. 
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5.X2  Results 


5.3.2.1       Surface  Roughness  Length 

One  potential  issue  with  the  AENV  model  is  the  way  surface  roughness  length  (Zo)  is 
calculated.  The  AENV  model  requires  that  hourly  data  to  be  used  in  calculation  of  Zo 
during  a  particular  month  be  initially  based  on  a  monthly  average  Zo  value  using  hourly 
data  where  the  wind  speed  is  >6  m/s  in  that  month  (WBK,  2008).  For  the  2003 
meteorological  data  used  in  this  study,  the  month  of  February  did  not  have  a  single  hourly 
wind  speed  value  >6  m/s  (shown  in  Table  16)  such  that  it  was  necessary  to  use  an 
average  Zo  value  of  the  two  surrounding  months.  In  general,  there  was  a  lack  of  hourly 
values  for  several  months  (e.g.,  February,  March,  and  November)  for  use  in  calculation  of 
Zo  given  the  6  m/s  minimum  wind  speed  cutoff  value.  The  number  of  hourly  values  used 
for  calculating  each  of  the  monthly  average  zo  values  ranged  from  0  to  42  (shown  in 
Table  16). 

Table  16        Monthly  average  surface  roughness  length  (Zo,  m),  standard  deviation, 
and  sample  size  (number  of  hourly  values)  in  which  wind  speed  >6 
m/s  for  calculation  of  monthly  average  surface  roughness  length  for  a 
coniferous  forest  cover. 


Month 

Zo 

Standard  Deviation 

Sample 
Size 

January 

1.24 

0.40 

19 

February 

0 

March 

1.50 

0.54 

6 

April 

0.67 

0.28 

42 

May 

1.46 

0.75 

29 

June 

1.26 

0.45 

26 

July 

1.20 

0.49 

10 

August 

1.11 

0.52 

13 

September 

1.61 

0.89 

9 

October 

1.30 

0.64 

10 

November 

0.72 

0.24 

2 

December 

0.50 

0.16 

7 

1.12 

Weighted  Annual  Average 

For  months  with  too  few  hours  in  which  to  calculate  Zo,  it  would  perhaps  be  more  suitable 
if  monthly  wind  speed  data  from  another  meteorological  year  could  be  used  to 
supplement  the  number  of  observations  used  in  the  initial  calculation  of  a  Zo  value  for  that 
that  month  at  a  site-specific  location.  Unfortunately,  these  data  would  not  be  applicable 
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for  a  location  undergoing  changes  in  surrounding  land  use.  The  approach  taken  here  was 
to  establish  a  requirement  for  the  minimum  number  of  hourly  observations  in  which  wind 
speed  >6  m/s  for  calculating  zq  for  a  particular  month: 

•  Seven  was  the  number  selected  as  the  minimum  number  of  observations  in  which 
to  use. 

•  The  average  Zo  value  of  two  surrounding  months  was  established  for  any  month  in 
which  an  insufficient  number  of  observations  existed. 

•  If  consecutive  months  data  have  an  insufficient  number  of  observations  for  a 
surface  cover  whose  characteristics  are  not  expected  to  undergo  substantial 
change  from  season  to  season  (e.g.,  a  coniferous  forest  canopy),  then  an  annual 
average  should  be  used  for  these  months.  (Alternatively,  if  consecutive  months 
have  an  insufficient  number  of  observations  for  a  surface  cover  whose 
characteristics  are  not  expected  to  undergo  substantial  change  within  a  season, 
then  an  average  for  the  month  within  the  season  could  be  used  for  to  represent 
these  months.) 

These  changes  were  incorporated  into  the  AENV  model  and  results  presented  in  Figure 
46.  Figure  46  is  a  plot  of  monthly  sunmied  deposition  flux  for  1 1  species  using  the 
AENV  model  showing  the  influence  resulting  from  use  of  surrogate  Zo  values  for  the 
months  of  February,  March,  and  November  versus  the  base  case.  Zo  used  for  November 
was  taken  as  the  average  of  October  and  December  and  the  Zo  value  used  for  February 
and  March  was  taken  as  the  annual  average  value  shown  in  Table  16.  Figure  46  shows 
that  these  substitutions  only  result  in  sUghtly  higher  deposition  fluxes  in  February  and 
March;  while  the  difference  in  November  is  negligible. 

Another  potential  issue  with  calculation  of  Zq  in  the  AENV  model  is  that  there  was  a  wide 
variation  computed  for  different  months  across  the  year.  In  theory,  a  mature  coniferous 
forest  should  have  a  similar  surface  roughness  length  across  the  year  because  the  foliage 
does  not  change  seasonally.  The  ENVC  model  (Appendix  A)  assumes  that  the  surface 
roughness  length  (Zq)  for  an  evergreen  needle  leaf  forest  is  consistent  throughout  the  year 
at  a  value  of  0.9  m.  In  the  AENV  model  in  this  study,  Zo  ranged  from  0.5  m  in  December 
to  1.6  meters  in  October  (Table  16). 

Even  within  a  particular  month,  Zq  was  variable  as  seen  in  standard  deviations  for  the 
AENV  model  (shown  in  Table  16).  It  has  been  observed  by  others  that  values  for  Zo  can 
change  for  different  wind  speeds  if  the  surface  cover  is  such  that  it  can  bend  with  the 
wind,  thereby  reducing  the  roughness  length  (Marht,  2001). 
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Figure  46      Monthly  summed  deposition  flux  for  11  species  using  AENV  model 

showing  influence  resulting  from  corrected  Zo  values  for  the  months  of 
February,  March,  and  November  for  a  coniferous  forest  cover. 


With  respect  to  the  dependence  of  Zq  upon  wind  speed,  it  was  found  that  if  the  value  of 
the  wind  speed  criterion  used  to  calculate  a  monthly  value  for  Zo  was  decreased  from  6 
m/s,  the  monthly  re-calculated  mean  Zq  value  and  corresponding  standard  deviation 
increased.  This  is  shown  in  Table  17. 

Weighted  annual  average  Zo  values  were  computed  and  plotted  in  Figure  47  for  each  wind 
speed  criterion  listed  in  Table  17.  This  figure  shows  a  linear  change  in  weighted  annual 
average  Zo  for  different  wind  speeds  cutoffs  used.  If  wind  speeds  are  further  separated 
into  specific  ranges  and  a  corresponding  Zq  value  calculated  for  each  range,  as  was  done 
in  Table  18,  it  is  seen  that  Zq  values  are  highly  variable  among  the  ranges  shown.  Wind 
direction  has  also  been  found  to  affect  calculation  of  Zo  as  a  result  of  differences  in 
landscape  features  in  different  directions  from  a  monitoring  station  (CEC,  1993). 


5.3.2.2       Surface  Wetness  and  Relative  Humidity 


Surface  wetness  is  difficult  to  measure  and  it  is  very  difficult  to  generalize  data  from  one 
area  to  an  entire  area  (Panofsky,  1974).  Surface  wetness  data  are  required  for  the  AENV 
model;  however  experience  has  shown  that  these  data  are  frequently  unavailable  due  to 
problems  with  surface  wetness  sensing  equipment  in  the  field.  Consequently,  the  AENV 
model  allows  a  minimum  relative  humidity  (RH)  measurement  cutoff  (87%)  to  be  used  as 
an  approximation  to  represent  a  surface  as  being  wet.  The  selection  of  the  RH  value 
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cutoff  was  further  examined  to  assess  the  sensitivity  of  this  AENV  model  criterion.  The 
current  cutoff  in  the  model  is  87%  and  two  alternative  values  -  80%  and  92%  -  were 
assessed. 


Table  17        Monthly  average  re-calculated  surface  roughness  length  (Zo)  values 
(m),  standard  deviation,  and  sample  size  (number  of  hourly  values) 
for  different  wind  speed  criteria  used  in  the  re-calculation  for  a 
coniferous  forest  cover. 


Wind  speed 
greater  than 

January 

February 

March 

April 

May 

June 

or  equal  to 

Mean 

St  Dev 

Sample 

Mean  St  Dev  Sample 

Mean 

St  Dev  Sample 

Mean  St  Dev  Sample 

Mean 

St  Dev 

Sample 

Mean  St  Dev  Sample 

(m/s) 

Size 

Size 

Size 

Size 

Size 

Size 

1.5 

1.76 

1.26 

347 

1.84 

1.27 

350 

1.93 

1.25 

517 

1.89 

1.34 

578 

2.32 

1.26 

601 

2.37 

1.17 

573 

2 

1.49 

1.09 

210 

1.68 

1.21 

215 

1.80 

1.19 

394 

1.77 

1.27 

476 

2.28 

1.19 

506 

2.41 

1.09 

447 

2.5 

1.18 

0.95 

133 

1.59 

1.10 

124 

1.66 

1.15 

287 

1.69 

1.25 

374 

2.20 

1.08 

398 

2.35 

1.02 

354 

3 

1.05 

0.85 

95 

1.66 

1.04 

68 

1.60 

1.11 

211 

1.45 

1.20 

286 

2.10 

1.01 

309 

2.17 

0.96 

259 

3.5 

1.19 

0.84 

64 

1.69 

1.01 

37 

1.48 

1.06 

161 

1.30 

1.13 

227 

2.10 

0.96 

236 

2.11 

0.93 

178 

4 

1.27 

0.74 

44 

1.47 

0.99 

16 

1.25 

0.95 

108 

1.00 

0.91 

169 

2.05 

0.96 

183 

1.99 

0.89 

122 

4.5 

1.42 

0.71 

36 

2.15 

0.94 

6 

1.03 

0.75 

58 

0.85 

0.71 

137 

1.97 

0.92 

125 

1.80 

0.85 

78 

5 

1.45 

0.66 

27 

2.00 

1.39 

2 

0.96 

0.67 

29 

0.70 

0.40 

96 

1.99 

0.93 

84 

1.59 

0.70 

54 

5.5 

1.46 

0.64 

25 

1.26 

0.54 

14 

0.68 

0.30 

66 

1.65 

0.91 

51 

1.45 

0.63 

38 

6 

1.24 

0.40 

19 

1.50 

0.54 

6 

0.67 

0.28 

42 

1.46 

0.75 

30 

1.26 

0.45 

25 

7 

1.20 

0.41 

14 

0.61 

0.23 

16 

0.97 

0.43 

10 

1.03 

0.30 

16 

8 

1.13 

0.14 

7 

0.50 

0.13 

4 

0.75 

0.11 

6 

0.95 

0.22 

10 

9 

1.22 

0.19 

3 

0.77 

0.12 

3 

Wind  speed 
greater  than 

July 

August 

September 

October 

November 

December 

or  equal  to 

Mean  St  Dev  Sample 

Mean 

St  Dev  Sample 

Mean 

St  Dev 

Sample 

Mean 

St  Dev  Sample 

Mean 

St  Dev  Sample 

Mean 

St  Dev  Sample 

(m/s) 

Size 

Size 

Size 

Size 

Size 

Size 

1.5 

2.22 

1.22 

433 

2.36 

1.16 

473 

1.77 

1.12 

437 

1.50 

1.00 

433 

1.29 

1.12 

442 

1.41 

1.35 

503 

2 

2.07 

1.10 

314 

2.31 

1.05 

341 

1.67 

1.05 

348 

1.47 

0.95 

311 

1.07 

0.97 

322 

1.14 

1.17 

325 

2.5 

1.94 

1.04 

226 

2.29 

1.04 

248 

1.59 

0.96 

274 

1.46 

0.90 

222 

0.92 

0.88 

238 

0.97 

0.96 

175 

3 

1.89 

1.03 

157 

2.21 

1.01 

167 

1.59 

0.90 

220 

1.36 

0.84 

155 

0.78 

0.70 

164 

0.98 

0.98 

97 

3.5 

1.74 

0.97 

110 

2.08 

0.98 

112 

1.59 

0.83 

162 

1.37 

0.80 

103 

0.76 

0.69 

111 

1.15 

0.95 

45 

4 

1.64 

0.85 

73 

1.96 

0.97 

60 

1.63 

0.77 

124 

1.43 

0.78 

66 

0.77 

0.59 

56 

1.04 

0.73 

22 

4.5 

1.56 

0.78 

49 

1.83 

1.00 

40 

1.66 

0.71 

85 

1.52 

0.77 

46 

0.91 

0.60 

31 

0.77 

0.47 

14 

5 

1.45 

0.76 

33 

1.76 

0.99 

28 

1.76 

0.72 

41 

1.54 

0.77 

25 

1.34 

0.55 

13 

0.55 

0.17 

9 

5.5 

1.32 

0.57 

19 

1.53 

0.91 

21 

1.84 

0.76 

20 

1.38 

0.69 

16 

0.95 

0.36 

7 

0.52 

0.16 

8 

6 

1.20 

0.49 

10 

1.11 

0.52 

13 

1.61 

0.89 

9 

1.30 

0.64 

10 

0.72 

0.24 

2 

0.50 

0.16 

7 

7 

1.22 

0.52 

5 

1.04 

0.46 

10 

0.40 

0.12 

4 

8 
9 

0.66 

0.20 

3 

Wetness  affects  surface  resistance  (Rc)  for  SO2  and  NO2  in  the  AENV  model.  Default 
values  used  to  represent  Rc  in  the  AENV  model  for  SO2  are  generally  lower  when  a 
surface  is  considered  wet  as  opposed  to  dry  (refer  to  Appendix  B).  Thus  wetness  will 
enhance  the  contribution  of  SO2  deposition  to  a  surface.  On  the  other  hand,  default 
values  used  to  represent  Rc  in  the  AENV  model  for  NO2  are  generally  higher  when  a 
surface  is  considered  wet  as  opposed  to  dry  (Appendix  B).  Thus  wetness  will  decrease 
the  contribution  of  NO2  deposition  to  a  surface. 
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Figure  47      Annual  average  re-calculated  surface  roughness  length  (Zo)  values  (m) 
for  different  wind  speed  criteria  used  in  the  re-calculation  for  a 
coniferous  forest  cover. 
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Table  18       Monthly  calculated  surface  roughness  length  (Zo)  values  (m)  and 

standard  deviation  for  different  wind  speed  ranges  for  a  coniferous 
forest  cover. 


Wind  speed 

January 

February 

March 

April 

May 

June 

range 

(m/s) 

Mean 

St  Dev 

Mean 

St  Dev 

Mean 

St  Dev 

Mean 

St  Dev 

Mean 

St  Dev 

Mean 

St  Dev 

2<u<3 

1.85 

1.14 

1.69 

1.29 

2.02 

1.25 

2.25 

1.24 

2.56 

1.39 

2.73 

1.18 

3  <u  <4 

0.86 

0.90 

1.72 

1.06 

1.97 

1.15 

2.10 

1.26 

2.17 

1.08 

2.34 

0.99 

4<u  <5 

0.99 

0.80 

1.40 

0.97 

1.35 

1.01 

1.38 

1.20 

2.09 

0.99 

2.30 

0.90 

5<u<6 

1.96 

0.88 

2.00 

1.39 

0.82 

0.64 

0.73 

0.48 

2.28 

0.89 

1.87 

0.76 

6<u<7 

1.36 

0.36 

1.50 

0.54 

0.70 

0.30 

1.71 

0.76 

1.67 

0.38 

7<u<8 

1.27 

0.58 

0.64 

0.25 

1.30 

0.54 

1.16 

0.38 

8<u<9 

1.07 

0.07 

0.50 

0.13 

0.74 

0.13 

0.95 

0.22 

9<u 

1.22 

0.19 

0.77 

0.12 

Wind  speed 

July 

August 

September 

October 

November 

December 

range 

(m/s) 

Mean 

St  Dev 

Mean 

St  Dev 

Mean 

St  Dev 

Mean 

St  Dev 

Mean 

St  Dev 

Mean 

St  Dev 

2<u<3 

2.24 

1.15 

2.41 

1.08 

1.79 

1.26 

1.57 

1.03 

1.36 

1.12 

1.21 

1.23 

3  <u  <4 

2.12 

1.11 

2.35 

1.01 

1.54 

1.04 

1.31 

0.88 

0.79 

0.75 

0.97 

1.05 

4<u  <5 

1.79 

0.90 

2.14 

0.92 

1.57 

0.79 

1.36 

0.79 

0.60 

0.49 

1.38 

0.78 

5<u<6 

1.55 

0.83 

2.33 

0.96 

1.81 

0.68 

1.70 

0.83 

1.45 

0.51 

0.71 

0.08 

6<u<7 

1.19 

0.53 

1.32 

0.76 

1.61 

0.89 

1.30 

0.64 

0.72 

0.24 

0.62 

0.13 

7<u<8 

1.22 

0.52 

1.21 

0.45 

0.40 

0.12 

8<u<9 

0.66 

0.20 

9<u 

The  effects  of  changing  the  surface  wetness  cutoff  for  all  the  land  use  categories,  except 
water,  on  monthly  summed  deposition  flux  for  1 1  species  are  shown  in  Figures  48  to  54 
below.  Water  had  no  change  associated  with  it.  These  figures  show  that  there  is  very 
httle  change  in  monthly  summed  deposition  flux  for  1 1  species  caused  using  a  different 
RH  cutoff  to  represent  a  surface  as  being  wet  (e.g.,  80%,  92%  ,  or  the  base  case  of  87%). 
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Figure  48      Monthly  summed  deposition  flux  for  11  species  using  AENV  model 
showing  influence  of  using  a  different  RH  cutoff  for  establishing  a 
surface  as  wet  (base  case  RH  =  0.87)  for  deciduous  forest  cover. 


Month 


RH  cutoff  0.87  (base  case)*-  RH  cutoff  0.80      RH  cutoff  0.92 


Figure  49      Monthly  summed  deposition  flux  for  11  species  using  AENV  model 
showing  influence  of  using  a  different  RH  cutoff  for  establishing  a 
surface  as  wet  (base  case  RH  =  0.87)  for  a  coniferous  forest  cover. 
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Figure  50      Monthly  summed  deposition  flux  for  11  species  using  AENV  model 
showing  influence  of  using  a  different  RH  cutoff  for  establishing  a 
surface  as  wet  (base  case  RH  =  0.87)  for  wetland/swamp  cover. 


Figure  51      Monthly  summed  deposition  flux  for  11  species  using  AENV  model 
showing  influence  of  using  a  different  RH  cutoff  for  establishing  a 
surface  as  wet  (base  case  RH  =  0.87)  for  grassland. 
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Figure  52      Monthly  summed  deposition  flux  for  11  species  using  AENV  model 
showing  influence  of  using  a  different  RH  cutoff  for  establishing  a 
surface  as  wet  (base  case  RH  =  0.87)  for  cropland. 


-—  RH  cutoff  0.87  (base  case)  -*- RH  cutoff  0.80 


RH  cutoff  0.92 


Figure  53      Monthly  summed  deposition  flux  for  11  species  using  AENV  model 
showing  influence  of  using  a  different  RH  cutoff  for  establishing  a 
surface  as  wet  (base  case  RH  =  0.87)  for  urban  cover. 
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Monthly  summed  deposition  flux  for  11  species  using  AENV  model 
showing  influence  of  using  a  different  RH  cutoff  for  establishing  a 
surface  as  wet  (base  case  RH  =  0.87)  for  snow/ice. 


5.3.2.3       Particulate  Size  Distributions 


Differences  exist  between  the  AENV  and  ENVC  models  in  terms  of  how  they  predict 
large  particle  deposition.  In  the  original  design  of  the  ENVC  model  -  which  only 
addressed  S04^'  and  NH4^  deposition  associated  with  small  particles  -  40  aerodynamic 
sizes  (referred  to  as  binsizes  by  ENVC)  were  used  to  represent  an  overall  particle  size 
distribution  in  calculation  of  the  deposition  velocity  for  species  associated  with  small 
particles.  The  largest  size  encompasses  particles  up  to  10  jim  in  diameter;  which  matches 
the  maximum  size  of  particles  classified  as  PMio-  This  number  of  binsizes  (40) 
encompasses  almost  100%  of  a  typical  small  particle  distribution  having  a  geometric 
mean  of  0.5  |am  and  a  standard  deviation  of  2.0  |im. 

When  the  ENVC  model  was  modified  to  include  calculation  of  deposition  velocity  for 
species  associated  with  large  particles  (i.e.,  NO3',  Na^,  K^,  Ca^^,  and  Mg^^),  the  number 
of  binsizes  remained  at  40  so  as  to  hmit  the  largest  particle  size  to  that  detected  as  PMio 
(i.e.,  10  i^m)  for  determining  deposition  flux  of  these  species.  The  size  distribution  of 
sodium  (Na^)  was  used  as  a  conservative  representation  for  the  distribution  of  all  large 
particles  (NO3",  Na^,  K^,  Ca^^,  and  Mg^^)  in  the  ENVC  model.  However  a  typical 
sodium  particle  distribution  has  a  geometric  mean  of  0.5  |im  and  a  standard  deviation  of 
2.0  \im,  and  only  75.86%  of  this  size  distribution  is  encompassed  within  the  first  40 
binsizes  to  a  diameter  of  10  fxm. 

Thus  >24%  of  the  particle  sizes  for  this  distribution  are  greater  than  10  ^m  and  are 
unaccounted  for  in  the  ENVC  model.  These  larger  sized  particles,  >10  |im  in  diameter, 
are  base  cation  species  with  increased  settling  properties  as  a  result  of  their  large  size 
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compared  to  fine  particles.  The  net  results  of  taking  larger  sizes  into  account  in  the 
calculation  would  be  a  greater  overall  buffering  effect;  or  more  importantly  a  net 
decrease  in  acid  deposition  flux  when  all  species  are  considered  in  the  calculation. 
One  possible  way  to  correct  for  this  limitation  is  to  divide  each  binsize  diameter  used  in 
the  ENVC  model  by  0.7586  (75.86%/100%)  to  standardize  and  better  represent  the  more 
likely  settling  properties  of  species  associated  with  larger  particles.  When  this 
adjustment  is  made,  corrected  ENVC  model  results  (Figure  55)  show  much  better 
agreement  with  AENV  model  results  than  previously  presented. 


^  o.o2oa 


1       2       3       4       5       6       7  9       10      11  12 

Month 

AENV        ENVC        corrected  ENVC 


Figure  55  Monthly  summed  deposition  flux  for  11  species  for  AENV  model  base 
case,  ENVC  model,  and  ENVC  model  corrected  for  large  particle  size 
distribution  for  a  coniferous  forest  cover. 


The  net  difference  in  annual  deposition  flux  for  1 1  species  between  the  AENV  model 
(base  case)  and  corrected  ENVC  model  was  1.2%;  as  opposed  to  a  7.7%  difference 
without  a  correction  (Table  2).  As  expected,  the  largest  difference  was  observed  in 
November.  A  breakdown  of  the  influence  of  the  ENVC  model  corrected  for  large 
particle  size  distribution  on  monthly  summed  deposition  flux  for  NO3',  Na^,  K^,  Ca~^,  and 
Mg^"^  for  coniferous  forest  cover  is  provided  in  Appendix  H. 

Increasing  the  number  of  bins  to  capture  larger  particle  sizes  in  the  distribution  (i.e.,  >10 
|im)  is  typically  not  done  because  of  uncertainty  in  extrapolating  the  amount  of  nitrate  in 
particles  larger  than  those  detected  by  PMio  measurements.  Current  use  of  a  truncated 
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distribution  for  large  particles  in  ENVC's  model  produces  a  conservative  estimate  of 
deposition  flux  when  considering  all  species.  Because  a  majority  of  base  cations  is 
associated  with  large  particle  deposition,  underestimating  deposition  of  large  particles 
results  in  an  underestimation  of  the  neutraUzation  capacity  of  base  cations. 
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6.0        FINDINGS  AND  RECOMMENDATIONS 


An  evaluation  of  the  analytical  basis  of  the  Alberta  Environment  (AENV)  and 
Environment  Canada  (ENVC)  inference  models  for  predicting  dry  deposition  of  1 1 
species  was  undertaken.  These  species  included  SO2,  NO2,  HNO3,  HNO2,  S04~ ,  NH/, 
NO3',  Na^,  K^,  Ca"^,  and  Mg~^.  Important  findings  and  recommends  include: 

1.  The  three  most  important  species  contributing  to  potential  acid  input  were  SO2, 
NO2,  and  HNO3. 

2.  Sensitivity  testing  of  resistance  factors  Ra,  Rb,  and  Rc  used  in  AENV's  model 
showed  that  step  changes  in  the  magnitude  of  Rc  were  much  more  important  in 
influencing  total  deposition  compared  to  step  changes  in  either  Ra  or  Rb  for  all 
species  examined.  The  influence  of  changes  in  Rc  was  independent  of  month  of 
the  year. 

3.  AENV  model  performance  in  terms  of  overall  deposition  flux  of  1 1  species  was 
observed  to  be  closely  dependent  on  deposition  velocity  modeling  for  SO2  and 
NO2;  particularly  with  respect  to  selection  of  Rc  values  used  to  parameterize  these 
gases.  Relative  to  performance  of  the  ENVC  model,  values  used  to  represent  Rc 
in  the  AENV  model  were  less  conservative  for  SO2  in  winter  and  spring  and  more 
conservative  in  summer  and  autumn.  For  NO2,  values  used  to  represent  Rc  in  the 
AENV  model  were  less  conservative  in  autumn  and  winter  and  more  conservative 
in  spring  and  summer.  For  HNO3,  deposition  in  the  AENV  model  was  driven  by 
Rb-  These  differences  do  not  necessarily  imply  problems  in  the  AENV  model. 
The  ENVC  inference  model  does  not  necessarily  represent  the  best  set  of 
assumptions  and  conditions  to  parameterize  dry  deposition.  In  addition,  AENV's 
parameterization  of  Rc  originates  from  previously  published  values  that  are 
seasonal-  and  species-specific  and,  therefore,  can  be  considered  reasonable. 

4.  Ra  is  an  indication  of  the  mixing  conditions  in  the  atmosphere  such  that  a  very 
low  Ra  implies  maximum  mixing  conditions  experienced  during  daytime 
conditions.  Setting  a  minimum  value  of  5  s/m  for  Ra  in  the  AENV  model  had 
only  a  small  positive  change  in  model  performance  relative  to  the  ENVC  model. 
Although  the  effect  of  this  change  is  small,  it  is  recommended  for  future  use  in  the 
AENV  model.  The  purpose  of  a  minimum  limit  for  Ra  acknowledges  that, 
inferentially,  there  is  always  some  limit  to  mixing  possible  in  a  natural 
atmosphere. 

5.  Setting  a  maximum  value  of  1000  s/m  for  Ra  in  the  AENV  model  resulted  in 
better  agreement  in  overall  deposition  flux  of  1 1  species  relative  to  the  ENVC 
model  for  the  conditions  tested.  High  Ra  values  would  be  associated  with  no 
mixing  in  the  atmosphere  (ideally  stable  night  time  conditions).  An  extremely 
high  Ra  value  would  not  necessarily  offer  a  reasonable  representation  of  natural 
conditions.  From  a  practical  point  of  view,  very  large  Ra  values  are  not  much 
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different  than  moderately  large  values  since  deposition  velocity  involves  the 
inverse  of  the  sum  of  all  three  resistances  (Ra  +  Rb  +  Rc).  When  hourly  deposition 
values  for  a  month  are  summed,  hourly  conditions  with  small  summed  values  for 
the  three  resistance  factors  (and  hence  correspondingly  large  deposition  velocity 
and  deposition  flux  values)  are  most  important  and  drive  the  overall  monthly 
summed  deposition  flux.  Hence  very  large  or  moderately  large  Ra  values  (and 
correspondingly  very  small  or  moderately  small  deposition  velocity  and 
deposition  fluxes)  are  unimportant.  This  limit  is  recommended  for  future  use  in 
the  AENV  model. 

6.  The  Monin-Obukhov  length  (L)  is  used  in  calculation  of  the  stabihty  function 
which  is  ultimately  used  to  calculate  Ra.  This  length  represents  the  height  above 
the  surface  in  which  convectively  driven  turbulence  dominates  over  mechanically 
driven  turbulence  and  it  is  indirectly  a  measure  of  the  convective  instability 
generated  by  the  vertical  heat  flux  through  the  surface  layer.  It  is  only  meaningful 
in  daytime  convectively  driven  boundary  layers.  L  is  limited  to  a  maximum  and 
minimum  value  of  -5  and  5  m  for  negative  and  positive  values,  respectively,  in  the 
ENVC  model.  A  negative  value  represents  atmospheric  instabiUty  and  a  positive 
value  represents  stability.  Li  general,  as  L  approaches  zero,  Ra  becomes  larger 
and  corresponding  deposition  velocity  is  smaller.  The  result  of  applying  these 
limits  to  L  is  a  slightly  lower  Ra  and  a  slightly  higher  corresponding  deposition 
velocity.  Since  it  is  higher  Ra  values  that  are  affected,  the  nature  of  the  change  is 
similar  to  that  of  setting  an  upper  limit  on  Ra.  The  effect  of  using  these  limits  to  L 
was  to  bring  AENV  model  results  closer  to  ENVC  model  results  (same  general 
effect  caused  by  setting  an  upper  Hmit  to  Ra).  These  limits  are  recommended  for 
future  use  in  the  AENV  model  as  small  absolute  L  values  represent  extremes 
which  are  not  representative  in  nature. 

7.  Wind  speed  values  in  the  AENV  model  are  allowed  to  vary  to  a  minimum  of  0.09 
m/s,  below  which  they  are  assigned  a  value  of  0  m/s.  Li  the  ENVC  model,  wind 
speed  values  are  only  allowed  to  vary  to  a  minimum  of  1  m/s,  below  which  they 
are  assigned  a  value  of  1  m/s.  The  direction  of  change  observed  by  using  a  1  m/s 
hmit  to  lowest  wind  speed  values  used  in  the  /VENV  model  was  to  bring  model 
results  in  closer  agreement  to  ENVC  model  results.  This  limit  is  recommended 
for  future  use  in  the  AENV  model. 

8.  One  potential  limitation  of  the  AENV  model  is  a  requirement  to  calculate  surface 
roughness  length  (Zo)  as  a  monthly  average  using  hourly  data  where  the  wind 
speed  is  >6  m/s.  At  issue  is  a  situation  where  only  a  small  number  of  hourly 
values  exist  during  a  month  in  which  the  wind  speed  is  >6  m/s.  Greater 
uncertainty  exists  in  the  calculated  average  Zo  value  when  computed  from  only  a 
small  number  of  hourly  values.  The  following  approach  is  recommended  for 
future  use  in  the  AENV  model  to  deal  with  this  limitation  in  most  cases  - 

•    Seven  hourly  observations  is  recommended  as  the  minimum  number  of 
observations  in  which  to  use  in  a  month  for  calculating  Zo. 
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•  The  average  Zq  value  of  two  surrounding  months  can  be  used  for  any  month  in 
which  an  insufficient  number  of  observations  (i.e.,  <7)  exist. 

•  If  two  or  more  consecutive  months  have  <7  observations  for  a  surface  cover 
whose  characteristics  are  not  expected  to  undergo  substantial  change  from 
season  to  season  (e.g.,  coniferous  forest  canopy),  then  an  annual  average  Zo 
value  can  be  used  to  represent  surface  roughness  length  for  these  months. 

9.  Wetness  alters  the  characteristics  of  a  surface  to  act  a  sink  for  depositing  species 
(primarily  gases).  The  AENV  model  requires  an  hourly  measurement  to  represent 
the  amount  of  wetness  on  a  surface  in  determination  of  surface  resistance  (R^). 
The  model  is  set  up  to  use  field  measurement  data  from  surface  wetness  sensing 
equipment.  However,  past  experience  has  indicated  that  this  field  equipment  does 
not  provide  reliable  measurement  data  in  an  Alberta  setting  due  to  the  wide  range 
of  climatic  conditions  experienced.  Alternatively,  the  model  allows  a  relative 
humidity  (RH)  measurement  cutoff  (87%)  to  be  used  as  a  gross  approximation  to 
represent  a  surface  as  being  wet.  An  evaluation  of  alternative  RH  cutoff  values 
indicated  little  change  in  monthly  summed  deposition  flux  for  1 1  species.  As  a 
result,  it  is  recommended  that  the  current  RH  cutoff  value  be  maintained  in  the 
model  for  future  use. 

10.  A  final  issue  was  examined  related  to  performance  of  the  ENVC  model.  The 
ENVC  model  has  an  apparent  limitation  in  the  way  it  treats  settling  properties  of 
species  associated  with  large  particles  (i.e.,  NO3',  Na^,  K^,  Ca~^,  and  Mg~^).  The 
model  forces  the  treatment  of  settling  for  >24%  of  these  particles  (by  mass)  as 
being  similar  to  settling  properties  of  fine  particles.  This  underestimates 
neutralization  capacity  offered  by  depositing  particulates  containing  base  cations; 
and  overestimates  total  deposition  flux  to  a  surface  when  all  depositing  species 
are  taken  into  account.  This  hnndtation  was  further  investigated  to  show  how  large 
the  effect  was  for  the  ENVC  model.  Because  the  ENVC  model  tends  to  estimate 
higher  total  deposition  fluxes  compared  to  the  AENV  model,  better  agreement 
among  these  two  models  was  achieved  when  corrections  were  made  to  the  ENVC 
model  to  address  this  Umitation. 

Table  19  summarizes  recommended  changes  to  the  AENV  model  for  future  use.  Figure 
56  shows  the  net  effect  of  implementing  recommended  changes  to  the  AENV  model 
relative  to  ENVC  model  results  for  a  coniferous  forest  cover  for  the  conditions  tested. 
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Table  19 


Recommended  changes  to  AENV  acid  deposition  inference  model. 


Parameter 

Original  AENV 
model  condition 

Revised  condition  evaluated 

Revised  condition 
recommended  for  future 
use  in  AENV  model 

aerodynamic 
boundary 
resistance  (Ra)  [m] 

upper  limit  on  Ra 
is  undefined 

upper  limit  capped  at  1000  s/m 
for  all  surface  covers  except 
water  (2000  s/m  for  water) 

yes 

upper  limit  on  Ra 
is  undefined 

lower  limit  capped  at  5  s/m  for 
all  surface  covers 

yes 

Monin-Obukhov 
length  (L)  [m] 

undefined  limits 

upper  limit  of  -5  m  on 
negative  side  and  lower  limit 
of  +5  m  on  positive  side 

yes 

wind  speed 
minimum  (u)  [m/s] 

set  to  0  m/s  if 
<0.09  m/s 

set  to  1  m/s  if  <1  m/s 

yes 

surface  roughness 
length  (Zo)  [m] 

computed  as 
monthly  average 
using  hourly  data 
where  the  wind 
speed  >6  m/s 

use  seven  hourly  observations 
as  minimum  number  of 
observations  in  a  month  to 
calculate  Zq  and  surrounding 
monthly  average  values  for  <7 
observations 

yes 

surface  wetness 
condition 

uses  relative 
humidity  (RH) 
measurement 
cutoff  of  87%  as 
approximafion  to 
represent  surface 
as  being  wet 

uses  alternative  (higher/lower) 
relative  humidity  (RH) 
measurement  cutoff  as 
approximafion  to  represent 
surface  as  being  wet 

no 
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Modified  AENV  Model 


Original  AENV  Model 


ENVC  Model 


Monthly  summed  deposition  flux  for  11  species  using  modified  AENV 
model  compared  to  original  AENV  model  and  ENVC  model  for  a 
coniferous  forest  cover. 
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APPENDIX  A     ENVIRONMENT  CANADA  CALCULATION 

METHODS  FOR  GASES  AND  PARTICULATES 


(after  Brook  et  al.,  1999;  Zhang  et  al.,  2003,  2002a,  2002b,  2001a,  and  2001b) 


Parameters  (for  gases): 

deposition  velocity  (m/s) 

Ra. 

aerodynamic  resistance  (s/m) 

boundary-layer  resistance  (s/m) 

Rc: 

surface  resistance  (s/m) 

k: 

von  Karman  constant  (0.4) 

u< 

friction  velocity  (m/s) 

u: 

horizontal  wind  speed  at  reference  height  (m/s) 

L: 

Monin-Obukhov  length 

H: 

sensible  heat  flux 

z: 

reference  height  (10m) 

zq: 

suridce  rougnness  lengui  yiiij 

siaDiiiiy  lunciion 

Ki: 

Bulk  Richardson  number 

gravitational  acceleration  (9.81  m/s^) 

Ti- 

temperature  at  reference  height  (  K) 

Tip' 

potential  temperature  at  reference  height  (°K) 

T'  . 

surface  temperature  (°K) 

T  • 

temperature  difference  (°K) 

P- 
r. 

surface  pressure  (kPa) 

T 

^avg- 

temperature  average,  7^^,^  -  (-'i/.  +i  J/2 

SR: 

solar  radiation  (w/m^) 

RH: 

relative  humidity  (%) 

Prec: 

precipitation  (mm/hour) 

FCLD: 

fraction  of  cloud  covering  (%) 

cosO : 

cosine  of  solar  zenith  angle 

cos  a : 

cosine  of  sun/leaf  angle,  set  as  0.5 

Di: 

molecular  diffusivity 

Dv: 

water  diffusivity 

Vi: 

thermal  diffusivity 

MW: 

Molecular  Weight,  MWair  =  29 

fraction  of  stomatal  blocking  under  wet  conditions 
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Rst'.  stomatal  resistance 

R^:  mesophyll  resistance  (dependent  only  on  chemical  species) 

Rac'.  in-canopy  aerodynamic  resistance  (not  chemical  species-dependent) 

Rg]  ground  resistance  (scaled  from  SO2  and  03's  Rg) 

Rcut'.  cuticle  uptake  resistance  (scaled  from  SO2  and  03's  R^ut) 

Gs(PAR):  unstressed  leaf  stomatal  conductance 

PAR:  Photosynthetic  Active  Radiation 

f(T) :  conductance-reducing  effects  of  air  temperature  T  (°C) 

/(D) :  conductance-reducing  effects  of  water- vapour-pressure  deficit  D 

D:  water-vapour-pressure  deficit 

f(y/) :  conductance-reducing  effects  of  water  stress  y/ 

y/ :  water  stress 

^^^1  and  i//'^.2 :  leaf-water-potential  dependency 

E^(T):  saturation  water  vapour  pressure  (kPa)  at  air  temperature  T  (°K) 

E:  ambient  water  vapour  pressure  (kPa) 

/■,„^„, :  snow  cover  fraction 

«/  snow 

Sd:  snow  depth  (cm) 

Sdmax'.  maximum  snow  depth  (cm) 

Rcutwo'-  reference  values  of  wet  cuticle  resistance,  see  Table  Al  for  Raano  (O3) 
Rcutwo(^02):     50  s/m  or  100  s/m  for  rain  or  dew  conditions,  respectively 

Rcutdo-  reference  values  of  dry  cuticle  resistance,  see  Table  Al 

LAI:  Leaf  Area  Index 

LUC:  Land  Use  Category  (26  LUCs  in  gas  model) 

SC:  SeasonalCategory  (5  in  particulate  model) 


Deposition  Velocity: 


Aerodynamic  Resistance  (Ra): 

ENVC  adapts  the  Acid  Deposition  and  Oxidant  Model  (ADOM)  formulated  by  Pleim  et. 
al.  (1984)  as  cited  in  Zhang  et  al.  (2001b)  and  further  investigated  by  Padro  et  al.  (1990) 
as  cited  in  Zhang  et  al.  (2001b). 
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R^^  -   [0.74  ln(z  /  Zo)-i/^h]^  with  a  lower  limit  of  5  and  an  upper  limit  of  1000  (upper 

limit  of  2000  if  water  or  lake  as  land  use) 

Friction  velocity: 


Bulk  Richardson  number: 

Ri  =         ,  with  a  fixed  value  of  1x10"'^  if  SR>0  and  Ri>0 

Temperature  difference: 

-^2p~  '  with  an  upper  limit  of  -10"^^  for  negative  values  and  a  lower  limit  of  10 
for  positive  values 

Stability  function: 


ku 


Stable,  Ri>0 


Unstable,  Ri<0 


stable,  0  <  —  <  1 
L 


0.74x21n^ 


+, 1-9- 
V  L 


 \      unstable,  -1  <  —  <  0 

2  L 


Monin-Obukhov  length: 
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kHg 
positive  side 


,  with  an  upper  limit  of  -5  on  the  negative  side  and  a  lower  limit  of  +5  on  the 


Sensible  heat  flux: 


H 


0.74J  [ln(z/zO)J   [(1  +  4.7/?/)- 
k   9ARi 


0.74J  [ln(z/zO)J        \  +  53[kl\n{zl  z^)]\  \Ri\ 


Stable,  Ri>0 


Unstable,  Ri<0 


Sub  Layer  Resistance  (Rh): 


Thermal  diffusivity: 


Vi  = 


145.8  X 10"' r 


+110.4 


Molecular  diffusivity: 


0.00  ir 


Di  = 


1  1 

+ 


V     air  gas  / 


D^^^  =0.369MW^.^+6.29 


Total  Surface  Resistance  (Rr); 


1  l-Wst 


+       ,  with  a  lower  limit  of  1 0  s/m 
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Fraction  of  stomatal  blocking  under  wet  conditions: 


Wst 


0  SR<  200Wm  "' 

(SR  -  200)  /  800       200  <SR<  600Wm 
0.5  SR  >  600Wm"^ 


Stomatal  resistance: 


"  G^{PAR)f{T)f{D)f{y/)DJD^ 
Water  diffusivity: 


D  =  — 


1  1 

+ 


MW  MW 


,]/3 

air     •  ~  water 


Unstressed  leaf  stomatal  conductance: 

(PAR)  =  F,„„  /  /-„  )  +        /  r„  {PAR,,^, ) 


F,„„=2cos4l-e- 

^shade  ~  ~  ^sun 


-0.5  LAI /cor  0 


rJPAR)  =  r^,^^i\  +  bJPAR) 


PAR 


PAR, 


_\Rdiffe^~'''^'  ^ +0.01R,.^i\A-0ALAI)e-''''  LAK  2.5  or  SPKlOOwm'' 
"  ~  1^^,//^^"''"^'"'^  +0.07/?,,,(l.l  -O.ILA/)^-^^^^    LAI  >  2.5    and    SR  >  200wm-' 

\ Rj-^  cosa/ cos 0+ PAR LAI  <  2.5    or    SolarRad  <  200wm~^ 
coscir/cos6>  +  PAR^^,^,^,^    LAI  >  2.5    and    SolarRad  >  200wm"' 


Conductance-reducing  effects  of  air  temperature  T: 
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Conductance-reducing  effects  of  water- vapour-pressure  deficit  D: 

f{D)  -  1  -  h^^^D  ,  with  a  lower  limit  of  0. 1  and  a  upper  limit  of  1 

Water-vapour-pressure  deficit: 
Saturation  vapor  pressure: 

17.27(7-273.16) 

£*(r)  =  0.6108^  (in  units  of  mb  in  code) 

Ambient  water  vapour  pressure: 

Conductance-reducing  effects  of  water  stress: 

1  1  ¥^¥.^ 

where  n/  -  -0.72  -  0.00135/? ,  with  a  lower  limit  of  0.1  and  upper  limit  of 


Mesophyll  resistance: 

Values  of  Rm  for  all  dry-depositing  species  in  AURAMS  gas-phase  chemical  mechanism 
listed  in  Table  Al. 


In-canopy  aerodynamic  resistance: 
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where 


Raco^  reference  value  for  in-canopy  aerodynamic  resistance  (Table  Al) 
LAI:  Leaf  area  index  (Table  A2) 

For  some  LUCs,  a  range  of  Raco  values  is  given  to  reflect  the  change  of  canopy  structure 
at  different  times  of  the  growing  season.  Raco  values  for  any  day  of  the  year  based  on 
minimum  and  maximum  LAI  values  given  as: 

^.0  (0  =       (min)  +  X  [^.co  (max)  -  R^^.,  (mm)] 

LAI  (Max)  -  LA/ (mm) 


Ground  resistance: 


a(i)    ^  J3(i) 


R(i)    R(SO,)  R(0,) 


where 


/?,(o,)  = 


2000  s  m  '    LUC  1-3  and  Snow  Surface 
200  s  m''        LUC  4-19,  25  and  26 
500  s  m''  LUC  20-24 


For  snow  surface,  Rg(03)  adjusted  by  including  a  snow  cover  fraction(  /,,„,,^,): 


1  _  i-2/:.„„  ^2/, 


RAO,)    RALUC)  ft 


where 

sd 
sd. 


note  -  both  /^^^^  and  2  /^^,^^  have  a  lower  limit  of  0  and  upper  limit  of  1 


max 

R  =  2000  s/m 


RASO,) 


20  s  m~' 
70(2 -r) 

50 

100 

R^,.  {Table  1) 


LUC  1  and  3 
LUC  2,  T  as  Surface  Temperature  (°C) 
LUC  4-26,  Rain 
LUC  4-26,  Dew 
LUC  4-26,  Not  Rain  or  Dew,  T>-\°C 


mm(2R^,,R^y-^^-''''^)    LUC  4-26,  Not  Rain  or  Dew,  T<-\°C 
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with  a  lower  limit  of  100  s/m  and  upper  limit  of  500  s/m  (limit  only  applies  to 
LUC2) 


For  snow  surface,  Rg(S02)  are  adjusted  by  snow  cover  fraction(  ): 


1  \-2f  2f 

J  siuw    _|_     J  snow 


where 

sd 

fsnow  '  ^^ote  -  both        and  2  /^^^^^^,  have  a  lower  limit  of  0  and  upper 

limit  of  1 

^„„=70(2-r),Tas°C 


Canopy  cuticle  resistance: 

R 


R 


where 


vutwQ 


R 


uidO 


Rain  or  Dew 

Not  Rain 
T  <-\°C 


RcutwoCOs)  listed  in  Table  Al. 

Rcutwo(S02)  =  50  s/m  or  100  s/m  for  rain  or  dew  conditions,  respectively. 
Rcutdo  listed  in  Table  Al. 

with  a  lower  limit  of  100  s/m  and  20  s/m  for  SO2  dry  and  wet  conditions, 
respectively 

For  snow  surface,  Rcut  are  adjusted  by  including  a  snow  cover  fraction  ( ) 


R. 


where 


R^ 


{   2000  O3 

l70(2-r)  so. 
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Roughness  length: 

^   .  ^      LAI  {t)- LAI  {mm)      r  -, 
Zo  (0  =  Zo  (mm)  +  .    X     (max)  -  Zq  (mm)J 

LAI  (Max)  -  LA/(mm) 


Rain  or  Dew  conditions: 

Condition 


Rain  T  >273A5°K  and  PreoOAmm/hr 
Dew  T>213A5°Kand  U   >  M  min 


where: 

\.5Coedew 


u 


nun 


DQ 

where: 


Coedew  - 


0.3        FCLD  <  0.25 
0.2    0.25  <  FCLD  <  0.75 
0.1      0.75  <  FCLD  <1 


and 

DQ  =  0.0622(1  -  RH)E^  (T^ ) ,  with  a  lower  limit  of  0.0001 

Parameters  (for  particulates): 

(Variables  for  computing  Ra  were  not  included) 


Vd'. 

deposition  velocity  (m/s) 

Ra: 

aerodynamic  resistance  (s/m) 

Rs: 

surface  resistance  (s/m) 

u*: 

friction  velocity  (m/s) 

u: 

horizontal  wind  speed  at  reference  height  (m/s) 

g'- 

gravitational  acceleration  (9.81  m/s^) 

T2: 

temperature  at  reference  height  (°K) 

LAL 

Leaf  Area  Index 

LUC: 

Land  Use  Category  (15  LUCs  in  particle  model) 

SC: 

Season  Category  (5  SCs  in  particle  model,  listed  in  Table  A5) 

P : 

density  of  particle 

dp: 

particle  diameter 

C: 

correction  factor  for  small  particles 
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jX :  dynamic  viscosity  of  air 

V :  kinematic  viscosity  of  air 

C:  Cunningham  slip  correction  factor 

Eb:  collection  efficiency  from  Brownian  diffusion 

EiM^  collection  efficiency  from  impaction 

Ein:  collection  efficiency  from  interception 

Rl:  correction  factor  representing  fraction  of  particles  that  sticks  to  surface 

Sc:  Schmidt  number 

D\  Brownian  diffusivity 

St.  Stokes  number 

K\  Boltzmann  constant  (gm  cm"/s^-°K) 

P:  surface  pressure  (kPa) 

X :  mean  free  path  of  air  molecules 

Deposition  Velocity: 

Assumes  SO4'''  and  NH/  have  the  same  Vd;  and  NO3",  Ca"^,  Mg'^,  Na^,  and  have  the 
same  Vd  as  described  in  Cheng  et  al.  (2001). 

Mass  median  diameter  (MMD)  and  geometric  standard  deviation  (GSD)  for  S04~"  chosen 
as  0.35  jam  and  2.0  |am,  respectively  as  reported  in  Wesely  et  al.  (1985). 

MMD  and  GSD  for  Na"^  are  taken  as  5.12  |jm  and  2.64  |am,  respectively  as  reported  in 
Ruijgrok  et  al  (1997). 

Gravitational  settling  velocity: 


Cunningham  slip  correction  factor: 


C  =  l  +— (1.257  +  0.4^-°""'") 
d  „ 
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Aerodynamic  Resistance: 

Ra  is  computed  with  the  same  approach  as  the  gas  model,  but  the  parameters  are  adapted 
from  Table  A4. 

Surface  Resistance: 

R  =  !  ,  with  a  lower  limit  of  5 

where  £q  is  an  empirical  constant  chosen  as  3  for  all  LUCs 

Collection  efficiency  from  Brownian  diffusion: 


where  y  lies  between  V2  and  ^/^  with  larger  values  for  rougher  surfaces.  (Table  A4) 


Schmidt  number: 


Sc  =  v/D 


Collection  efficiency  from  impaction: 


where 

a:  constant  varying  with  LUC  (Table  A4) 
b:  constant,  chosen  as  2 


Stokes  number: 


smooth  surface,  LUC  8,  9,  12,  13,  and  14 


V 


vegetated  surface,  LUCs  other  than  above 


g 
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where  A  is  the  characteristic  radius  varying  with  LUC  and  season  (Table  A4) 


Collection  efficiency  from  interception: 


with  a  upper  Hmit  of  0.6 


where  A  is  the  characteristic  radius  varying  with  LUC  and  season  (Table  A4) 

Correction  factor  representing  fraction  of  particles  that  sticks  to  the  surface: 

=  expi-St^'^) ,  with  a  low  limit  of  0.5 


Mean  free  path  of  air  molecules: 


;L  =  6.54x10" 


ju  101.3 


1.818x10"'  P 


v293.15y 


Brownian  diffusivity: 

CKT, 


D 


3miu 


Dynamic  viscosity  of  air: 


145.8xl0"'r/^' 
+110.4 


Kinematic  viscosity  of  air: 


V  = 


Bin  division: 

For  I  =  1  to  40, 
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binsize  =  0.001'' I 


*  ^2-5 


Units  required  for  input  parameters  to  Visual  Basic  Program 


hourly  pollutant  concentration: 
hourly  cloud  fraction: 
hourly  pressure: 
hourly  snow  depth: 
hourly  precipitation: 
hourly  solar  radiation: 
hourly  relative  humidity: 
hourly  wind  speed: 
hourly  temperature: 


tenths  (converted  to  decimal  within  program) 

kPa  (converted  to  mb  within  program) 

cm 

mm 

W/m^ 

%  (converted  to  decimal  within  program) 
km/hr  (converted  to  m/s  within  program) 
°C  (converted  to  K  within  program) 
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Table  AILUC  and  all  related  parameters  (all  resistance  have  units  of  s/m;  na  =  not 
applicable;  f(u)  =  function  of  wind  speed). 


r\acO 

Rcut  dO 
03 

Rcu,wO 
03 

RculdO 

SO2 

RgU 
SO2 

rslmin 

(S/m) 

brs 

(Wm-2) 

Tn,n 

(°C) 

Tniaji 

(X) 

Topi 
(°C) 

b,po 
(kPa-') 

(Mpa) 

¥c2 

(Mpa) 

Zo 
(m) 

If 

1 

water 

0 

na 

na 

na 

20 

na 

na 

na 

na 

na 

na 

na 

na 

f(u) 

na 

2 

ice 

0 

na 

na 

na 

Eq,  (8a) 

na 

na 

na 

na 

na 

na 

na 

na 

0.01 

3 

inland  lal<e 

0 

20 

f(u) 

na 

4 

evergreen  needleleaf  trees 

100 

4000 

200 

2000 

200 

250 

44 

-5 

40 

15 

0.31 

-2 

-25 

0.9 

200 

5 

evergreen  broadleaf  trees 

250 

6000 

400 

2500 

100 

150 

40 

0 

45 

30 

0.27 

-1 

-50 

2,0 

400 

6 

deciduous  needleleaf  trees 

60-100 

4000 

200 

2000 

200 

250 

44 

-5 

40 

15 

0.31 

-2 

-25 

04-09 

200 

7 

deciduous  broadleaf  trees 

100-250 

6000 

400 

2500 

200 

150 

43 

0 

45 

27 

0.36 

-1.9 

-2.5 

0  4-1,0 

200 

8 

tropical  broadleaf  trees 

300 

6000 

400 

2500 

100 

150 

40 

0 

45 

30 

0.27 

-5,0 

2  5 

400 

9 

drought  deciduous  trees 

100 

8000 

400 

6000 

300 

250 

44 

0 

45 

25 

0,31 

-1 

-4,0 

0.6 

200 

10 

evergreen  broadleaf  shrubs 

60 

6000 

400 

2000 

200 

150 

40 

0 

45 

30 

0.27 

-2 

-4,0 

0.2 

50 

11 

deciduous  shrubs 

20-60 

5000 

300 

2000 

200 

150 

44 

-5 

40 

15 

0.27 

-2 

-4,0 

0.05- 
02 

50 

12 

thorn  shrubs 

40 

5000 

300 

2000 

200 

250 

44 

0 

45 

25 

0.27 

-2 

-35 

0.2 

50 

13 

short  grass  and  forbs 

20 

4000 

200 

1000 

200 

150 

50 

5 

40 

30 

0 

-1.5 

-2,5 

0,04 

5 

14 

long  grass 

10^0 

4000 

200 

1000 

200 

100 

20 

5 

45 

25 

0 

-1.5 

-2,5 

0,02- 

0,1 

20 

15 

crops 

10-40 

4000 

200 

1500 

200 

120 

40 

5 

45 

27 

0 

-1.5 

-2.5 

0,02- 
0,1 

10 

16 

rice 

10^0 

4000 

200 

1500 

50 

120 

40 

5 

45 

27 

0 

-1.5 

-25 

0,02- 
0,1 

10 

17 

sugar 

10^0 

4000 

200 

2000 

200 

120 

50 

5 

45 

25 

0 

-1.5 

-2,5 

0,02- 
0,1 

10 

18 

maize 

10-50 

5000 

300 

2000 

200 

250 

65 

5 

45 

25 

0 

-1,5 

-2.5 

0,02- 
01 

10 

19 

cotton 

10^0 

5000 

300 

2000 

200 

125 

65 

10 

45 

30 

0 

-1.5 

-2,5 

0,02- 
0,2 

10 

20 

irrigated  crops 

20 

4000 

200 

2000 

50 

150 

40 

5 

45 

25 

0 

-1.5 

-2,5 

0.05 

10 

21 

urban 

40 

6000 

400 

4000 

300 

200 

42 

0 

45 

22 

0.31 

-1.5 

-3 

1,0 

50 

22 

tundra 

0 

8000 

400 

2000 

300 

150 

25 

-5 

40 

20 

0.24 

0 

-1,5 

0,03 

2 

23 

swamp 

20 

5000 

300 

1500 

50 

150 

40 

0 

45 

20 

0.27 

-1.5 

-2,5 

0.1 

10 

24 

Desert 

0 

na 

na 

na 

700 

na 

na 

na 

na 

na 

na 

na 

na 

0,04 

2 

25 

mixed  wood  forests 

100 

4000 

200 

2500 

200 

150 

44 

-3 

42 

21 

0.34 

-2 

-2.5 

0,6-0,9 

200 

26 

Transitional  forest 

100 

4000 

200 

2500 

200 

150 

43 

0 

45 

25 

0.31 

-2 

-3 

0,6-0,9 

200 
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Table  A2  LAI(i)  values  dependent  upon  Land  Use  Category  (LUC)  ("i"  represents 

month  number,  "i=14"  and  "i=15"  represent  minimum  and  maximum  value, 
respectively). 


1>0. 

1 
1 

Z 

J 

J 

O 

7 

o 
o 

Q 

1  fl 

1  1 
1  1 

1 1 
1  Z 

1  2 
1  J 

1  A 
14 

1  c 

1 

water 

0 

0 

0 

0 

2 

ice 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

inland  lake 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

evergreen  needleleaf  trees 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

evergreen  broadleaf  trees 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

deciduous  needleleaf  trees 

0.1 

0.1 

0.5 

1.0 

2.0 

4.0 

5.0 

5.0 

4.0 

2.0 

1.0 

0.1 

0.1 

0.1 

5.0 

7 

deciduous  broadleaf  trees 

0.1 

0.1 

0.5 

1.0 

2.0 

4.0 

5.0 

5.0 

4.0 

2.0 

1.0 

0.1 

0.1 

0.1 

5.0 

8 

tropical  broadleaf  trees 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

9 

drought  deciduous  trees 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

10 

evergreen  broadleaf  stirubs 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

11 

deciduous  shrubs 

0.5 

0.5 

1.0 

1.0 

1.5 

2.0 

3.0 

3.0 

2.0 

1.5 

1.0 

0.5 

0.5 

0.5 

3.0 

12 

thorn  shrubs 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

13 

short  grass  and  forbs 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

14 

long  grass 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

1.0 

2.0 

2.0 

1.5 

1.0 

1.0 

0.5 

0.5 

2.0 

15 

0.1 

0.1 

0.1 

0.5 

1.0 

2.0 

3.0 

3.5 

4.0 

0.1 

0.1 

0.1 

0.1 

0.1 

4.0 

16 

rice 

0.1 

0.1 

0.1 

0.5 

1.0 

2.5 

4.0 

5.0 

6.0 

0.1 

0.1 

0.1 

0.1 

0.1 

6.0 

17 

sugar 

0.1 

0.1 

0.1 

0.5 

1.0 

3.0 

4.0 

4.5 

5.0 

0.1 

0.1 

0.1 

0.1 

0.1 

5.0 

18 

maize 

0.1 

0.1 

0.1 

0.5 

1.0 

2.0 

3.0 

3.5 

4.0 

0.1 

0.1 

0.1 

0.1 

0.1 

4.0 

19 

cotton 

0.1 

0.1 

0.1 

0.5 

1.0 

3.0 

4.0 

4.5 

5.0 

0.1 

0.1 

0.1 

0.1 

0.1 

5.0 

20 

irrigated  crops 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

21 

urban 

0.1 

0.1 

0.1 

0.1 

0.5 

1.0 

1.0 

1.0 

1.0 

1.0 

0.4 

0.1 

0.1 

0.1 

1.0 

22 

tundra 

1.0 

1.0 

0.5 

0.1 

0.1 

0.1 

0.1 

1.0 

2.0 

1.5 

1.5 

1.0 

1.0 

0,1 

2.0 

23 

swamp 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

24 

Desert 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

25 

mixed  wood  forests 

3.0 

3.0 

3.0 

4.0 

4.5 

5.0 

5.0 

5.0 

4.0 

3.0 

3.0 

3.0 

3.0 

3.0 

5.0 

26 

Transitional  forest 

3.0 

3.0 

3.0 

4.0 

4.5 

5.0 

5.0 

5.0 

4.0 

3.0 

3.0 

3.0 

3.0 

3.0 

5.0 
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Table  ASSpecies-specific  chemical  and  physical  parameters.^ 


No. 

Symbol 

Name 

Depn  Mmts 

MW 

H  (M  atm  '') 

H*  (M  atm-1) 

pe^'  (W) 

Rm 

a 

b 

1 

SO2 

Sulphur  dioxide 

Yes 

64 

(1.1-1.5)x10o 

2.65x105 

75.5  to  -7.6 

0 

1 

0 

2 

H2SO4 

Sulphuric  acid 

No 

98 

2.1  X  105 

>2.1  x105 

4.9  to  4.3 

0 

1 

1 

F 

NO 

Nitric  oxide 

Yes 

30 

(1.4-1.9)  x1 0-3 

H 

27.8  to  3.0 

Not  considered 

3 

NO2 

Nitrogen  dioxide 

Yes 

46 

(0.74.1)  X  10-2 

H 

28.4  to  8.2 

0 

0 

0.8 

4 

O3 

Ozone 

Yes 

48 

(0.9-1.3)  X  10-2 

H 

28.1  to  18.6 

0 

0 

1 

5 

H2O2 

Hydrogen  peroxide 

Yes 

34 

(0.7-1.4)  X  105 

H 

24.8  to  9.7 

0 

1 

1 

6 

HNO3 

Nitric  acid 

Yes 

63 

(0.1-2.6)  X  106 

3.2x1013 

14.1  to  8.9 

0 

10 

10 

7 

MONO 

Nitrous  acid 

Yes 

47 

(3.7-5.0)  X  101 

2.6x105 

17.5  to  14.8 

0 

2 

2 

8 

HNO4 

Pernitric  acid 

No 

79 

(0.1-1)  X  105 

>1  x107 

No  data 

0 

5 

5 

9 

NH3 

Ammonia 

Yes 

17 

(1.0-7.8)  X  101 

1.1  X  104 

Not  applicable 

0 

1 

0 

10 

PAN 

Peroxyacetylnitrate 

Yes 

121 

(2.8-5.0)  X  10" 

H 

30.2  to  -1.5 

0 

0 

0.6 

11 

PPN 

Peroxypropylnitrate 

No 

135 

2.9x10° 

H 

37.8  to  -2.3 

0 

0 

0.6 

12 

APAN 

Aromatic  acylnitrate 

No 

183 

No  data 

5 

46.9  to  11.2 

0 

0 

0.8 

13 

MPAN 

Peroxymethacrylic  nitric  anhydride 

No 

147 

1.7x10" 

H 

3.1 

0 

0 

0.3 

14 

HCHO 

Formaldehyde 

Yes 

30 

(0.3-1 .4)  X  10^ 

4.9x103 

3.0  to  -0.1 

0 

0.8 

0.2 

15 

MCHO 

Acetaldehyde 

No 

44 

(1.0-1.7)  x  101 

>15 

-1.0  to  -3.9 

100 

0 

0.05 

16 

PALD 

C3  carbonyls 

No 

58 

(2.4-3.7)  X  10" 

H 

-1.3  to -1.8 

100 

0 

0.05 

17 

C4A 

C4-G5  carbonyls 

No 

72 

(0.9-1 .8)  X  101 

H 

-1.3  to -1.8 

100 

0 

0.05 

18 

C7A 

06-08  carbonyls 

No 

128 

(0.4-11)  X  101 

H 

-1.5 

100 

0 

0.05 

19 

ACHO 

Aromatic  carbonyls 

No 

106 

(3.5-4.2)  X  101 

H 

-1.0  to  -2.5 

100 

0 

0.05 

20 

MVK 

Methyl-vinyl-ketone 

No 

70 

(2.1-4.4)  X  101 

H 

0.2 

0 

0 

0.05 

21 

MACR 

Methacrolein 

No 

70 

(4.3-6.5)  X  10° 

H 

-1 .2 

100 

0 

0.05 

22 

MGLY 

Methylgloxal 

No 

72 

(0.4-3.2)  xlO^ 

H 

-0.7 

0 

0.01 

0 

23 

MOH 

Methyl  alcohol 

No 

32 

(1.4-2.3)  X  102 

>H 

3.0 

0 

0.6 

0.1 

24 

ETOH 

Ethyl  alcohol 

No 

46 

(1.2-2.3)  X  102 

>H 

-1.3  to -2.9 

0 

0.6 

0 

25 

POH 

03  alcohol 

No 

60 

(0.9-1 .7)  X  102 

>H 

-0.3 

0 

04 

0 

26 

ORES 

Oresol 

No 

104 

8.3x102 

H 

-2.5 

0 

0.01 

0 

27 

FORM 

Formic  acid 

Yes 

46 

(0.9-8.9)  X 103 

9.8x105 

1.9  to -64 

0 

2 

0 

28 

ACAC 

Acetic  acid 

Yes 

60 

(0.8-9.3)  x  103 

9.6x105 

-3.1  to  -9.6 

0 

1.5 

0 

29 

ROOH 

Organic  peroxides 

Yes 

48 

(0.1-3.1)  x  102 

H 

4.2  to  3.6 

0 

0.1 

0.8 

30 

ONIT 

Organic  nitrates 

No 

77 

2.0x10° 

H 

10.5  to -5.0 

100 

0 

0.5 

31 

INIT 

Isoprene  nitrate 

No 

147 

2.0x10° 

H 

No  data 

100 

0 

0.5 

3  Values  for  Henry's  Law  constant  (H)  were  obtained  from  Howard  and  Meylan  (1997)  and  Sander  (1999). 
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Table  A4Parameters  related  to  LUC  and  SC  in  particle  model. 


NO. 

1  and  ii'^p  p^itpnnrip^  Dp^rrintinn 
(LUC) 

SC  1 

SC2 

Zo  (m) 
SC3 

SC4 

SC5 

SC  1 

SC2 

A  (rnm) 
SC3 

SC4 

SC5 

1 

Evergreen}needlele3f  trees 

0.8 

0.9 

0.9 

0.9 

0.8 

2 

2 

2 

2 

2 

1 

0.56 

2 

Evergreen  brosdiesf  trees 

2.65 

2.65 

2.65 

2.65 

2.65 

5 

5 

5 

5 

5 

0.6 

0.58 

3 

Deciduous  needleleaf  trees 

0.85 

0.85 

0.8 

0.55 

0.6 

2 

2 

5 

5 

2 

1.1 

0.56 

4 

Deciduous  broadleaf  trees 

1.05 

1.05 

0.95 

0.55 

0.75 

5 

5 

10 

10 

5 

0.8 

0.56 

5 

Mixed  broadleleaf  and  needleaf  trees 

1.15 

1.15 

1.15 

1.15 

1.15 

5 

5 

5 

5 

5 

0.8 

0.56 

6 

Grass 

0.1 

0.1 

0.05 

0.02 

0.05 

2 

2 

5 

5 

2 

1.2 

0.54 

7 

Crops,  mixed  farrninQ 

0.1 

0.1 

0.02 

0.02 

0.05 

2 

2 

5 

5 

2 

1.2 

0.54 

8 

Desert 

0.04 

0.04 

0.04 

0.04 

0.04 

na 

na 

na 

na 

na 

50 

0.54 

9 

Tundra 

0.03 

0.03 

0.03 

0.03 

0.03 

na 

na 

na 

na 

na 

50 

0.54 

10 

Shrubs  and  interrupted  wood-lands 

0.1 

0.1 

0.1 

0.1 

0.1 

10 

10 

10 

10 

10 

1.3 

0.54 

11 

Wet  land  with  plants 

0.03 

0.03 

0.02 

0.02 

0.03 

10 

10 

10 

10 

10 

2 

0.54 

12 

Ice  cap  and  glacier 

0.01 

0.01 

0.01 

0.01 

0.01 

na 

na 

na 

na 

na 

50 

0.54 

13 

Inland  water 

f(u) 

f(u) 

f(u) 

f(u) 

f(u) 

na 

na 

na 

na 

na 

100 

0.5 

14 

Ocean 

f{u) 

f(u) 

f(u) 

f(u) 

f{u) 

na 

na 

na 

na 

na 

100 

0.5 

15 

Urban 

1 

1 

1 

1 

1 

10 

10 

10 

10 

10 

1.5 

0.56 

Table  A5        Latitude  values  applied  for  the  determination  of  Seasonal  categories  (SC) 
used  in  particle  model. 


Refinement  Study  of  Dry  Deposition  Inference 
Methods  Used  in  Alberta 


102 


APPENDIX  B     ALBERTA  ENVIRONMENT  CALCULATION 

METHODS  FOR  GASES  AND  PARTICULATES 

(after  Cheng  et  al.,  2001  and  WBK,  2008) 


Parameters: 

B 

defined  equation  described  herein 

C 

concentration  (|ig/m"^) 

D 

diffusion  coefficient  of  the  substance  of  interest  (cm"/s) 

F 

dry  deposition  flux  (|ig/m^/s) 

g 

gravitational  acceleration  (9.81  m/s  ) 

H 

sensible  heat  flux  (w/m^) 

k 

von  Karman  constant  (0.4) 

L 

Monin-Obukhov  length  scale 

PAI 

Potential  Acid  Input  (kg  H^/ha/yr) 

Pr 

Prandtl  number  for  air  (0.72) 

Ra 

aerodynamic  resistance  (s/m) 

Rb 

boundary-layer  resistance  (s/m) 

Rc 

surface  resistance  (s/m) 

Ri 

Bulk  Richardson  number 

RH 

relative  humidity 

Sc 

Schmidt  number 

sw 

soil  wetness 

Td 

temperature  difference  between  10  and  2  m  (Tio  -  T2) 

T2 

temperature  at  2  m  (Kelvin) 

u 

wind  speed  (m/s) 

friction  velocity  (m/s) 

Vd 

deposition  velocity  (m/s) 

[X] 

concentration  of  X  chemical  species  deposited  (kg/ha/yr) 

z 

reference  height  (10  m) 

zo 

surface  roughness  length  (m) 

dynamic  viscosity  of  air  (18.0  x  10"^  N-s/m^  at  1  atm  and  25 °C) 

p 

density  of  air  (1.18  kg/m^  at  1  atm  and  25 °C) 

standard  deviation  of  wind  direction  (rad) 

integrated  stability  correction  term 
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Equations: 


PAI 


dry 


^  [SO2]  ,  [^^2]  ^  [hno,]  ^  [hno,]  ^  ^\sol]  ^  \no^^  [nh:] 

64         46  47  63  96  62  18 


64         46  47  63  96 

39        23  40  24 


deposition 
concentrations 
units  in  kg/ha/yr 


F  =  VdC 


Vd  = 


(R^+R,+R, 


Summary  of  Species  Specific  Deposition  Velocity  Formulae: 

1 


V 


d{S02. 


d  (NO. 


V 


d  [HNO^ ) 


d  [HNCh  \ 


+R 


h{S(h] 


.    Rc  is  treated  as  being  negligible  for  nitric  and 
nitrous  acid. 


d[SOf,NH;]  Dip 


dm,Ca^^,Ms^\K^,Nan  +  (0.5  X /?^(^^,.  j ) 


Aerodynamic  Resistance  (Ra); 


Ra  =  — -\\r\  —  - —  I }> (Ra  is  infinite  and  Vd  =  0  when  u  and  Td  are  zero) 

ku    \  Zn 


Friction  yelocity: 


u  =  -^^^  (this  relationship  is  used  as  an  initial  estimate  of  u"^  to  calculate  zo,  a  more 
precise  value  of  u  is  calculated  after  zq  is  obtained  -  refer  to  below) 


Refinement  Study  of  Dry  Deposition  Inference 
Methods  Used  in  Alberta 


104 


Zo  =  z  ^  ^ "  ^  (calculated  as  a  monthly  average  using  data  where  the  wind  speed  is  >6 
m/s) 


A  more  precise  value  of  u  calculated  after  Zo  is  obtained  (based  on  atmospheric 
conditions): 


u*  =  ku< 

In 

z 

(l  +  4.7/?/)  \  Stable  conditionsR-,  >  0 


ku 


In 


z 


9ARi 


i  +  (iab) 


ku 


u  - 


In 


r  \ 
z 


Unstable  conditionsRi  <  0 


Neutral  conditions  Rj  =  0 


Bulk  Richardson  number 

Ri  = 


T.u 


Stability  function  (based  on  atmospheric  conditions): 

5z 


i|/  = 


Stable  conditions 


\|/=  21n< 
i|/  =  0 


1  +  Jl 


\5z 
L 


Unstable  conditions 


Neutral  conditions 


Monin-Obukhov  length; 

kHg 
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Sensible  heat  flux  (based  on  atmospheric  conditions): 


H 


0.74 


In 


z 


(l  +  4.1  Ri] 


Neutral  and  Stable  conditions 


H 


0.74 


In 


z 


1  _   ^-^^^    [    Unstable  conditions 
(1  +  5.3B] 


Where  B  =  9.4 


In 


Boundary  Layer  Resistance  (Rh): 
Rb  (gases): 

R. 


_2_ 

ku. 


7]  1 

—  X  — 


pD  Pr 


Ri 


^  for  SO2  and  HNO3 


Rh  = 


6.18 


forN02 


Rh  = 


6.09 


for  HNO, 


Rb  (particulates): 

Rb  values  for  particulate  sulphate  are  obtained  from  scientific  literature  for  daytime  and 
nighttime  as  a  function  of  surface  type  and  weighted  according  to  average  day  length  for 
each  month  at  a  mid- Alberta  latitude  location  (54°N  latitude)  after  Cheng  and  Angle 
(1993)  as  cited  in  Cheng  et  al.  (2001). 
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Table  Bl        Boundary-Layer  Resistance  (s/cm)  for  Fine  Particulate  (S04~",  NH4'^),  Day 
Length  Weighted  Averages  at  54°N  Latitude  for  the  Middle  of  Each 
Month. 


Surface  Type 


Winter 
(Dec,  Jan,  Feb) 


Spring 
(Mar,  Apr,  May) 


Summer 
(Jun,  July,  Aug) 


Autumn 
(Sep,  Oct,  Nov) 


Dry 

Wet 

Dry 

Wei 

Dry 

Wcl 

Dry 

Wcl 

Deciduous  Forest 

16.9 

0 

5.4 

0 

1.3 

0 

3.2 

0 

Coniferous  Forest 

2.5 

0 

2.7 

0 

1.9 

0 

2.3 

0 

Wetland/Swamp* 

20.4 

0 

3.8 

0 

2.6 

0 

3.2 

0 

Grassland* 

20.4 

0 

5.6 

0 

3.9 

0 

4.7 

0 

Cropland* 

20.4 

0 

9.0' 

0 

3.9 

0 

7.9^^ 

0 

Urban^ 

33.9 

0 

10.9 

0 

2.6 

0 

6.3 

0 

Open  Water 

0 

0 

0 

0 

0 

0 

0 

0 

Snow/Ice 

20.4 

0 

*  in  winter,  wetland,  grassland,  and  cropland  treated  as  a  snow  surface.'  bare  soil  and  active  growth. 

*  bare  soil  and  senescent  growth.^  consists  of  a  mixture  of  deciduous  forest  and  buildings. 


Table  B2        Boundary-Layer  Resistance  (s/cm)  for  Coarse  Particulate  (NO3 ,  Na,  K, 
Ca,  Mg)  Day  Length  Weighted  Averages  at  54°N  Latitude  for  the  Middle 
of  Each  Month. 


Surface  Type 


Winter 
(Dec,  Jan,  Feb) 


Spring 
(Mar,  Apr,  May) 


Summer 
(Jun,  July,  Aug) 


Autumn 
(Sep,  Oct.  Nov) 


Dry 

Wet 

Dry 

Wet 

Dry 

Wet 

Dry 

Wet 

Deciduous  Forest 

8.45 

0 

2.7 

0 

0.65 

0 

1.6 

0 

Coniferous  Forest 

1.25 

0 

1.35 

0 

0.95 

0 

1.15 

0 

Wetland/Swamp* 

10.2 

0 

1.9 

0 

1.3 

0 

1.6 

0 

Grassland* 

10.2 

0 

2.8 

0 

1.95 

0 

2.35 

0 

Cropland* 

10.2 

0 

4.5' 

0 

1.95 

0 

3.95^ 

0 

Urban^ 

16.95 

0 

5.45 

0 

1.3 

0 

3.15 

0 

Open  Water 

0 

0 

0 

0 

0 

0 

0 

0 

Snow/Ice 

10.2 

0 

*  in  winter,  wedand,  grassland,  and  cropland  treated  as  a  snow  surface.'  bare  soil  and  active  growth. 
^  bare  soil  and  senescent  growth.^  consists  of  a  mixture  of  deciduous  forest  and  buildings. 


Day  length  Weighted  Seasonal  Average  Rb 


day  length 
lAhrs 


1+1 


day  length 
lAhrs 


nighl 


Day  length  =  0. 133{cos  '  (-  tan(55°  )x  tm{Solar  Declination)} 


c.  ,  ^...1  .  \  360x{2S4  + Julian  Day)^ 

Solar  Decimation  =  23.45^  sinl  ^  — 

365  J 
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Surface  Resistance  (Rr): 


Bulk  surface  resistance  values  are  used  from  literature  as  a  function  of  surface  type, 
surface  wetness,  and  incident  radiation.  Day  length  weighted  average  Rc  values  for  SO2 
and  NO2  are  used  from  Voldner  et  al.  (1986),  Arrit  et  al.  (1987)  and  Walcek  et  al.  (1986) 
as  cited  in  Cheng  et  al.  (2001): 

Table  B3        Day  Length  Weighted  Averages  Bulk  Surface  Resistance  (s/cm)  for 

 Sulphur  Dioxide  (SO2):  

Winter  Spring  Summer  Autumn 

Surface  Type  (Dec,  Jan,  Feb)         (Mar,  Apr,  May)         (Jun,  July,  Aug)         (Sep,  Oct,  Nov) 


Dry 

Wet 

Dry 

Wet 

Dry 

Wet 

Dry 

Wet 

Deciduous  Forest 

10 

10 

4.7 

0 

3.5 

0 

7.9 

0.4 

Coniferous  Forest 

5 

5 

4.1 

0 

3.5 

0 

4.9 

0.2 

Wetland/Swamp* 

7 

1 

0.5 

0 

0.7 

0 

1 

0.1 

Grassland* 

7 

1 

0 

1.3 

0 

2 

0.1 

Cropland* 

7 

1 

0+ 

0 

2 

0 

2* 

0.1 

Urban^ 

10 

2 

10 

0 

10 

0 

10 

0.1 

Open  Water 

0 

0 

0 

0 

0 

0 

0 

0 

Snow/Ice 

7 

1 

*  in  winter,  wetland,  grassland,  and  cropland  treated  as  a  snow  surface.^  bare  soil  and  active  growth. 
"  bare  soil  and  senescent  growth.^  consists  of  a  mixture  of  deciduous  forest  and  buildings. 


Table  B4        Day  Length  Weighted  Averages  Bulk  Surface  Resistance  (s/cm)  for 

 Nitrogen  Dioxide  (NO2):  

Winter  Spring  Summer  Autumn 

Surface  Type  (Dec,  Jan,  Feb)         (Mar,  Apr,  May)         (Jun,  July,  Aug)         (Sep,  Oct,  Nov) 


Dry 

Wet 

Dry 

Wet 

Dry 

Wet 

Dry 

Wet 

Deciduous  Forest 

20.0 

70.0 

3.3 

70.0 

2.2 

70.0 

4.7 

70.0 

Coniferous  Forest 

10.0 

70.0 

2.7 

70.0 

2.2 

70.0 

3.3 

70.0 

Wetland/Swamp* 

50.0 

70.0 

12.1 

70.0 

11.5 

70.0 

12.9 

70.0 

Grassland* 

50.0 

70.0 

3.3 

70.0 

3.3 

70.0 

6.6 

70.0 

Cropland* 

50.0 

70.0 

3.3' 

70.0 

4.6 

70.0 

7.9* 

70.0 

Urban^ 

10.0 

70.0 

10.0 

70.0 

10.0 

70.0 

10.0 

70.0 

Open  Water 

70.0 

70.0 

70.0 

70.0 

70.0 

70.0 

70.0 

70.0 

Snow/Ice 

50.0 

70.0 

*  in  winter,  wetland,  grassland,  and  cropland  treated  as  a  snow  surface.^  bare  soil  and  active  growth. 
^  bare  soil  and  senescent  growth.^  consists  of  a  mixture  of  deciduous  forest  and  buildings. 


Rc  (HNO3):  10  s/m  (for  all  seasons  and  all  surfaces) 

Rc  (HNO2):  10  s/m  (for  all  seasons  and  all  surfaces) 

Rc  (NH3):  28  s/m  (dry  conditions) 

9  s/m  (wet  conditions) 

201  s/m  (when  T2<0°C) 
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Rc  (particulates):        0  s/m 


Rc  is  calculated  based  on  surface  wetness  criteria,  such  that  it  either  represents  a  "total 
dry  condition,"  "total  wet  condition,"  or  "weighted  wet  condition"  using  the  following 
flowchart,  and  relative  humidity  (RH)  and  surface  wetness  (SW)  criteria: 


1 

RH  available 

Figure  Bl Flowchart  for  calculating  Rc  (Adapted  from  Bates  (1996) 

Default  =  Rc  value  for  dry  conditions 
Wet  Rc  =  Rc  value  for  wet  conditions 
Weighted  Wet  Rc  =  Time  weighted  wet  Rc 


Time  weighted  wet  Rc 


SW 
100 


X  wet  Rc 


1  \xdry  Rc 

100  ' 


Calculation  of  Vd  in  the  absence  of  meteorological  data: 

Missing  hourly  meteorological  data  are  treated  in  the  following  manner: 
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•  1  hour  of  meteorological  data  missing     the  average  resistance  of  the  hours  before 
and  after  are  used  to  represent  the  missing  hour 

•  consecutive  hours  of  meteorological  data  missing     each  hour's  calculated  median 
resistance  for  the  month  is  used  to  represent  the  missing  hours 

Units  required  for  input  parameters  to  Visual  Basic  Program 

hourly  pollutant  concentration: |ig/m^ 

hourly  wind  speed:  km/hr  (converted  to  m/s  within  program) 

hourly  wind  direction  std.  dev.:        degrees  (converted  to  radians  within  program) 
hourly  temperature:  °C  (converted  to  K  in  program) 

hourly  delta  temperature:  °C 

hourly  relative  humidity:  %  (converted  to  decimal  within  program) 

hourly  surface  wetness:  % 
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APPENDIX  C     MONTHLY  AVERAGE  SPECIES 
CONCENTRATIONS 


Table  CI  Monthly  average  chemical  concentrations. 


gaseous  chemicals 

(pg/m  ) 

particulate  chemicals 

(pg'm  ) 

S02 

N02 

HN02 

HN03 

S04 

NH4 

N03 

Na 

K 

Ca 

M  s: 

Month  1 

2.9 

19.1 

0.76 

1.46 

1.03 

0.65 

0.30 

0.14 

0.09 

0.19 

0,04 

Month  2 

3.7 

18.6 

0.14 

0.62 

1.01 

0.46 

0.15 

0.10 

0.01 

0.26 

0.04 

Month  3 

3.3 

10.9 

0.07 

0.61 

2.35 

0.69 

0.17 

0.12 

0.01 

0.35 

0.03 

Month  4 

4.3 

9.4 

0.11 

0.74 

0.95 

0.49 

0.14 

0.07 

0.03 

0.81 

0.04 

Month  5 

3.2 

6.1 

0.13 

0.41 

0.67 

0.31 

0.05 

0.04 

0.02 

0.38 

0.05 

Month  6 

2.3 

3.7 

0.01 

0.19 

0.53 

0.25 

0.04 

0.02 

0.03 

0.47 

0.04 

Month  7 

3.9 

5.0 

0.02 

0.25 

0.50 

0.23 

0.08 

0.05 

0.06 

0.79 

0.07 

Month  8 

4.2 

7.0 

0.06 

0.31 

0.70 

0.34 

0.07 

0.05 

0.04 

0.61 

0.07 

Month  9 

1.7 

7.1 

0.11 

0.21 

0.59 

0.29 

0.05 

0.06 

0.04 

0.23 

0.04 

Month  10 

2.6 

10.7 

0.27 

1.24 

0.97 

0.43 

0.07 

0.05 

0.04 

0.88 

0.06 

Month  1 1 

1.0 

14.4 

0.14 

0.19 

0.79 

0.66 

0.48 

0.10 

0.07 

1.92 

0.06 

Month  12 

4.9 

21.7 

0.80 

0.34 

0.68 

0.34 

0.38 

0.10 

0.09 

0.55 

0.05 

25.0 


0.0  A — ^ — ^ — ^ — I         T    f    7    r  I — T — ^ — ^ — I 

0       1       2       3      4       5       6      7       8       9      10     11      12     13  14 


month 

I     S02      N02      HNQ3      S04 1 

Figure  CI       Monthly  averaged  chemical  concentrations  for  the  largest  N  and  S 
contributors. 
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2.50 


2.00 


0       1       2       3       4       5       6       7       8       9      10     11      12     13  14 

month 

I  —  K      Na ->-Ca -*-Mg| 
Figure  C3Monthly  averaged  chemical  concentrations  for  base  cations. 
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APPENDIX  D 


TABLES  FOR  MONTHLY  DEPOSITION 
VELOCITIES  OF  COMPONENT  CHEMICALS 


Table  Dl  Monthly  average  deposition  velocities  for  SO2  -  coniferous  forest  (cm/s). 


month 

AENV 

ENVC 

difference 

Month  1 

0.09 

0.14 

-0.06 

Month  2 

0.10 

0.12 

-0.03 

Month  3 

0.13 

0.23 

-0.11 

Month  4 

0.17 

0.34 

-0.16 

Month  5 

0.25 

0.37 

-0.12 

Month  6 

0.49 

0.36 

0.13 

Month  7 

0.40 

0.30 

0.11 

Month  8 

0.35 

0.29 

0.06 

Month  9 

0.38 

0.38 

0.01 

Month  10 

0.36 

0.35 

0.01 

Month  1 1 

0.12 

0.19 

-0.07 

Month  12 

0.13 

0.18 

-0.05 

/erage  deposition  velocities  for  NO2  - 

-  coniferous 

month 

AENV 

ENVC 

difference 

Month  1 

0.05 

0.08 

-0.04 

Month  2 

0.05 

0.07 

-0.02 

Month  3 

0.09 

0.15 

-0.06 

Month  4 

0.18 

0.22 

-0.04 

Month  5 

0.20 

0.24 

-0.04 

Month  6 

0.25 

0.24 

0.01 

Month  7 

0.24 

0.20 

0.04 

Month  8 

0.24 

0.19 

0.05 

Month  9 

0.15 

0.23 

-0.08 

Month  10 

0.11 

0.20 

-0.09 

Month  1 1 

0.08 

0.11 

-0.03 

Month  12 

0.07 

0.10 

-0.03 
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Table  D3        Monthly  average  deposition  velocities  for  HNO2  -  coniferous  forest 
(cm/s). 


month 

AENV 

ENVC 

difference 

Month  1 

0.94 

0.32 

0.62 

Month  2 

1.15 

0.27 

0.87 

Month  3 

1.68 

0.50 

1.18 

Month  4 

1.78 

0.68 

1.10 

Month  5 

2.17 

0.73 

1.45 

Month  6 

2.04 

0.69 

1.35 

Month  7 

1.83 

0.58 

1.24 

Month  8 

1.73 

0.57 

1.17 

Month  9 

2.14 

0.76 

1.38 

Month  10 

1.51 

0.69 

0.82 

Month  1 1 

1.29 

0.42 

0.87 

Month  12 

1.07 

0.39 

0.67 

Table  D4       Monthly  average  deposition  velocities  for  HNO3  -  coniferous  forest 
(cm/s). 


month 

AENV 

ENVC 

difference 

Month  1 

0.88 

0.86 

0.02 

Month  2 

1.07 

0.82 

0.25 

Month  3 

1.57 

1.25 

0.32 

Month  4 

1.68 

1.59 

0.09 

Month  5 

2.05 

1.58 

0.47 

Month  6 

1.92 

1.51 

0.41 

Month  7 

1.71 

1.33 

0.38 

Month  8 

1.63 

1.30 

0.32 

Month  9 

2.01 

1.59 

0.42 

Month  10 

1.42 

1.32 

0.09 

Month  1 1 

1.21 

1.11 

0.10 

Month  12 

1.00 

1.05 

-0.04 
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Table  D5       Monthly  average  deposition  velocities  for  S04^7NH4^  -  coniferous  forest 
(cm/s). 


month 

AENV 

ENVC 

difference 

Month  1 

0.17 

0.09 

0.08 

Month  2 

0.19 

0.09 

0.10 

Month  3 

0.20 

0.12 

0.08 

Month  4 

0.25 

0.15 

0.10 

Month  5 

0.23 

0.14 

0.09 

Month  6 

0.35 

0.13 

0.22 

Month  7 

0.34 

0.10 

0.24 

Month  8 

0.33 

0.10 

0.23 

Month  9 

0.29 

0.11 

0.18 

Month  10 

0.25 

0.10 

0.16 

Month  1 1 

0.26 

0.10 

0.16 

Month  12 

0.24 

0.09 

0.15 

Table  D6       Monthly  average  deposition  velocities  for  N03VNa''/KVCa^^/Mg^^  - 
coniferous  forest  (cm/s) 


month 

AENV 

ENVC 

difference 

Month  1 

0.30 

0.30 

0.00 

Month  2 

•  0.35 

0.26 

0.09 

Month  3 

0.38 

0.38 

-0.01 

Month  4 

0.46 

0.51 

-0.05 

Month  5 

0.44 

0.48 

-0.04 

Month  6 

0.64 

0.41 

0.23 

Month  7 

0.63 

0.29 

0.33 

Month  8 

0.60 

0.32 

0.28 

Month  9 

0.54 

0.36 

0.17 

Month  10 

0.45 

0.30 

0.16 

Month  1 1 

0.46 

0.31 

0.15 

Month  12 

0.42 

0.28 

0.15 
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Table  DVSmall  and  large  particle  deposition  velocities  for  AENV  LUCs. 


deciduous  forest 

coniferous  forest 

wetland/swamp 

grassland 

cropland 

urban 

water 

snow/ice 

small 

small 

small 

small 

small 

small 

small 

small 

Month  1 

0.03 

0.17 

0.03 

0.03 

0.03 

0.02 

2.85 

0.03 

Month  2 

0.03 

0.19 

0.03 

0.03 

0.03 

0.02 

3.49 

0.03 

Month  3 

0.10 

0.20 

0.14 

0.10 

0.06 

0.05 

5.97 

Month  4 

0.13 

0.25 

0.18 

0.13 

0.08 

0.07 

4.84 

Month  5 

0.12 

0.23 

0.17 

0.12 

0.08 

0.06 

8.21 

Month  6 

0.49 

0.35 

0.26 

0.18 

0.18 

0.26 

6.70 

Month  7 

0.48 

0.34 

0.26 

0.18 

0.18 

0.26 

5.77 

Month  8 

0.46 

0.33 

0.25 

0.17 

0.17 

0.25 

5.28 

Month  9 

0.21 

0.29 

0.21 

0.15 

0.09 

0.11 

8.05 

Month  10 

0.19 

0.25 

0.19 

0.13 

0.08 

0.10 

4.67 

Month  1 1 

0.19 

0.26 

0.19 

0.14 

0.09 

0.11 

3.22 

Month  12 

0.04 

0.24 

0.04 

0.04 

0.04 

0.02 

2.33 

0.04 

large 

large 

large 

large 

large 

large 

large 

large 

Month  1 

0.06 

0.30 

0.05 

0.05 

0.05 

0.03 

2.85 

0.05 

Month  2 

0.06 

0.35 

0.05 

0.05 

0.05 

0.03 

3.49 

0.05 

Month  3 

0.20 

0.38 

0.28 

0.19 

0.12 

0.10 

5.97 

Month  4 

0.25 

0.46 

0.34 

0.24 

0.15 

0.13 

4.84 

Month  5 

0.23 

0.44 

0.32 

0.23 

0.14 

0.12 

8.21 

Month  6 

0.88 

0.64 

0.49 

0.34 

0.34 

0.49 

6.70 

Month  7 

0.86 

0.63 

0.48 

0.34 

0.34 

0.48 

5.77 

Month  8 

0.82 

0.60 

0.46 

0.32 

0.32 

0.46 

5.28 

Month  9 

0.40 

0.54 

0.40 

0.28 

0.17 

0.22 

8.05 

Month  10 

0.34 

0.45 

0.34 

0.25 

0.16 

0.19 

4.67 

Month  1 1 

0.35 

0.46 

0.35 

0.25 

0.16 

0.20 

3.22 

Month  12 

0.08 

0.42 

0.07 

0.07 

0.07 

0.04 

2.33 

0.07 
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APPENDIX  E 


MONTHLY  AVERAGED  DIURNAL  DEPOSITION 
VELOCITIES  FOR  SO2,  NO2,  AND  HNO3 


SO2 


January 


Iw'  0.8 

1  0.6 

>  0.4 

CM 

o 

(f>  0.2 


0  4—1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — ' — I — I — r— 


♦  ENVC 
■  AENV 


3    5    7    9   11  13  15  17  19  21  23 
Hour  of  day  (Jan.) 

February 


^  0.8 
I  0.6 


0.4  -- 


0.2 
0 


♦  ENVC 
■  AENV 


3    5    7    9  11  13  15  17  19  21  23 
Hour  of  day  (Feb.) 

March 


0.8  - 

To 

0.6  -- 

> 

CM 

0.4  -- 

o 

<n 

0.2^^ 

0  + 

♦  ENVC 
■  AENV 


3    5    7    9  11  13  15  17  19  21  23 
Hour  of  day  (Mar.) 
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April 


o  0.6 


0.4 
0.2 

0  4 


■■■■■II 


♦  ENVC 
■  AENV 


1    3    5    7    9  11  13  15  17  19  21  23 
Hour  of  day  (Apr.) 


May 


♦  ENVC 
■  AENV 


1    3    5    7    9   11  13  15  17  19  21  23 
Hour  of  day  (May) 


June 


^  0.8 
o  0.6 


0.2 

0  -\ 


0.4  


 1- 


■  ■ 


♦  ENVC 
■  AENV 


1    3    5    7    9  11  13  15  17  19  21  23 
Hour  of  day  (Jun.) 
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July 


0.8  + 
1  0.6 


0.4 
0.2  t 
0 


♦  ENVC 
■  AENV 


1    3    5    7    9  11  13  15  17  19  21  23 
Hour  of  day  (Jul.) 


August 


♦  ENVC 
■  AENV 


3    5    7    9   11  13  15  17  19  21  23 
Hour  of  day  (Aug.) 


September 


♦  ENVC 
■  AENV 


1    3    5    7    9   11  13  15  17  19  21  23 
Hour  of  day  (Sep.) 
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October 


„  0.8 
i  0.6 


♦  ^4 


0.4 

0.2 
0 


♦ 


♦  ENVC 
■  AENV 


3    5    7    9   11  13  15  17  19  21  23 
Hour  of  day  (Oct.) 


November 


^  0.8 
1  0.6 


0.4 
0.2 
0 


♦  ♦♦ 


,  t iirt Uitl ■  ■  ■  ■  ■  ■  ■  Jtii t r f 


3    5    7    9   11  13  15  17  19  21  23 
Hour  of  day  (Nov.) 
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^  0.8 
w 
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0.2 
0 
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1    3    5    7    9   11  13  15  17  19  21  23 
Hour  of  day  (Dec.) 
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APPENDIX  F      COUNTS  FOR  RESISTANCE  FACTORS 


Table  Fl        Number  of  times  hourly  calculated  Ra,  Rb,  or  Re  was  the  largest 
(dominant)  resistance  factor  in  AENV  model  for  SO2. 


J  an 


feb 


apr 


may 


jun 


jul 


aug 


sep 


oct 


dec 


value  of  total  R  (s/m) 


dep  vel  Vd  (cm/s) 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


<10 


>10.0 


10-100 


.0-10.0 


100-300 


0.3-1.0 


300-1000 


0.1-0.3 


346 
0 
0 

346 


351 


351 


421 

0 
0 
421 


520 
4 
0 

516 


479 
5 
0 

474 


475 
22 
0 
453 


457 
19 
0 

438 


461 
18 
0 
443 


414 
27 
0 
387 


353 
25 
0 
328 


483 
0 
0 

483 


525 
0 
0 

525 


>1()00 


<0.1 
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Table  F2        Number  of  times  hourly  calculated  Ra,  Rb,  or  Rc  was  the  largest 
(domdnant)  resistance  factor  in  ENVC  model  for  SO2. 


value  of  total  R  (s/m) 

<10 

10-100 

1 00-300 

300- 1 000 

>1000 

dep  vel  Vd  (cm/s) 

>10.0 

1.0-10.0 

0.3-1.0 

0.1-0.3 

<0.1 

jsn 

total  resistance 

Q 

3 

68 

300 

372 

Ra 

n 

w 

0 

0 

0 

0 

Rb 

0 

0 

0 

0 

0 

Rc 

Q 

3 

68 

300 

372 

feb 

total  resistance 

u 

0 

8 

339 

325 

Ra 

0 

0 

0 

0 

0 

Rb 

Q 

0 

0 

0 

0 

Rc 

n 

KJ 

0 

8 

339 

325 

mar 

total  resistance 

Q 

8 

157 

359 

220 

Ra 

0 

0 

0 

0 

Rb 

Q 

0 

0 

0 

0 

Rc 

8 

157 

359 

220 

apr 

total  resistance 

0 

10 

301 

312 

97 

Ra 

0 

0 

0 

0 

0 

Rb 

0 

0 

0 

0 

0 

Rc 

0 

10 

301 

312 

97 

may 

total  resistance 

0 
u 

30 

311 

271 

132 

Ra 

0 

0 

0 

0 

0 

Rb 

0 

0 

0 

0 

0 

Rc 

0 

30 

31 1 

271 

132 

jun 

total  resistance 

Q 

19 

299 

306 

96 

Ra 

fl 

VJ 

0 

0 

0 

0 

Rb 

0 

0 

0 

0 

0 

Rc 

w 

19 

299 

306 

96 

iul 

total  resistance 

Q 

0 

275 

356 

1 13 

Ra 

0 

0 

0 

0 

Rb 

0 

0 

0 

0 

0 

Rc 

0 

0 

275 

356 

1 13 

total  resistance 

Q 

5 

258 

342 

139 

Ra 

0 

0 

0 

0 

Rb 

0 

0 

0 

0 

0 

Rc 

Q 

5 

258 

342 

139 

Sep 

total  resistance 

Q 

22 

320 

270 

108 

Ra 

0 

0 

0 

0 

0 

Rb 

Q 

0 

0 

0 

0 

Rc 

0 

22 

320 

270 

108 

Oct 

total  resistance 

u 

12 

317 

264 

151 

Ra 

0 

0 

0 

0 

0 

Rb 

0 

0 

0 

0 

0 

Rc 

0 

12 

317 

264 

151 

nov 

total  resistance 

0 

1 

117 

355 

247 

Ra 

0 

0 

0 

0 

0 

Rb 

0 

0 

0 

0 

0 

Rc 

0 

1 

117 

355 

247 

dec 

total  resistance 

0 

1 

89 

385 

269 

Ra 

0 

0 

0 

0 

0 

Rb 

0 

0 

0 

0 

0 

Rc 

0 

1 

89 

385 

269 
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Table  F3        Number  of  times  hourly  calculated  Ra,  Rb,  or  Rc  was  the  largest 
(dominant)  resistance  factor  in  AENV  model  for  NO2. 


jan 


feb 


apr 


may 


jun 


jul 


aug 


sep 


dec 


value  of  total  R  (s/m) 


dep  vel  Vd  (cm/s) 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


total  resistance 

Ra 

Rb 

Rc 


<10 


>10.0 


10-100 


1.0-10.0 


100-300 


0.3-1.0 


159 
0 
0 

159 


300 
0 
0 

300 


397 
0 
0 

397 


403 
0 
0 

403 


395 
0 
0 

395 


300-1000 


0.1-0.3 


213 
13 
0 

200 


136 
16 
0 
120 


390 
7 
0 

383 


271 
6 
0 

265 
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Table  F4        Number  of  times  hourly  calculated  Ra,  Rb,  or  Rc  was  the  largest 
(dominant)  resistance  factor  in  ENVC  model  for  NO2. 


vflliip  of  tnfal  R  (ii/m\ 

<10 

10-100 

100-300 

I  \J\J  ■^\J\J 

300- 1 000 

^\J\J  l\J\J\J 

>1000 

dep  vel  Vd  (cm/s) 

>10.0 

1.0-10.0 

0.3-1.0 

0.1-0.3 

<0.1 

total  resistance 

0 

0 

17 

199 

527 

Ra 

0 

0 

0 

0 

0 

Rb 

0 

0 

0 

0 

0 

Rc 

0 

0 

17 

199 

527 

feb 

total  resistance 

0 

0 

2 

184 

486 

Ra 

0 

0 

0 

0 

0 

Rb 

0 

0 

0 

0 

0 

Rc 

0 

0 

2 

184 

486 

rnsr 

total  resistance 

0 

0 

72 

344 

328 

Ra 

0 

0 

0 

0 

0 

Rb 

0 

0 

0 

0 

0 

Rc 

0 

0 

72 

344 

328 

apr 

total  resistance 

0 

0 

161 

390 

169 

Ra 

0 

0 

0 

0 

0 

Rb 

0 

0 

0 

0 

0 

Rc 

0 

0 

161 

390 

169 

may 

total  resistance 

0 

0 

197 

364 

183 

Ra 

0 

0 

0 

0 

0 

Rb 

0 

0 

0 

0 

0 

Rc 

0 

0 

197 

364 

183 

inn 

JUU 

total  resistance 

0 

0 

188 

379 

153 

Ra 

0 

0 

0 

0 

0 

Rb 

0 

0 

0 

0 

0 

Rc 

0 

0 

188 

379 

153 

iul 

total  resistance 

0 

0 

97 

463 

184 

Ra 

0 

0 

0 

0 

0 

Rb 

0 

0 

0 

0 

0 

Rc 

0 

0 

97 

463 

184 

dug 

total  resistance 

0 

0 

92 

434 

218 

Ra 

0 

0 

0 

0 

0 

Rb 

0 

0 

0 

0 

0 

Rc 

0 

0 

92 

434 

218 

Sep 

total  resistance 

0 

0 

170 

383 

167 

Ra 

0 

0 

0 

0 

0 

Rb 

0 

0 

0 

0 

0 

Rc 

0 

0 

170 

383 

167 

Oct 

total  rpsititanrp 

0 

0 

144 

398 

202 

Ra 

0 

0 

0 

0 

0 

Rb 

0 

0 

0 

0 

0 

Rc 

0 

0 

144 

398 

202 

nov 

total  resistance 

0 

0 

29 

302 

389 

Ra 

0 

0 

0 

0 

0 

Rb 

0 

0 

0 

0 

0 

Rc 

0 

0 

29 

302 

389 

dec 

total  resistance 

0 

0 

15 

277 

452 

Ra 

0 

0 

0 

0 

0 

Rb 

0 

0 

0 

0 

0 

Rc 

0 

0 

15 

277 

452 
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Table  F5        Number  of  times  hourly  calculated  Ra,  Rb,  or  Rc  was  the  largest 
(dominant)  resistance  factor  in  AENV  model  for  HNO3. 


vnliip  of  total  R  <'<;/iTi\ 

<10 

10- 

100 

100-300 

300- 

1000 

>l  ()()() 

>10.0 

1 .0- 

10.0 

0.3-1.0 

0. 1 

0.3 

<().  I 

jan 

tntnl  ri^Qi*it iinpp 

0 

257 

73 

41 

372 

Ra 

0 

22 

58 

40 

1  75 

Rb 

0 

210 

15 

1 

0 

Rc 

0 

25 

0 

0 

0 

feb 

total  resistance 

0 

291 

46 

21 

314 

Ra 

0 

8 

29 

21 

98 

Rb 

0 

261 

17 

0 

0 

Rc 

0 

22 

0 

0 

0 

total  resistance 

0 

383 

29 

24 

308 

Ra 

0 

16 

26 

24 

1 13 

Rb 

0 

267 

3 

0 

0 

Rc 

0 

100 

0 

0 

Q 

apr 

total  resistance 

0 

469 

44 

16 

191 

Ra 

0 

35 

36 

16 

87 

Rb 

0 

351 

8 

0 

0 

Rc 

0 

83 

0 

0 

0 

may 

tntfll  rpsiQtanrp 

0 

438 

44 

30 

232 

Ra 

0 

48 

43 

30 

152 

Rb 

0 

21 1 

1 

0 

0 

Rc 

0 

179 

0 

0 

0 

jun 

total  resistance 

0 

461 

54 

37 

168 

Ra 

0 

49 

54 

37 

149 

Rb 

0 

319 

0 

0 

0 

Rc 

0 

93 

0 

0 

0 

iiil 

tntfil  rPQiQtj^nr*p 

0 

463 

59 

35 

187 

Ra 

0 

7 

35 

35 

122 

Rb 

0 

377 

24 

0 

0 

Rc 

0 

79 

0 

0 

0 

aug 

total  resistance 

0 

447 

63 

27 

207 

Ra 

0 

22 

51 

27 

158 

Rb 

0 

374 

12 

0 

0 

Rc 

0 

51 

0 

0 

0 

Sep 

total  resistance 

0 

447 

50 

36 

187 

Ra 

0 

18 

46 

36 

139 

Rb 

0 

284 

4 

0 

0 

Rc 

0 

145 

0 

Q 

Q 

Oct 

total  rp^ii^tanpp 

0 

377 

77 

5 1 

239 

Ra 

0 

61 

74 

51 

219 

Rb 

0 

252 

3 

0 

0 

Rc 

0 

64 

0 

0 

0 

nov 

total  resistance 

0 

373 

94 

33 

220 

Ra 

0 

57 

87 

33 

178 

Rb 

0 

305 

7 

0 

0 

Rc 

0 

11 

0 

0 

0 

dec 

total  resistance 

0 

367 

139 

47 

191 

Ra 

0 

360 

139 

47 

141 

Rb 

0 

0 

0 

0 

0 

Rc 

0 

7 

0 

0 

0 
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Table  F6        Number  of  times  hourly  calculated  Ra,  Rb,  or  Rc  was  the  largest 
(dominant)  resistance  factor  in  ENVC  model  for  HNO3. 


VdlUC  Ul  lUlal  Ix  ^^Vlll^ 

<10 

10  inn 

ion  "^00 

•^00  1000 

J WW- 1 www 

^1  nno 
^1  www 

UCp  VCl  V  U 

>10.0 

1  0-100 

0  3-10 

0  1-03 

^0  1 

LUldl  ICSlslallLC 

0 
u 

9QQ 

1 1 

Ql 

y  I 

Ra 

u 

12 

45 

53 

68 
uo 

Rb 

3 

0 

A 
W 

n 

w 

u 

222 

ZJt 

yjj 

9^ 

feb 

lUlal  IColslullCC 

A 
U 

205 

336 

67 

f\A 
ut 

Ra 

u 

3 

4 

jy 

Al 
tz 

Rb 

n 

0 

0 

0 

0 

Rc 

202 

332 

98 

22 

ivjlal  iCololallCC 

u 

400 

198 

80 

66 

Ra 

0 
u 

4 

16 

51 

^7 

J  / 

Rh 

J  1 

A 

w 

A 

w 

A 
W 

Rp 

U 

1  89 

9Q 

9Q 
Ay 

apr 

total  rpQwtanpp 

478 

161 

56 

25 

Rn 

u 

99 

J  J 

so 

JW 

9S 

Rb 

n 
u 

Jt 

Q 

A 
U 

A 
W 

Rr 

u 

409 

1  96 

A 
0 

A 
W 

may 

lUldl  ICSlcildUcC 

n 

W 

462 

154 

80 

48 

Ra 

43 

UJ 

ot 

Al 

t  / 

Rb 

u 

56 

0 

0 

Rc 

n 

363 

89 

16 

n 

w 

jun 

tot^il  rpcictanpp 
lUldl  ICdJ^ldllCC 

0 
U 

433 

202 

69 

16 

Ra 

U 

70 

lOQ 

ut 

Rb 

VJ 

16 

0 
w 

A 
W 

A 
W 

Rc 

n 
U 

347 

93 

5 

0 

ill! 

total  rpciQtanpp 

n 
U 

431 

197 

101 

15 

Rq 
Ixd 

n 
U 

Jt 

8Q 

101 

IWl 

1  < 
1 J 

Rb 

U 

97 

n 

A 
W 

A 
W 

Rr 

n 

U 

108 
IWo 

A 
W 

A 
W 

aug 

total  rpcictanr^p 

A 

U 

4^7 

160 

135 

12 

Ra 
Ixd 

A 

U 

to 

87 
0  / 

LJJ 

19 
iz 

Rb 

A 

19 

0 

0 

0 

Rc 

A 

u 

^70 

1  J 

n 

w 

n 

w 

sep 

total  rpcictanpp 

A 
V 

446 

134 

129 

1 1 

Ra 

A 

62 

103 

129 

1  [ 

Rb 

A 
U 

43 

0 

0 

0 

Rc 

A 
U 

341 

3 1 

0 

Oct 

total  fpcictanr*p 
lUldl  ICalSldllCC 

A 
U 

"^80 

JOK) 

IQQ 

1 1 7 

48 
to 

Ra 

0 

152 

163 

1 14 

48 

Rb 

0 

19 

0 

0 

0 

Rc 

0 

209 

36 

3 

0 

nov 

total  resistance 

0 

334 

229 

88 

69 

Ra 

0 

13 

62 

80 

69 

Rb 

0 

6 

0 

0 

0 

Rc 

0 

315 

167 

8 

0 

dec 

total  resistance 

0 

339 

288 

75 

42 

Ra 

0 

7 

42 

69 

42 

Rb 

0 

9 

0 

0 

0 

Rc 

0 

323 

246 

6 

0 
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Table  F7        Number  of  times  hourly  calculated  Ra,  Rb,  or  Rc  was  the  largest 

(dominant)  resistance  factor  in  AENV  model  for  HNO3  for  total  resistance 
<100  s/m. 


value  of  total  R  (s/m) 

<15 

15-25 

25-40 

40-65 

65- 100 

dep  vel  Vd  (cm/s) 

>6.7 

4.0-6.7 

2.5-4.0 

1.5-2.5 

1 .0- 1 .5 

jan 

total  resistance 

0 

19 

67 

103 

68 

Ra 

0 

0 

0 

3 

19 

Rb 

0 

0 

61 

100 

49 

Rc 

0 

19 

6 

0 

0 

feb 

total  resistance 

0 

7 

1 1 1 

121 

52 

Ra 

0 

0 

0 

1 

7 

Rb 

0 

0 

96 

120 

45 

Rc 

0 

7 

15 

0 

0 

mar 

total  resistance 

0 

74 

173 

103 

33 

Ra 

0 

0 

0 

3 

13 

Rb 

0 

0 

147 

100 

20 

Rc 

0 

74 

26 

0 

0 

apr 

total  resistance 

0 

29 

197 

162 

81 

Ra 

0 

0 

0 

3 

32 

Rb 

0 

0 

143 

159 

49 

Rc 

0 

29 

54 

0 

0 

may 

total  resistance 

0 

146 

189 

69 

34 

Ra 

0 

0 

0 

17 

31 

Rb 

0 

0 

156 

52 

3 

Rc 

0 

146 

33 

0 

0 

jun 

total  resistance 

0 

58 

242 

123 

38 

Ra 

0 

0 

4 

17 

28 

Rb 

0 

0 

203 

106 

10 

Rc 

0 

58 

35 

0 

0 

Jul 

total  resistance 

0 

51 

183 

171 

58 

Ra 

0 

0 

0 

0 

7 

Rb 

0 

0 

155 

171 

51 

Rc 

0 

51 

28 

0 

0 

aug 

total  resistance 

0 

31 

197 

168 

51 

Ra 

0 

0 

0 

2 

20 

Rb 

0 

0 

177 

166 

31 

Rc 

0 

31 

20 

0 

0 

sep 

total  resistance 

0 

124 

179 

105 

39 

Ra 

0 

0 

0 

1 

17 

Rb 

0 

0 

158 

104 

22 

Rc 

0 

124 

21 

0 

0 

oct 

total  resistance 

0 

42 

155 

1 12 

68 

Ra 

0 

0 

0 

12 

49 

KD 

0 

n 
U 

1 1 1 

Rc 

0 

42 

22 

0 

0 

nov 

total  resistance 

0 

1 

121 

164 

87 

Ra 

0 

0 

0 

15 

42 

Rb 

0 

0 

111 

149 

45 

Rc 

0 

1 

10 

0 

0 

dec 

total  resistance 

0 

3 

25 

177 

162 

Ra 

0 

0 

21 

177 

162 

Rb 

0 

0 

0 

0 

0 

Rc 

0 

3 

4 

0 

0 
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Table  F8 


Number  of  times  hourly  calculated  Ra,  Rb,  or  Rc  was  the  largest 
(dominant)  resistance  factor  in  ENVC  model  for  HNO3  for  total  resistance 
<100  s/m. 


value  of  total  R  (s/m) 

<15 

15-25 

25-40 

40-65 

65-100 

dep  vel  Vd  (cm/s) 

>6.7 

4.0-6.7 

2.5-4.0 

1.5-2.5 

1.0-1.5 

jan 

total  resistance 

0 

11 

19 

73 

134 

Ra 

0 

0 

0 

4 

8 

Rb 

0 

0 

0 

2 

1 

Rc 

0 

11 

19 

67 

125 

feb 

total  resistance 

0 

0 

1 

75 

129 

Ra 

0 

0 

0 

2 

1 

Rb 

0 

0 

0 

0 

0 

Rc 

0 

0 

1 

73 

128 

mar 

total  resistance 

0 

9 

75 

155 

161 

Ra 

0 

0 

0 

0 

4 

Rb 

0 

0 

15 

18 

4 

Rc 

0 

9 

60 

137 

153 

apr 

total  resistance 

0 

15 

138 

170 

155 

Ra 

0 

0 

0 

5 

17 

Rb 

0 

0 

34 

15 

5 

Rc 

0 

15 

104 

150 

133 

may 

total  resistance 

0 

27 

138 

166 

131 

Ra 

0 

0 

5 

17 

21 

Rb 

0 

0 

38 

13 

5 

Rc 

0 

27 

95 

136 

105 

jun 

total  resistance 

0 

11 

135 

150 

137 

Ra 

0 

0 

18 

22 

30 

Rb 

0 

0 

12 

4 

0 

Rc 

0 

11 

105 

124 

107 

Jul 

total  resistance 

0 

9 

78 

186 

158 

Ra 

0 

0 

1 

14 

39 

Rb 

0 

0 

6 

11 

10 

Rc 

0 

9 

71 

161 

109 

aug 

total  resistance 

0 

4 

81 

199 

153 

Ra 

0 

0 

0 

17 

31 

Rb 

0 

0 

9 

5 

5 

Rc 

0 

4 

72 

177 

117 

sep 

total  resistance 

0 

35 

132 

162 

117 

Ra 

0 

0 

1 

23 

38 

Rb 

0 

0 

15 

14 

14 

Rc 

0 

35 

116 

125 

65 

oct 

total  resistance 

0 

9 

114 

150 

107 

Ra 

0 

0 

17 

76 

59 

Rb 

0 

0 

18 

0 

Rc 

0 

9 

79 

73 

48 

nov 

total  resistance 

0 

0 

58 

145 

131 

Ra 

0 

0 

0 

7 

6 

Rb 

0 

0 

5 

1 

0 

Rc 

0 

0 

53 

137 

125 

dec 

total  resistance 

0 

1 

42 

131 

165 

Ra 

0 

0 

0 

5 

2 

Rb 

0 

0 

6 

3 

0 

Rc 

0 

1 

36 

123 

163 
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APPENDIX  G  BREAKDOWN  OF  CHEMICALS  CONTRIBUTING 
TO  DIFFERENCES  IN  DEPOSITION  CAUSED  BY 
PERCENT  CHANGES  OF  RESISTANCES  (Ra,  Rb 
AND  Rc) 

Light  yellow,  gold,  and  orange  highlighting  in  tables  shown  below  indicate  small, 
medium,  and  large  differences  in  deposition  (in      equivalents),  respectively. 

Table  Gl        Lifluence  of  ±10%  and  ±50%  change  in  Ra  on  total  deposition  (in  H"^ 
equivalents)  as  a  result  of  sensitivity  testing  of  AENV  model. 


base  case  total 

absolute  difference  from  base  case  (| 

y  H-H/ha) 

Ra 

0.5 

0.9 

1.1 

1.5 

Month  1 

16.9 

2.6 

0.4 

-0.4 

-1.7 

Month  2 

11.7 

1.0 

0.2 

-0.2 

-0.7 

Month  3 

17.0 

1.3 

0.2 

-0.2 

-0.9 

Month  4 

15.9 

1.7 

0.3 

-0.3 

-1.2 

Month  5 

14.0 

1.5 

0.2 

-0.2 

-1.0 

Month  6 

7.5 

0.7 

0.1 

-0.1 

-0.5 

Month  7 

7.6 

0.3 

0.0 

0.0 

-0.2 

Month  8 

13.6 

1.1 

0.2 

-0.2 

-0.7 

Month  9 

9.8 

0.7 

0.1 

-0.1 

-0.5 

Month  10 

16.8 

2.4 

0.4 

-0.4 

-1.6 

Month  1 1 

0.2 

-0.1 

0.0 

0.0 

0.1 

Month  12 

16.3 

2.3 

0.4 

-0.3 

-1.5 

Refinement  Study  of  Dry  Deposition  Inference 
Methods  Used  in  Alberta 


137 


Table  G2       Influence  of  ±10%  and  ±50%  change  in  Ra  on  HNO2  deposition  (in  H"" 
equivalents)  as  a  result  of  sensitivity  testing  of  AENV  model. 


base  case 

absolute  difference  from  base  case  {\ 

I  H+/ha) 

month 

(g  H+/ha) 

0.5 

0.9 

1.1 

1.5 

Month  1 

3.2 

0.8 

0.1 

-0.1 

-0.5 

Month  2 

0.8 

0.2 

0.0 

0.0 

-0.1 

Month  3 

0.6 

0.1 

0.0 

0.0 

-0.1 

Month  4 

1.0 

0.3 

0.0 

0.0 

-0.2 

Month  5 

1.7 

0.3 

0.1 

0.0 

-0.2 

Month  6 

0.2 

0.0 

0.0 

0.0 

0.0 

Month  7 

0.3 

0.1 

0.0 

0.0 

0.0 

Month  8 

0.5 

0.1 

0.0 

0.0 

-0.1 

Month  9 

1.1 

0.2 

0.0 

0.0 

-0.1 

Month  10 

2.5 

0.5 

0.1 

-0.1 

-0.4 

Month  1 1 

1.2 

0.3 

0.1 

0.0 

-0.2 

Month  12 

4.8 

1.5 

0.2 

-0.2 

-0.9 

Table  G3        Influence  of  ±10%  and  ±50%  change  in  Ra  on  HNO3  deposition  (in  H"" 
equivalents)  as  a  result  of  sensitivity  testing  of  AENV  model. 


base  case 

absolute  difference  from  base  case  (^ 

I  H+/ha) 

month 

(g  H+/ha) 

0.5 

0.9 

1.1 

1.5 

Month  1 

4.8 

1.1 

0.2 

-0.2 

-0.7 

Month  2 

2.5 

0.5 

0.1 

-0.1 

-0.3 

Month  3 

4.1 

0.7 

0.1 

-0.1 

-0.5 

Month  4 

4.6 

1.0 

0.2 

-0.2 

-0.7 

Month  5 

3.8 

0.6 

0.1 

-0.1 

-0.4 

Month  6 

1.7 

0.3 

0.1 

0.0 

-0.2 

Month  7 

1.9 

0.4 

0.1 

-0.1 

-0.3 

Month  8 

2.0 

0.4 

0.1 

-0.1 

-0.3 

Month  9 

1.4 

0.2 

0.0 

0.0 

-0.2 

Month  10 

8.3 

1.5 

0.3 

-0.2 

-1.1 

Month  1 1 

1.0 

0.2 

0.0 

0.0 

-0.2 

Month  12 

1.5 

0.4 

0.1 

-0.1 

-0.3 

Refinement  Study  of  Dry  Deposition  Inference 
Methods  Used  in  Alberta 


138 


Table  G4        Influence  of  ±10%  and  ±50%  change  in  Rb  on  total  deposition  (in  H"^ 
equivalents)  as  a  result  of  sensitivity  testing  of  AENV  model. 


base  case  total 

absolute  difference  from  base  case  (g  H+/ha) 

Kb 

0.5 

u.y 

1 . 1 

1 .5 

Montn  1 

iD.y 

3.Z 

U.3 

-0.4 

- 1 .9 

Montn  L 

W.I 

1  o 

1 .8 

0.3 

-O.z 

- 1 .0 

Montn  3 

1  '7  n 
1  /.U 

3.7 

0.5 

-0.4 

1  o 

- 1 .8 

Month  4 

1  <  n 
1  j.y 

U.  i 

(\  1 
U.  1 

n  1 
-U.  1 

-U.u 

Montn  J 

14. U 

0.8 

0.2 

-0.2 

-0.8 

Month  6 

7.5 

-0.7 

0.0 

0.0 

0.0 

Month  7 

7.6 

-3.4 

-0.4 

0.3 

1.0 

Month  8 

13.6 

-0.6 

0.0 

0.0 

-0.2 

Month  9 

9.8 

0.6 

0.1 

-0.1 

-0.5 

Month  10 

16.8 

0.5 

0.2 

-0.2 

-1.1 

Month  1 1 

-5.8 

-0.7 

0.6 

2.2 

Month  12 

16.3 

0.6 

0.2 

-0.2 

-0.7 

Table  G5       Influence  of  ±10%  and  ±50%  change  in  Rb  on  HNO3  deposition  (in 
equivalents)  as  a  result  of  sensitivity  testing  of  AENV  model. 


base  case 

absolute  difference  from  base  case  (g  H+/ha) 

month 

(g  H-H/ha) 

0.5 

0.9 

1.1 

1.5 

Month  1 

4.8 

1.3 

0.2 

-0.2 

-0.8 

Month  2 

2.5 

0.7 

0.1 

-0.1 

-0.4 

Month  3 

4.1 

1.0 

0.2 

-0.2 

-0.7 

Month  4 

4.6 

1.1 

0.2 

-0.2 

-0.7 

Month  5 

3.8 

0.8 

0.1 

-0.1 

-0.6 

Month  6 

1.7 

0.4 

0.1 

-0.1 

-0.3 

Month  7 

1.9 

0.5 

0.1 

-0.1 

-0.3 

Month  8 

2.0 

0.5 

0.1 

-0.1 

-0.3 

Month  9 

1.4 

0.3 

0.1 

-0.1 

-0.2 

Month  10 

8.3 

1.9 

0.3 

-0.3 

-1.3 

Month  1 1 

1.0 

0.2 

0.0 

0.0 

-0.2 

Month  12 

1.5 

0.4 

0.1 

-0.1 

-0.2 
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Table  G6       Influence  of  ±10%  and  +50%  change  in  Rb  on  SO/  deposition  (in 
equivalents)  as  a  result  of  sensitivity  testing  of  AENV  model. 


base  case 

absolute  difference  from  base  case  (g  H+/ha) 

month 

(g  H+/ha) 

0.5 

0.9 

1.1 

1.5 

Montn  1 

l.U 

0.7 

0.1 

-0.1 

-0.3 

Montn  z 

l.U 

0.8 

0.1 

f\  1 

-0.1 

-0.3 

Montn  3 

2.4 

0.3 

-0.2 

-0.9 

Montn  4 

1.3 

1  1 
1.1 

n  1 
U.  1 

-U.  1 

-U.4 

jvionin  J 

0.8 

0.1 

-0.1 

-0.3 

Month  6 

1.0 

0.8 

0.1 

-0.1 

-0.3 

Month  7 

1.0 

0.8 

0.1 

-0.1 

-0.3 

Month  8 

1.3 

1.0 

0.1 

-0.1 

-0.4 

Month  9 

0.9 

0.8 

0.1 

-0.1 

-0.3 

Month  10 

1.3 

1.0 

0.1 

-0.1 

-0.4 

Month  1 1 

1.2 

0.9 

0.1 

-0.1 

-0.3 

Month  12 

1.0 

0.7 

0.1 

-0.1 

-0.3 

Table  G7        Influence  of  ±10%  and  ±50%  change  in  Rb  on  NH/  deposition  (in 
equivalents)  as  a  result  of  sensitivity  testing  of  AENV  model. 


base  case 

absolute  difference  from  base  case  (] 

I  H-h/ha) 

month 

(g  H+/ha) 

0.5 

0.9 

1.1 

1.5 

Month  1 

1.7 

1.1 

0.1 

-0.1 

-0.4 

Month  2 

1.3 

1.0 

0.1 

-0.1 

-0.4 

Month  3 

2.2 

1.8 

0.2 

-0.2 

-0.7 

Month  4 

1.8 

1.5 

0.2 

-0.1 

-0.5 

Month  5 

1.2 

1.0 

0.1 

-0.1 

-0.4 

Month  6 

1.2 

1.0 

0.1 

-0.1 

-0.4 

Month  7 

1.2 

0.9 

0.1 

-0.1 

-0.3 

Month  8 

1.7 

1.3 

0.2 

-0.1 

-0.5 

Month  9 

1.2 

1.0 

0.1 

-0.1 

-0.4 

Month  10 

1.5 

1.2 

0.1 

-0.1 

-0.4 

Month  1 1 

2.7 

2.0 

0.2 

-0.2 

-0.8 

Month  12 

1.2 

0.9 

0.1 

-0.1 

-0.3 
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Table  G8        Influence  of  ±10%  and  ±50%  change  in  Rb  on  Ca"^  deposition  (in 
equivalents)  as  a  result  of  sensitivity  testing  of  AENV  model. 


base  case 

absolute  difference  from  base  case  (g  H-f-/ha) 

month 

(g  H+/ha) 

0.5 

0.9 

1  1 
1 . 1 

1  .J 

Montn  1 

- 1 .3 

-0.5 

C\  1 

-0. 1 

0. 1 

A  O 

O.Z 

Montn  z 

-z.U 

-0.9 

-0. 1 

0. 1 

A  Q 

0.3 

Montn  3 

-3.3 

1  A 

-1.4 

-O.z 

A  1 

0. 1 

A  C 

0.5 

Montn  4 

VHHp 

-3.  / 

A  /I 

-0.4 

A  A 

U.4 

1  A 

1 .4 

Monin  J 

A  1 

-1.9 

-0.2 

0.2 

0.7 

Month  6 

-6.6 

-2.8 

-0.3 

0.3 

1.1 

Month  7 

"^imiP^  -12.0 

-5.1 

-0.6 

0.5 

2.0 

Month  8 

-8.5 

-3.6 

-0.4 

0.4 

1.4 

Month  9 

-3.1 

-1.4 

-0.2 

0.1 

0.5 

Month  10 

-9.0 

-3.9 

-0.5 

0.4 

1.5 

Month  1 1 

-22.0 

-9.0 

-1.1 

1.0 

3.7 

Month  12 

-5.5 

-2.2 

-0.3 

0.2 

0.9 

Table  G9       Influence  of  ±10%  and  ±50%  change  in  Rc  on  total  deposition  (in 
equivalents)  as  a  result  of  sensitivity  testing  of  AENV  model. 


base  case  total 

absolute  difference  from  base  case  {\ 

I  H+/ha) 

Rc 

0.5 

0.9 

1.1 

1.5 

Month  1 

16.9 

6.7 

0.9 

-0.7 

-2.9 

Month  2 

11.7 

6.5 

0.8 

-0.7 

-2.6 

Month  3 

17.0 

8.6 

1.1 

-0.9 

-3.5 

Month  4 

15.9 

10.6 

1.3 

-1.1 

-4.4 

Month  5 

14.0 

8.5 

1.1 

-0.9 

-3.6 

Month  6 

7.6 

5.8 

0.7 

-0.6 

-2.4 

Month  7 

7.6 

9.4 

1.1 

-1.0 

-3.6 

Month  8 

13.6 

10.4 

1.3 

-1.1 

-4.2 

Month  9 

9.8 

6.6 

0.8 

-0.7 

-2.6 

Month  10 

16.8 

8.9 

1.2 

-1.0 

-4.1 

Month  1 1 

0.2 

6.2 

0.8 

-0.6 

-2.4 

Month  12 

16.3 

9.8 

1.2 

-1.0 

-4.0 
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Table  GIO      Influence  of  ±10%  and  ±50%  change  in  Rc  on  SO2  deposition  (in 
equivalents)  as  a  result  of  sensitivity  testing  of  AENV  model. 


base  case 

absolute  difference  from  base  case  (g  H+/ha) 

month 

0.5 

0.9 

1.1 

1.5 

Montn  I 

1  0 
1  .z 

0.9 

A  1 

0.1 

A  1 

-0.1 

A  ^ 

-0.3 

Montn  z 

1.0 

1.3 

O.z 

A  1 

-0. 1 

A  C 
-O.J 

iviontn  D 

z.z 

1  0 
1.0 

A  0 

O.Z 

A  0 

-O.Z 

A  ^ 

-0.6 

JVlontn  4 

3. 1 

Z.  / 

A  Q 

A  ^ 

-U.3 

1  A 

-l.U 

ivionin  J 

Z.D 

2.3 

0.3 

-0.2 

-0.8 

Month  6 

3.0 

1.8 

0.2 

-0.2 

-0.6 

Month  7 

4.5 

3.6 

0.4 

-0.3 

-1.3 

Month  8 

5.1 

3.2 

0.4 

-0.3 

-1.2 

Month  9 

2.1 

1.3 

0.2 

-0.1 

-0.5 

Month  10 

2.1 

1.3 

0.2 

-0.1 

-0.5 

Month  1 1 

0.6 

0.5 

0.1 

-0.1 

-0.2 

Month  12 

2.7 

2.1 

0.3 

-0.2 

-0.8 

Table  Gl  1      Influence  of  ±10%  and  ±50%  change  in  Rc  on  NO2  deposition  (in 
equivalents)  as  a  result  of  sensitivity  testing  of  AENV  model. 


base  case 

absolute  difference  from  base  case  (| 

I  H+/ha) 

month 

(g  H+/ha) 

0.5 

0.9 

1.1 

1.5 

Month  1 

6.4 

4.7 

0.6 

-0.5 

-1.8 

Month  2 

6.1 

4.8 

0.6 

-0.5 

-1.8 

Month  3 

7.3 

5.8 

0.7 

-0.6 

-2.2 

Month  4 

9.5 

7.0 

0.9 

-0.7 

-2.7 

Month  5 

6.8 

5.0 

0.6 

-0.5 

-1.9 

Month  6 

5.1 

3.7 

0.5 

-0.4 

-1.5 

Month  7 

7.0 

5.4 

0.7 

-0.5 

-2.1 

Month  8 

9.2 

6.8 

0.8 

-0.7 

-2.7 

Month  9 

5.7 

4.7 

0.6 

-0.5 

-1.7 

Month  10 

7.0 

5.4 

0.7 

-0.5 

-2.0 

Month  1 1 

6.9 

5.3 

0.6 

-0.5 

-2.0 

Month  12 

8.9 

7.1 

0.8 

-0.7 

-2.6 
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APPENDIX  H 


DEPOSITION  CONTRIBUTIONS  FOR  SPECIES 
ASSOCIATED  WITH  LARGE  PARTICLES 
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Figure  HI       Monthly  summed  deposition  flux  for  NO3  for  AENV  model  base  case, 
ENVC  model,  and  ENVC  model  corrected  for  large  particle  size 
distribution  for  a  coniferous  forest  cover. 
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Figure  H2      Monthly  summed  deposition  flux  for  Na^  for  AENV  model  base  case, 
ENVC  model,  and  ENVC  model  corrected  for  large  particle  size 
distribution  for  a  coniferous  forest  cover. 
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Figure  H3       Monthly  summed  deposition  flux  for     for  AENV  model  base  case, 
ENVC  model,  and  ENVC  model  corrected  for  large  particle  size 
distribution  for  a  coniferous  forest  cover. 
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Figure  H4      Monthly  summed  deposition  flux  for  Ca^^  for  AENV  model  base  case, 
ENVC  model,  and  ENVC  model  corrected  for  large  particle  size 
distribution  for  a  coniferous  forest  cover. 
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Figure  H5       Monthly  summed  deposition  flux  for  Mg^^  for  AENV  model  base  case, 
ENVC  model,  and  ENVC  model  corrected  for  large  particle  size 
distribution  for  a  coniferous  forest  cover. 
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